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The Slip, Twinning, Cohesion, Growth and Boundaries 
of. Crystals 


By H. WILMAN 
Imperial College, London 


MS. received 26th Fuly 1949, and in amended form 15th September 1950 


ABSTRACT. The nature and stability of intercrystalline boundaries are discussed and 
conclusions as to preferred relative orientations are illustrated by ball-bearing and bubble- 
raft models. 

Electron-diffraction evidence from crystals of layer-lattice, ionic, metallic and van der 
Waals (organic molecular) types is presented, showing that crystal pairs and symmetrical 
or unsymmetrical triplet or multiplet crystal groupings occur, having a common lattice 
row, but in relative azimuthal orientations at intervals agreeing with the predictions. The 
component crystals were thus in coherent metastable contact of the type discussed, and in 
some cases strong secondary elastic scattering showed that the crystals were extensively 
superposed on the plane which was normal to the common row. 

The seven types of electron-diffraction pattern observed from such groupings, and 
the nature of the specimen preparation, suggest that the origin of these groupings is a 
mechanical deformation of the crystals during or after growth, by a process of ‘rotational 
slip’ which has not hitherto been explored or defined in nature. Such slip is directly 
demonstrated macroscopically as a deformation process in potassium ferrocyanide trihydrate 
and gypsum, the slip being produced on the planes parallel to the highly perfect cleavage. 

The relation between rotational slip and the known types of translational slip and 
dislocation theory is indicated, and in particular the ‘ deformation bands’ investigated 
especially by Barrett and his collaborators are concluded to be essentially rotational slip 
bands initially at low deformation, though modified by translational types of slip at higher 
deformation. - 

A valuable new insight is also provided into many observations on crystal growth and 
properties, especially the nature, determining conditions, stability and disorientation of 
crystals growing epitaxially on single crystal substrates, and into the nature of the deforma- 
tion caused by unidirectional abrasion of crystals. 


$1. INTRODUCTION 
HE following electron-diffraction and other results, already partly 
reported (Wilman 1950), lead to new considerations of the nature and 
origin of certain types of intercrystalline boundaries and to the macroscopic 
demonstration of a hitherto unexplored type of deformation of crystals, best 
called ‘rotational slip’. Such slip can be defined as the slipping of one part 
of a crystal on a neighbouring part so that the two atomic or molecular sheets 
which slide over one another are parallel to densely populated net planes, as in 
translational slip, but rotationally displaced about an axis normal to their plane 
(with or without simultaneous translational displacement) to one of a series of 
PROC. PHYS. SOC. LXIV, 4—A 22 
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definable new azimuths where metastable equilibrium can occur. ‘These results 
provide an interpretation of several types of electron-diffraction patterns, an 
explanation of many observations on the growth and properties of crystals, and 
especially a new and fundamental insight into (i) the nature of epitaxial crystal 
growth (see also Evans and Wilman 1950), (ii) the nature of ‘ deformation bands’ 
in crystals, and (iii) the deformation caused by unidirectional abrasion of a 
crystal. 


§2. THEORETICAL CRYSTAL ORIENTATIONS CORRESPONDING 
TO METASTABLE INTERCRYSTALLINE BOUNDARIES 
The lattices of two crystals of the same species, if not in the same orientation, 
are related by a rotation through an angle 6 about a common lattice direction. 
The present work concerns primarily those boundaries where two crystals meet 
at a plane which is (a) normal to the common lattice direction, or (b) parallel to it, 


(a) Superposition of two crystals in different azimuths on a similar boundary face. 

One of the simplest types of intercrystalline boundary is a plane such that the 
net planes parallel to it in the two crystals are of the same form, though the crystals 
are rotated relative to each other about the common lattice direction normal to 
this plane. In this case a proportion of the atoms (or ions or molecules) in the 
contacting planes lie in potential troughs similar to those they would occupy if 
the two planes were superposed in the same azimuth, as in the single crystals, 
but most of the atoms are raised partly out of such potential depressions 
(Figures 1 and 2). 

As Figure 2 illustrates, the atoms in one sheet which coincide with the 
potential troughs in the neighbouring sheet lie on a two-dimensional network. 
If the atomic lattice is rectangular with axes a,b, and if the two superposed 
sheets are rotated in opposite directions by 5/2 from the parallel azimuth so that 
the [u, v,] and [w, v,] rows coincide along the initial a axis direction (Figure 3 (a)), 
then 

tan6/2=v,b/u,a;  Dy=(u;2a"4-0,7b2)2 =0,6 coseco/Z., i .cnice (1) 
At the same time two rows [u.v,] and [u, v2] will coincide along the initial 5 
direction, thus 


tan0/Z=1,0/0,b3" Dy=(v9*b* usa") 2 =u, acoseco/ ze akan (2) 
The axial ratio of the coincidence lattice is therefore 
D2/D =(ta[0,) (a/b); eee (3) 


thus if v;=w,=1 this ratio becomes a/b, i.e. the reciprocal of that of the atomic 
lattice, so that the coincidence lattice is similar to the atomic lattice in each sheet 
but enlarged by a factor cosec 8/2 and rotated by 90° from the initial azimuth 
of superposition. In other cases the axes D, and D, are further multiplied by 
V, and uw, respectively. The condition for the two types of coincidences to 
occur simultaneously is that u, w/v, V2 =b?/a?. 

It will be seen that the two lattices are reflection twins of each other across 
their (v#0) or (vuw0) planes respectively, though the crystals themselves will only 
be twins if their symmetry of atomic arrangement is appropriately high. 

Since the number of coincidences per unit area is a measure of the relative 
stability (potential energy of fitting together of the two sheets) it may be concluded 
that metastable equilibrium will be possible with a high order of stability at 
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i discrete relative azimuths 5 corresponding to small (integral) values of 


| Uy, V1, Uz and vy. ‘The potential energy of fitting together of the two crystals is 


of course partly determined by the local atomic distributions in addition to the 


»§ factor expressing the effect of the relative lattice periodicities. 


{7 Figure 1. Fitting between an atom row forming part of an upper layer, on a similar lower layer 


which is in an azimuthally rotated position in which many atoms no longer lie at potential 
minima. 


Figure 2. A hexagonal sheet of atoms (++) rotated through 6=9° 26’ over the two sets of potential 
troughs (@ and ©) of a similar atomic sheet. 


The above expressions apply to plane lattices of square, simple rectangular, 
rhomboid (centred rectangular) and hexagonal (centred rectangular with b/a = 1/3) 
type. In the remaining most general lattice type, with unequal axes not at 90°, 
the azimuths at which relatively stable atomic lattice superposition can occur are 

22-2 
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less easily defined. If any two rows in the plane lattice are at 90° to each other or 
very close to it, the above expressions can be applied. If not, then even if exact 
coincidences are rare there may for certain values of 5 be rows of approximate 
coincidences moderately closely distributed in the interface. At small angles 6 
such approximate coincidences will occur (cf. Figure 3(c)) at a spacing 
(n,ab sin y)/T, in a direction through the rotation axis and normal to it in the 
plane of contact, along the normal to a lattice translation T,, when 


6 =n, T,/{(n,ab sin y)/T,}=(ny/ne)T?/absiny,  —...... (4) 
where n, and ng are integers, small and large respectively. If the axes a or b are 


the smallest lattice translations, the distance between coincidence points is lowest 
when 7, is taken to be a or J, so that, respectively, 


d—n,b/(n,asiny), or d—n,a/(n.bsiny). eer 
Green and Weigle (1948) have stated somewhat similar expressions for the mozre 


pattern of two superposed sets of parallel lines, differing slightly from the above 
results applicable to lattices. 


(b) Contact of two crystals on a plane parallel to the common lattice row about which 
they are rotationally displaced relative to each other. 

Provided the common axis c is a moderately densely populated row with none 
more densely populated at a small angle to it, the azimuths 6 defined by 
equations (1) and (2) are also characteristic of metastable boundaries between 
crystals in this type of contact, which may be called ‘laterally adjacent’ instead 
of superposed. At any one value of 6 the plane of contact could be that through ¢ 
and the bisector of the two a or two b axes normal to ¢ in the two crystals, or any 
other plane (or set of plane elements) parallel to c, the stability of the boundary 
being proportional to the density of coincidence points per unit area of boundary 
plane, and thus dependent on the direction of the contact plane or plane elements 
of contact. Superposed and adjacent contact faces can clearly occur 
simultaneously between crystals in which the axes satisfy the above conditions. 


(c) Illustrations of the above types of boundary by models. 

The discrete azimuths 5 for metastable ‘superposed’ contact are well shown 
by a ball-bearing model. If a horizontal layer of Zin. or 4in. diameter balls, 
hexagonally close packed in a hexagon of about 10-15 balls per side, is held as 
a unit by a close-fitting ‘lid’ with turned-down edges and is rotated over a larger 
similar layer, it is found that although much the deepest collective potential 
trough (due to gravity) occurs at zero azimuth difference, a large number of stable 
positions of rest are possible in shallower potential minima at the values of 6 
defined by (1) and (2), that at 6 =27° 48’ (or 32° 12’) being especially deep. When 
the upper layer was reduced to only five balls per side these minima were scarcely 
noticeable. When 6 =0 translational slip was strongly limited to the zig-zags of 
4 (21) type, making up a net (10) glide (cf. Seitz and Read 1941, Barrett 1943, 
Zener 1948); but when 6 was even only a degree or two (using the larger layers) 
this was replaced by practically equally easy translational slip in any direction, 
with scarcely a trace of uneven passage over collective potential maxima and 
minima except at a few azimuths such as 6=27°48’ where the regions of 
coincidence were relatively closely spaced. 

Bubble-raft models (Bragg 1942, Bragg and Nye 1947, Bragg and Lomer 
1949) provide not only examples of ‘ dislocations’ and a means of studying the 
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ite i effect of these on electron-diffraction patterns (Dyson 1949) but also, as Figure 4 
‘ol shows, illustrate well the preferred discrete azimuths at which the more stable 
ings] “laterally adjacent’ intercrystal boundaries occur, and the way in which these 
(so are compatible with curved or irregular contact surfaces. In Figure 4 (ii), for 
‘nee example, 6 =22° (or 38° or 82°), and owing to the hexagonal symmetry there are 
\/i@) coincident orthohexagonal (21) lattice translations in a collinear series IJ, JK, 
KL, LM and also CD of (51) or (13) type, RS of (21) and OQ of (10:2) type. 
( y In Figure 4 (iii) 8=28° (or 32° etc.) and CE, HJ, KL and OP are of (12) type 
-@ and GI, HK and MO of (61) type at the more strongly curved general boundary. 


wlll 
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. (ii) 


nue Figure 4. (i) A plane ‘surface of misfit’, suggested by G. I. Taylor; (ii) and (iii) bubble-raft 

i | ‘erystal’ junctions, drawn from Dyson’s (1949) Figures 13 (a) and 11 (a) respectively 
(8 = 384°, [uv]=[21], [51], [91] or [13]; and 628°, [wv]=[71], [61] or [12] respectively, 
relative to orthohexagonal axes). In (iii) the lattices have a slight elastic shear causing 
curvature of the rows. A very small shear of this kind may occur in ‘soft’ crystals 
combined with partial relative rotational slip between successive layers. 


| §3. ELECTRON-DIFFRACTION EVIDENCE OF PREFERRED 
INTERCRYSTALLINE BOUNDARIES AND ROTATIONAL SLIP 


ve Typical examples recorded (except one of graphite) by the author and others 
(in this laboratory since 1934, are listed in Table 2 and illustrated in Figures 5 
m) to 55*. A Finch-type diffraction camera was used (Finch and Wilman 1937) at 
* =about 60kv. and 50cm. camera length. Details of specimen preparation and 
i) discussion of the effects of greater or less crystal thickness on the electron- 
()} diffraction patterns have been supplied but are not reproduced here owing to 
limits of space. It must be stated, however, that patterns obtained with the 
electron beam along the common lattice row which was the axis of the rotational 
displacement were almost exclusively used, and that of such patterns the 
following types occur: 

(i) Spot patterns from crystal pairs at 5 given by equations (1) and (2). 

] (ii) ‘Symmetrical’ triplets, quadruplets, etc., at equal intervals of 5 given 
“by (1) and (2). 

y (iii) ‘Unsymmetrical’ triplets, quadruplets, etc., usually at regularly 
decreasing intervals of 5 given by (1) and (2), and usually with the spots becoming 


sl * Most of these Figures are reproduced in the Plates, inserted at the end of the issue. 
i 
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progressively fainter apparently indicating thinner or laterally smaller crystal 
regions, or both. 

(iv) As (i) or (ii) but with additional spots due to strong secondary scattering, 
showing that the component crystal sheets are extensively superposed 
(Figures 26(a) and 31) (cf. Wilman 1940). 

(v) As (i) to (aii) but short or long arcs instead of spots, in (ili) often ending 
in a ‘ tailed’ arc. 


‘Table 1. Lattice constants of the materials investigated (in a. where in thick 
type, otherwise left in kx. where accuracy is low or doubtful). Sources 
of these data for most of the organic compounds are stated by Warren (1936) 
and ‘T'yson (1938) 
C=cubic; FC=face centred; BC=body-centred; ‘Tl =tetragonal; 
H = hexagonal ; R= rhombohedral; O = orthorhombic; M = monoclinic; 
Tri =triclinic i 


Material Symm. a b c B 
Cdl, H 4-24 — 6:84 = 
Graphite H 2°461 = == — 
NaCl FCC 5°640 -- —— a 
PbS FCC 5°929 a _— — 
MgO FCC 4-203 os — — 
NaNO, FCR 6-48 a—102° 40’ 
As,Og FCC 11-046 — — 
CuO M 4-563 5-910 5-108 99° 29’ 
Pb ferrocyanide (e - 8:28 — — — 
Cu, Ag, Au, Al, Pd, Pt FCC — = — = 
n—Cs.He¢ O 7-416 4-958 ~80 = 
N—Cyo Hyp O a TS ee > 
n—-C,,H;,;0H M 8-80 4-90 44-2 56° 40’ 
n—-C,,H;,0H M?; O? — — — — 
n-C,,H,,1 O? — — — — 
A long-chain wax O? _ “= — “= 
Methyl stearate O 7-48 5:02 — — 
Stearic acid M — — — — 
Elaidic acid ae 26°5 — 10:3 — 
Ferric stearate a — — = = 
Ba stearate H 4°85 4°85 — — 
Urotropine BEE 7-02: — = = 
Anthracene M 8°58 6:02 11-18 1253 
Naphthacene yen S 7-94 6:02 13:5 67° 50% 
(2.3 benzanthracene) ‘| = 80° 407 y=92° 40’ 
1 : 2 benzanthracene a 792 ore as 99 
M 8:21 6°53 48-4 90° 
ace he ace O 8-88 12:08 24-8 
5:6 cyclo-penteno-1:2 B.A. M 12:02 8:61 13-78 OSs 
1:2,5:6 dibenz- M 6:59 7-84 14-17 1035 
anthracene a 
Chrysene M 8°34 6:18 25-0 115-8° 
Oo pep ene 10 4:52 20-32 13-47 97-4° 
: O 7:59 7-69 22:38 — 
Me-cholanthrene a Mm é : ; a oe 1 ae 
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(vi) Single-crystal type of pattern but with tailed arcs extending from the 
| spots or replacing them, showing a range of orientation of the crystals about an 
axis parallel to the electron beam. 

(vii) Patterns of long arcs which decrease in intensity from a well-defined 
of head and are interrupted by clear-cut gaps which correspond to a main 6 interval 
“S3) (satisfying equations (1) or (2)) between crystal sheets meeting at a eu 

*) of the above type (Figures 45 and 49). 


Figure 26 6. C3,H¢_ doublet pattern and the spots due to secondary scattering (those not lying 
on the two large rectangular spot patterns). 


‘ © r e © 
: 8 e a) & 6 
mone e ‘ 
1 SoS me i} . . . 
e q=— == a | 8 e e 
Figure 45 (b). Sequence of spot pattern, gap Figure 49(b). Main centr: d-rectangle spot pattern 
6,, short arc, gap 82, and long arc tailing and typical members of the rotationally displaced 
off tadpole-wise, in 1:2,5:6 dibenz- arcs such as AA’, BB’, tadpole-shaped arcs such 
anthracene pattern of Figure 45 (a) as beyond A” and B’, and displaced pattern of 
(inner part only, shown). spots such as A’” and A’’”. 


As would be expected (cf. the familiar case of diffraction by twinned crystals), 
the patterns from the two crystals in contact do not differ much from the 
superposition, of the patterns the two would yield independently, with the 
addition of diffracted beams attributable to secondary scattering. Even in the 
superposed type of contact the resulting small atomic displacements in the 
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contact region are likely to be mainly normal to the contact plane (Figure 1 
and §2), i.e. parallel to the electron beam in the cases investigated here, thus 
causing no appreciable change in the diffraction pattern near the undeflected-beam 
spot. Extensive superposition of the component crystal sheets is proved in at 
least two cases (Figures 26(a) and 31) by the unusually strong spots also present 
due to secondary scattering, but when such spots are not detected the contact is 
still not necessarily (or may not be entirely) of the laterally adjacent type, and it 
may be partly of superposed type or even entirely so, the lack of secondary 
scattering being then merely due to the thinness of the two layers or their small 
proportion of overlap. Diffuse regions of increased intensity surrounding the 
arcs in intense patterns from beaten metal foils also appear to be largely due to 
such secondary scattering in superposed component crystal sheets. 

The nature of the patterns and the observed values of 6 prove that crystal 
groupings in contact, of one or both types considered in §2, often occur in the 
case of crystals grown from solution (CdI,, NaCl, NaNO, and the organic 
materials), cleaved crystal flakes (graphite and CuO), single-crystal films floated 
off soluble substrates (As,O,, Ag and Pd, prepared by sublimation on to an NaCl 
cleavage face), beaten metal foils (Cu, Ag, Au, Al, Pt, Pd), and rolled or beaten 
metal surfaces after recrystallization at bright red heat (Cu, Ag, Al). 

In all the examples observed the crystals had been subjected to deforming 
forces either deliberately (graphite, CuO, beaten or rolled metals) or 
inadvertently during growth (e.g. by surface tension and currents in solution, 
acting on tabular or dendritic crystals or growth layers on crystal faces) or after 
growth (for example, by forces on thin films due to contact with liquids, or by 
shearing of metal gauze supports, in As,O3, Ag and Pd). ‘These facts and the 
occurrence of the above seven types of pattern all point to the conclusion that the 
observed crystal groupings are most simply explained as more or less extensively 
superposed and cohering lamellae azimuthally displaced by a rotational slip of 
either symmetrical or unsymmetrical type, with some simultaneous or subsequent 
translational slip generally also present. 

Such a rotational slip, demonstrated macroscopically in §4, provides also a 
natural explanation of the origin of the nuclei which give rise at about 900°c. to 
secondary recrystallization in heavily cold-rolled sheets of metals like copper 
and aluminium, as observed in the present results and also those of Kronberg and 
Wilson (1949), who showed, indeed, that such secondary recrystallization in 
preferred orientations was associated with mechanical deformation at the sheared 
edges of their specimen sheets (cf. Evans, Layton and Wilman 1951). The 
cooperative small movements of atoms round the undisturbed coincident atoms 
which they suggested as origin (somewhat similar to those in the a—$ quartz 
transformation which have been not very aptly referred to as a ‘ Drehgleitung’ 
(Tertsch 1949)), seem much less probable than a rotational slip of lamellae as a 
whole. 

Further support for the present rotational-slip explanation of the crystal 
groupings illustrated is provided by the adhesion and continued growth of two 
pre-existing crystals in preferred azimuths when they come together on plane 
faces, as was observed by Gaubert (1896) and especially by Schaskolsky and 
Schubnikow (1933). The last two authors shook small octahedra or cubes of 
alum into saturated alum solution in which a large octahedral alum crystal was 
immersed with one face horizontal, and found that the small octahedra sometimes 
adhered and continued to grow parallel to the substrate crystal and sometimes 
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in the 60°-rotated {111}-twin position, while the small cubes tended to grow 
with a diagonal of the contact face parallel to an octahedral edge of the substrate 
crystal. Deicha (1948) reported similar groupings in NaCl crystallized from 
supersaturated solution. ; 


§4. DIRECT DEMONSTRATION OF ROTATIONAL SLIP 
AS A DEFORMATION PROCESS IN CRYSTALS 


w) (a) Macroscopic demonstration in potassium ferrocyanide trihydrate and gypsum. 
The results of §3 not only indicate the common occurrence of rotational slip 
in many types of crystals—especially by the observations of symmetrical and un- 
symmetrical triplets and multiplets, the clear evidence of superposition of the 
layers in some of these cases, and the often effectively continuous sequences 
(sometimes interrupted) in azimuth—but they also show that it occurs on densely 
populated planes which are often cleavage planes or else those planes which 
develop in contact with and parallel to the substrate during growth of crystals from 
solution, and which then also usually develop as main growth faces on the crystals. 
Cleavage and orientation faces are largely an expression of the cohesional 
anisotropy, in that the binding of the atoms or molecules is particularly strong 
laterally within each sheet parallel to these planes, and is weaker between atoms or 
molecules in neighbouring sheets. In rotational slip, as in cleavage, there is a 
step-by-step disengagement (cf. a ‘zip-fastener’) of the atoms of one sheet from 
the potential depressions of the neighbouring sheet, beginning with the atoms 
furthest from the rotation axis, but the neighbouring sheets are not pulled apart 
(or only slightly so) and retain mutual coherence such that many of the atoms again 
lie in or partly in potential troughs of the adjacent sheet. 
Potassium ferrocyanide trihydrate K,Fe(CN),.3H,O and gypsum, CaSQ,. 
2H,0 were therefore chosen as especially suitable for a macroscopic demonstration 
of rotational slip deformations, since highly perfect crystals 1-2 cm. in diameter and 
several millimetres thick were available and both have only one nearly perfect : 
cleavage face (cf. Buckley 1934), while the easy translational slip and flexure 
parallel to these cleavage planes is already known. The transparent yellow 
K,Fe(CN),.3H,O is monoclinic, a:b: c=0-3936: 1:0-2943, B=90° 2’, pseudo- 
tetragonal with base (010), and exceptionally perfect (010) cleavage (Groth 1906). 
Gypsum is also well known to be monoclinic with (010) as highly perfect cleavage, 
and its structure has been determined by x-ray diffraction (Wooster 1936, cf. 
de Jong and Bouman 1938). 
A shear couple was applied in the suggested optimum direction, 1.e. about an 
axis normal to the cleavage face. The couple was exerted by a frictional grip over 
the opposite faces held between two metal flats (paper or pitch covered) mounted in 
a lathe normal to the lathe axis. Figures 56 and 57 show typical results. 
The transparency of the crystals was practically unimpaired except where some 
fractures or surface abrasion occurred locally, and where the rotation had caused a 
large overhang of some parts of the crystal and led to slight curving over of these 
regions, which tended to develop small cleavage areas between layers having 
different flexure. The main area of the cleavage faces remained substantially flat 
i and the detail of initial small steps remained sharp. In the plan views of Figures 

i 56 and 57, taken witha black glass or paper background, the step edges on the upper 
and lower faces can be seen inclined at the angle of rotation of the faces, but both sets 
| are notin focus at once. Initially smooth side faces remained practically as glassy 
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as before, to the eye, but had developed a regularly twisted shape (Figure 56(a’)). 
This crystal shape cannot be accounted for except by a uniform relative rotational 
slip over each other of successive layers parallel to the cleavage plane, since the 
crystals were still cleavable easily on the planes parallel to the upper and lower 
faces, these cleavages being as smooth as initially though slightly undulating in 
crystals subjected to large shear (~40°/mm.). 

Crystals cut to roughly cylindrical shape to avoid overhang developing during 
shear also showed, after the rotation, a regular plastic rotational slip, so that shallow 
scratches initially parallel to the cylinder axis became inclined lines at practically 
constant inclination to the vertical (Figures 56 (c) and 56(d)). No appreciable 
change from this form of slip was noticed over a wide range of speed of shear. 

A striking difference of ease of rotation was found between the ferrocyanide and 
the gypsum. Inthe ferrocyanide, pieces of the order of 1 cm. square by 2-4 mm. 
thick could be rotated easily with the thumb, applying only about 10 kg/cm? 
pressure to obtain the necessary frictional torque. In gypsum a much higher 
pressure approximately 350 kg/cm? had to be applied before the frictional force 
was sufficient to cause appreciable rotation. In still harder crystals such as mica, 
calcite, zinc blende, quartz, sapphire, etc. rotational slip must be correspondingly 
more difficult to produce. 


(b) Interpretation of ‘deformation bands’ as initially rotational-slip lamellae. 


The lamellar ‘deformation bands’ or ‘etch bands’ observed microscopically 
in metals are less sharply defined at low deformation than normal translational- 
slip bands or twin bands and are best made visible by very slow etching. They have 
been investigated in detail by Elam (1935, 1936), Barrett (1939, 1940, 1943, 
1949), Barrett and Levenson (1939, 1940, 1941), Greninger (1938), Brick (1940), 
Brick and Williamson (1941), Crussard (1949) and Honeycombe (1950). 

Barrett, and Barrett and Levenson demonstrated the form of these bands 
especially clearly and showed that (1) they are parallel to densely populated net 
planes, e.g. (100) and (111) in Fe, (100) and (110) in Al and Cu, which are not the 
normal transitional slip planes—(110), (221) or (321) in Fe, (111) in Al and Cu 
—nor the twin planes—(221) in Fe, (111) in Cu; (ii) when the bands are first 
observed their lattice is rotated by several degrees about an axis normal to the plane 
of the lamellae (see Barrett and Levenson 1940*); (iii) at high deformation the 
bands show much flexure and their relative orientations become more complex; 
(iv) the sub-division of crystals into such lamellae, and the changes of orientation in 
these, are processes which tend to lead to development of preferred orientations 
in heavily cold-worked metal wires or sheets. 

Although Barrett and subsequent authors appear to have discussed these 
deformation bands only in terms of multiple translational slip on systems of the 
crystallographically equivalent planes (cf. Wood 1940) the author’s results above 
(§§ 2, 3, 4(@)) lead to the logical conclusion that these bands originate as rotational- 
slip bands before translational slip has caused much orientation change, though 
subsequent flexure and other complex forms of translational slip modify their 
form and orientation extensively at high deformation. Clearly, a crystal in a 
polycrystalline aggregate will be liable to be acted on by non-homogeneous shear 


* They state that at 28% reduction of an Al single crystal by compression they observed ‘‘ bands 
parallel to (010) planes, and the orientation difference across a band boundary was 3°”, increasing 
to 8° at 50% reduction. ‘“‘ For both compressions the two orientations have the [010] axis in common; 
that is, the pole of the plane of the bands.” 


| 
| 
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- couples, and if these cause rotational slip any subsequent translational slip, being 
_ stopped at the boundary rotational-slip planes of the deformation bands, must 
| cause lattice rotation relative to these boundaries and soon hinder further sub- 
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. division by rotational slip ; thus such bands would be expected to be formed mainly 
in the initial stages of deformation, as is observed. 


§5. THE RELATION OF ROTATIONAL SLIP TO THE FORMS OF 
TRANSLATIONAL SLIP AND TO DISLOCATION THEORY 


From §§3 and 4 rotational slip is seen to be analogous in form to translational 


_ slip in that (i) it occurs on densely populated net planes, not necessarily only on 


the usual translational-slip planes but particularly easily on cleavage planes 
(cf. also Evans and Wilman 1950, Evans, Layton and Wilman 1951); (ii) as in 
flexural translational slip, at the end of the slip a proportion of the atoms at the 
interface are again in potential troughs similar to those they occupied before; 
(iii) the slip lamellae, as indicated by the above electron-diffraction patterns 
are often only about 100a. thick (cf. also Heidenreich and Shockley 1947, 1948) ; 
and (iv) the primary slip planes tend to be regularly spaced and the slip 
displacements equal, i.e. ‘symmetrical slip’ (cf. Andrade 1935, Treuting and 
Brick 1942, Barrett 1944, 1947), though (v) unsymmetrical slip sometimes occurs 
under conditions which are at present uncertain. 

With respect to translational slip, as Barrett (1939) pointed out, the 
development of striations on some crystal faces is not in itself a criterion for 


_ pure translational slip, since it may be due to ‘deformation bands’. The 


unsymmetrical type of translational slip has not yet received much attention but 
seems to be well shown (with flexure finally, however) in a single-crystal rod 
of tin shown in Figure 58, illustrated by Masing and Polanyi (1923) and Polanyi 


a aaa 


Figure 58. Form of slip bands in tin in an example observed by Polanyi (1928). 


(1928, Figure 7). Schmid and Boas (1935) state that the extension of the crystal 
was very slow. The regular decrease of thickness of the slip blocks and the extent 
of slip are strikingly analogous to the typical decrease of spot intensity and 
azimuth intervals in the asymmetrical rotational slip in §3 (though in those of 
the latter examples which relate to growth from solution an explanation may also 
be possible in terms of dendritic growths subjected to forces which deflect the 
narrower growing tip progressively by rotational slip). 

It is of interest that in the frictional analogy of slip in a pile of cards an 
unsymmetrical distribution of translational or rotational slip displacements can 
be produced in a pile of similar superposed glazed cards (or even one of crépe 
paper sheets) subjected to only light normal pressure and slow displacement by 
tangential frictional force on the end faces; but in this case the thinner ‘ slip 
blocks’ near the moving face are displaced furthest, which is opposite to the 
results of §3 (unless the most displaced layers are also smaller in lateral extent) 
and the above example of tin. Larger pressures and speeds of slip led to a linear 
variation of slip displacement with height in the pile of cards, analogous to that 
in Figures 56 and 57 (cf. also lamellar viscous flow in liquids). 

Flexural translational slip is usually described as being characteristically 
about an axis F normal to the translational-slip direction T, and thus parallel 
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to the slip plane. It must obviously be possible, however, to bend a crystal 
about other axes parallel to the slip plane or inclined to it, and it seems likely that 
when the elastic limit is exceeded there may be some plastic rotational slip as 
well as translational slip. ‘This possibility has not yet been tested closely (but 
see Cahn 1949). 

Twist slip as described by Miigge (1898), Johnsen (1914) and Buerger (1930) 
with the translational-slip direction ‘Tas torsion axis, appears to be a special 
combination of translational slip and flexure of the slip planes. Closer 
investigation may well show that some rotational slip is produced either on the 
translational slip planes or on the planes normal to T’. 

Kinking (Miigge 1898, Friedel 1926, Brilliantow and Obreimow 1934, 1937, 
Orowan 1942, 1947, Crussard 1949, Barrett 1949) appears to be purely of a 
flexural nature and occurs when the axis of compression is roughly along the 
translational slip planes. ‘The present writer finds indeed that a close analogy to 
Orowan’s (1942, Figure 1 (a)) kink pairs in cadmium can be easily produced by 
lengthwise compression of a pile of paper strips formed by the remaining top end 
of a tear-off writing pad. Crussard (1949) has suggested that the change of 
direction at a kink may occur sharply on an atomic scale so that preferred kink 
angles occur and that the kink planes in cadmium shown by Orowan’s figure 
may be {1011}. However, though Orowan only stated (0001) was parallel to the 
wire axis, it appears that the axis must have been near (100), since that is the 
direction of fastest growth in (0001) and is the translational-slip direction. Thus 
the kink planes could more possibly be {1122}, giving a kink angle 6=117° 56’ 
(Figure 59), which agrees reasonably well with the estimated angle 116° 
(Crussard’s estimate was ¢ = 113° on one side, 118—120° on the other). 


Figure 59. Possible diagrammatic section of a kink, the kink angle ¢ being then given by 
tan ¢/2=n,d/ngt, where n, and ng are integers. 

Though kinking thus seems capable of giving rise in some cases to crystals 
in the ‘laterally adjacent’ type of contact, it seems improbable that thin but 
laterally extensive crystal sheets growing from solution would kink about an axis 
normal to their plane, except possibly in their initial stages when they were still 
small. Even then in many of the cases observed in §3 the rotation axes are not 
those which would correspond to kinking related as above to the normal 
translational-slip direction. Inhomogeneous shear and break up of thin crystal 
sheets followed by re-coalescence during growth (cf. the soap-bubble sheets, 
Dyson 1949) might also give ‘laterally adjacent’ contacts at the preferred angles 
defined in §2 and thus lead to electron diffraction spot patterns like some in § 3, 
but seem unlikely to give rise to the practically continuous azimuthal sequences 
which yield the arc patterns observed in §3. 
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Like the experimental knowledge of slip, the theoretical treatment in terms 
of ‘ dislocations ’ also appears as yet to be rather imperfect. According to Frank 
(1948, 1949) the ‘proper description for a rotation of two lattices on a like contact 
plane is an accumulation of screw dislocations (Cottrell 1949) in a crossed-grid 
array in that plane. 

Van der Merwe (1950) has illustrated the necessary slight atomic 
displacements (laterally) from their lattice-point positions in two superposed 
plane atomic lattices of square type, according to such a set of dislocations. His 
calculation, on this basis, of the potential energy of two such superposed atomic 
planes in relation to the angle of displacement 6 of one atomic sheet over the 
other, shows no minima such as were indicated by the author’s ball-bearing 
model (§ 2) and electron-diffraction data (§3), and this shows that the calculations 
are incomplete in detail. If van der Merwe’s calculated smooth variation of 
potential energy with 5 were correct it is evident that two crystals in such 
superposed contact must automatically rotate into parallelism, but this is shown 
not to be so by the present results. It is thus suggested that two densely populated 
plane atomic lattices rotated over each other in contact undergo chiefly only 
adjustments of atomic positions normal to their initial plane array. 


§6. THE RELATION OF ROTATIONAL SLIP TO THE GROWTH 
AN DEP ROPER ELE SMO CR Y.S tAViS 


(a) General nature of growth. 


The results of §3 show that rotational slip often occurs in crystals grown from 
solution, and is generally along the atomic or molecular layers which are first formed 
on the substrate and later covered by similar successive layers as main boundary 
faces of the growing crystals. The optically visible thick layers (cf. Bunn and. 
Emmett 1949) growing across the faces of many crystals may well owe their origin 
to occasional small rotational slips in their early stages of development. The 
experiments of Bernauer (1926, 1927, 1929) and Rossmann (1934) on growth of 
‘twisted’ and ‘spiral’ crystals suggest that vigorous stirring or agitation causes 
such a disturbance in the superposition of layers on the growing crystals (cf. §3 
above), which then possess abnormal optical polarization properties (as was found 
in paraffin crystals by Hubbard (1945), who, however, overlooked its true origin). 
Their illustrations appear to represent dendritic crystal forms in many cases, and it 
will be evident that repeated rotational slip on planes parallel to the substrate, 
especially across the narrow tips of the growing dendrites, could lead to the 
observed crystal shapes, which may correspond to the electron-diffraction patterns 
of ‘tailed’ arcs or spot-groups (§3). 

In close relation to such a process, the observation by Balarew and Kolarew 
(1939) that a variation of rate of dissolution at crystal faces (gypsum) occurs 
according to the direction and speed of stirring seems explicable in terms of slip 
processes caused thereby (rotational or translational) in surface regions. 

Some striking macroscopic twisted crystals (Spencer 1921, Frondel 1936), 
for example of quartz and pyrites, probably originate by a rotational-slip deforma- 
tion after growth, this being made possible by the crystal being embedded ina 
strong retaining matrix hindering fragmentation during the torsional deformation. 

Vicinal faces are explicable (Bunn and Emmett 1949) as a more or less regular 
terracing of superposed growth layers on main crystal faces (cf. Wyckoff 1948), 
but some reported cases could correspond to mosaic components related by 
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rotational-slip contacts, e.g. (1.0.10) of galena (Hintze 1898), which is at an angle 
of tan-! (1/10) 2 tan“ (1/20) to the nearest cube face, could be really a (100) face 
of amosaic component which was first formed from the main crystal bya rotational 
slip on (010) to the position defined by u,/v,=20 in equation (1). 

In crystals having well-defined faces other than the rotational-slip plane, a 
rotational slip would cause light reflections to be elongated into streaks or, if slip 
were only at relatively large intervals, rows of separate reflection spots (Brilliantow 
and Obreimow 1934). In x-ray diffraction, similarly, the asterisms and associated 
trains of diffraction spots observed from some mosaic crystals (e.g. NaNO, grown 
from the melt, Leonhardt and 'Tiemeyer 1936) and in some cases of plastically 
deformed crystals are almost certainly due in many cases to rotational slip in the 
crystals. 

Many other observations and properties of crystals will no doubt be clarified 
when considered in the light of possible presence of rotational slip. 


(6) Rotational slp and epitaxial crystal growth. 


The results of §§ 2-4 provide a fundamental new insight into the conditions 
which determine epitaxial crystal growth and its stability and degree of perfection. 

(i) The effect of cracks or steps tn the substrate. ‘The rotational-slip results show 
in the clearest possible way that the existence of angles or edges of cleavage steps, 
whether these are macroscopic or only monatomic in height, is quite unnecessary 
for epitaxial orientation to occur as suggested by ‘Thomson and Cochrane (1939) 
and Thomson (1948). When such steps or cracks are present they may, however, 
often determine the location or habit (cf. Drabble 1949) of the deposit crystals. 

(ii) The origin and nature of epitaxial growth. ‘The epitaxial growth of one 
crystal on another is thus determined by the tendency of the atoms (or ions or 
molecules), on arrival on the substrate, to be held in the potential troughs, which 
exist in a periodic array on a flat substrate crystal face, though the deposit atoms 
have on arrival an initial kinetic energy and they retain on the average a certain 
mobility over the substrate, characteristic of the nature of the deposit and substrate 
and the temperature of the surface. It is then to be expected that, if the tempera- 
ture is not so high that recrystallization simultaneously occurs, the deposit crystals 
must be built up at first in a more or less strained or ‘ pseudomorphic’ form as Finch 
and Quarrell (1933) found and named it. The orientation is determined by the 
arrangement of the first layer or so of deposit atoms or molecules in positions of 
lowest attainable potential energy on the substrate (Elleman and Wilman 1948) 
and this in turn is determined mainly but not quite entirely by that of the substrate 
surface layer of atoms if this is a densely populated atomic plane and other 
conditions are appropriate. 

(iii) Lattice relationships necessary for epitaxy. Firstly, though crystal growth 
is in general a process of addition of atoms (or molecules) separately on to the first 
formed small groups of atoms, 1.e. the crystal nuclei, whereas rotational slip involves 
movements of sheets of atoms (or molecules) as a whole, it is quite clear from §3 
that from the /attice-fitting point of view there is no hindrance to epitaxial orien- 
tations being developed which have not even one densely populated lattice row 
parallel to one in the substrate. However, such cases must clearly be much rarer 
than those with at least one densely populated row in the deposit crystal parallel to 
one in the substrate, as is found in practice. As Elleman and Wilman (1949) 
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| emphasized, the epitaxial orientation is in general determined by two distinct 
| factors—initially by the fitting of the deposit atom groups with the local atomic 
| distribution in the substrate surface so that low mutual potential energy is attained, 

_ then by the tendency to satisfy as far as possible the mutual fitting together of the 
) longer range lattice periodicities. ‘The above rarer case will thus mainly be expected 
_ only when the substrate and deposit crystal lattice periodicities are large and small 


respectively, so that the deposit crystal orientation is initiated as usual by the local 
substrate-surface atomic distribution and the crystals can grow big enough to be’ 
relatively immobile before the stresses caused by the lattice misfits tend to modify 
the initial nucleus orientation. 

(iv) Expitaxy of pre-existing crystals. ‘The rotational-slip results are largely 


_ analogous to, and further supported by, the observations of preferred azimuths 


of adhesion and further growth (from solution) of pre-existing crystals which 
had plane faces in contact with a like or unlike plane face of a larger crystal (see § 3). 
The necessity of considering such an origin in some cases of epitaxial growth, 
such as (001) KCl growth on (001) of mica, reported by Deicha (1946), has been 
pointed out by Royer (1948), though the present rotational-slip results show that 
it is in general unlikely that only one azimuthal type of orientation would result 
by this process. 

(v) Stability of expitaxially orientated crystals: lattice dimension relationships. 
The degree of strain and stability of orientation of an epitaxial growth depends 
partly on the nature and degree of the lattice fitting and partly on the local atomic 
fitting (see (iii) above). As in the case of rotational slip, the larger the density of 
coincidence points (or close approximations to such) of the substrate and deposit 
lattices, per unit area of the interface, the stronger the adhesion of the crystal 
to the substrate crystal and also the more stable its orientation since there will be 
less strain (or more evenly distributed strain) due to the atomic misfits, and thus 
less tendency to disorientation by plastic deformation or by recrystallization. 

Clearly from this point of view, high stability and perfection of the orientation 
(without recrystallization occurring) are not associated with the existence of a 
lattice row in the deposit crystal with a periodic spacing which is roughly the 
same as (say within 15°% of) that of the parallel lattice row in the substrate. It 
must, on the other hand, be associated with the occurrence of rows such that 
x unit translations of one row closely approximate to y units of the parallel row, 
x and y being small integers. Numerous cases of this last type have been observed 
in this laboratory, and others have been noted occasionally elsewhere. 

(vi) Release of expitaxial strain by slp deformations. When the deposit-crystal 
thickness and extent reach a certain value, dependent on the lattice constant 
difference, the strain must be relieved by slip at the junction region, as Frank 
and van der Merwe (1949) and van der Merwe (1950) pointed out and discussed. 
The normal crystal structure thus tends to be attained though, as is now seen more 
clearly, in some cases at the expense of perfection of orientation or the development 
of mechanical twinning, the crystals, then usually still quite small, continuing 
growth in the modified orientations. : 

The {111}-twinning in face-centred-cubic metals such as silver, formed from 
the vapour on alkali-halide {001} substrates of NaCl type and on metal substrates 
by electro-deposition, may thus arise in a sudden release of strain in the thin 
initially pseudomorphic crystals if these were in parallel {001} orientation. In 
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PbS deposits on {001} NaCl (Elleman and Wilman 1948) the {111} and {332} 
twinnings observed were also types which have been observed to be producible 
by mechanical deformation and may have arisen in this way. 

Since this paper was first submitted for publication the first example of 
rotational slip (again on cleavage planes, cf. §4) due to expitaxial stresses has 
been found by Evans and Wilman (1950) in ZnO formed on zinc-blende {110} 
cleavage faces at 520°c. in air. At red heat no such disorientation was observed, 
probably because at this temperature the ZnO crystals formed underwent slip 
more easily, thus giving earlier releasing of the strain (and/or a recrystallization 
readjustment of atomic positions to more normal structure in the pseudomorphic 
nucleus) before the crystal size became large and much stress accumulated. 
A similar lack of disorientation was observed in Cu,O formed anodically at room 
temperature in parallel growth on a Cu (110) face (Goswami 1950). 

Other examples of directed disorientation about well-defined axes have now 
been found in this laboratory with the collaboration of Dr. A. Goswami (1950), 
Dr. M. K. Gharpurey (1950), Dr. D. M. Evans and Dr. K. Leu. In many other 
cases the electron-diffraction patterns show arcing at all azimuths because of the 
high symmetry of the substrate or deposit crystals, which leads to many different 
but equivalent azimuthal orientations in the deposit, or several equivalent but 
difterently orientated possible rotational slip planes. 


(c) Deformation of crystals by unidirectional abrasion. 

Results indicating prominent rotational slip caused by unidirectional abrasion 
of crystals have been noted (Wilman 195(C, Evans and Wilman 1950). A fuller 
description of the new results from copper and iron is given by Evans, Layton 


and Wilman (1951). 
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ABSTRACT. Strong strip specimens are prepared from CdO powder such that they 
possess an electrical conductivity of the order of 300 ohm~tcm™! at room temperature. 
Simultaneous measurements of the Hall constant R and the electrical conductivity o are 
made throughout the temperature range 0-500° c., and the thermoelectric power dE/dT 
is measured immediately afterwards on the same specimen. The properties are found to 
resemble those of a metal: Ris found to have a constant value independent of the tempera- 
ture. dE/dT is approximately 60 uv/deg. c. at room temperatures, whilst o decreases with 
temperature, its value at 500° c. being roughly half that at 0° c. 

The theory of degenerate semiconductors is developed, and since the values of R indicate 
that there is a constant number of free electrons in each specimen, this result is incorporated 
into the comparison of the theory with the experimental results. Satisfactory agreement 
is obtained between theory and experiment. Agreement is also obtained when the theory 
is compared with results obtained by another investigator for the thermoelectric power of 
CdO. The approximate form of the theory predicts that the thermoelectric power and the 
resistivity should be directly proportional to the absolute temperature. The experimental 
results are found to confirm this. CdO behaves as an excess degenerate semiconductor, 
with properties strongly resembling those of a metal. 


§1. INTRODUCTION 
EMICONDUCTORS may be defined as ‘“‘ materials whose conductivity is 
due to impurities or non-stoichiometric excesses of one or other of the 
components of the crystal lattice’’. Purification (removal of defects from 
the lattice) causes the conductivity to decrease. This definition covers a very wide 
range of substances which possess very varied properties. It therefore seems 


* The major part of the subject matter of this paper formed part of a thesis which was accepted 
by London University for the degree of Ph.D. 
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necessary, in this paper, to subdivide electronic semiconductors, as defined 
above, into two groups: (a) semiconducting materials which have so few 
electrons that the classical distribution of energy amongst the electrons may be 
used with an accuracy closer than 5°% will be termed ‘ classical semiconductors ’ ; 
(6) semiconducting materials which have so many free electrons that the 
Fermi—Dirac statistics must be applied to the electrons, will be termed ‘ degenerate 
semiconductors’. : 

The simple theory of semiconductors as developed by Wilson (1931 a, b, 1939) 
and Fowler (1933 a,b, 1936) is applicable only to classical semiconductors. 
It has been shown that substances such as germanium, silicon, tellurium, silicon 
carbide, lead sulphide, and cadmium oxide, under certain conditions must be 
included in amongst the degenerate semiconductors. Busch and his collaborators 
(1946-1947) have shown that where the free electron density in their specimens 
of silicon carbide causes appreciable degeneracy, the employment of the 
Fermi—Dirac statistics brings theory into line with experiment. Similarly, 
Putley’s measurements (1948, 1949) on very degenerate specimens of germanium 
and silicon were found to be consistent with the theory, provided that the 
approximations of the classical theory were not introduced. 

Because of its interesting electrical properties, which are totally different 
from those of the majority of investigated semiconductors, cadmium oxide has 
already been the subject of several investigations. In the majority of these only 
one electrical property, either conductivity or thermoelectric power, has been 
measured (Fischer et al. 1932, Andrews 1947, Andrews and Hogarth 1949). 
Hogarth (1948) studied both these properties, but not simultaneously, nor upon 
the same specimen. General conclusions drawn from such incomplete data 
must always remain doubtful, since the properties may vary greatly from 
specimen to specimen, depending upon their origin and heat treatment, and 
discussions of measurements on two specimens, whether of the same property 
or of two separate properties, must remain inconclusive. 

The need for a more comprehensive investigation of cadmium oxide led to 
the present work, in which the Hall constant, the specific conductivity, and the 
thermoelectric power, with their variations with temperature, were measured, 
all upon the same specimen, and as far as possible simultaneously. In this way 
it was hoped that correlations between the properties might be attempted, and 
also that the theory might be developed to account for the variations. Measure- 
ments of the Hall constant especially were aimed at, since they give direct 
information about the conduction, and no previous account of any measurements 


could be found. 


§2. PREPARATION OF THE SPECIMEN 

Cadmium oxide crystals were not available for the investigation ; furthermore, 
cadmium forms a protective layer of oxide when oxidized, preventing further 
oxidation. It is possible to evaporate a thin film of cadmium upon a surface, 
and then oxidize it, but such a process would give all the added complications 
of working with thin films. It was therefore decided that a compressed and 
sintered specimen of CdO powder should be the subject of this investigation. ‘The 
specimen was made rectangular in shape, approximately 7 cm. x 1:25 cm. x 0-1.cm., 
by compression in an appropriately shaped steel mould. ‘The compression was 
performed by a hydraulic press, which would produce a force of about 7 tons. 
The powder, the purest obtainable, bound together well though the specimens 
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produced were necessarily very brittle due to their thinness. The specimens 
were then sintered in the atmosphere, inside an electrical furnace, at about 
720° c., for one or two days. The resultant strips of the oxide were quite strong, 
and gave a definite metallic ring when dropped. Furthermore, as will be shown 
later, the specimens, formed as above, consisted of some stable structure Cd,O, 
where x is slightly greater than unity. The densities of the specimens were of 
the order of 6gm/cm’, whilst the crystal density is 8 gm/cm?. 


§3. METHOD OF MEASURING HALL CONSTANT AND 
CON DWC my Ay: 


The apparatus used was essentially the same as that described in a previous 
paper (Wright and Andrews 1949) and found very successful for use in similar 
work upon nickel oxide. The flat specimen was held against a plate of 
‘ Syndanio ’ (a type of asbestos slate) by means of the platinum electrical contacts, 
which consisted of two current leads at the ends of the specimen, two potential 
probes about 5cm. apart on either side of the centre of the specimen along the 
length of the specimen, and two Hall £.M.F. point contacts on either side laterally 
of the centre of the specimen, about 1:25cm. apart. The ratio of the distances 
apart of the potential probes and the Hall effect contacts, 4 : 1, is sufficient to 
eliminate the effect of size upon the Hall E.mM.F. measured. The temperature 
was measured by means of two thermocouples, one in contact with each potential 
probe. ‘The system, fixed upon a holder, between nichrome ribbon heaters, 
was less than lcm. thick, so that the poles of the water-cooled electromagnet 
between which it was inserted could be close together, and a magnetic field of 
12-5 kilo-oersteds, virtually constant over a region 3cm. wide, could be 
maintained. 

The current through the specimen was produced by a source of E.M.F. in 
series with the specimen, a variable resistance, and a standard 0-01 ohm resistance. 
The current could be varied from 0-1 to 2-0amp., but normally a current of 
0-8 amp. was used, which caused very little heating in the specimen. 

All air currents were carefully shielded from the specimen by enclosing 
the whole apparatus. ‘To perform the measurements, the system is first allowed 
to attain thermal equilibrium at the required temperature, then the potential 
across the standard resistance is measured with the potentiometer, giving a value 
for the current flowing in the circuit; the potential across the potential probes 
is also measured and the specific conductivity calculated. ‘The magnetic field is 
switched on, and the E.M.F. across the Hall effect probes is measured using the 
potentiometer. ‘This is composed in the main of two components, one due to 
the Hall effect, and the other due to the fact that it is impossible to set the probes 
exactly upon an equipotential surface; the E.M.F. is therefore again measured 
with the magnetic field reversed, when the former contribution is reversed, but 
not the latter. The Hall constant may be calculated from the readings and the 
dimensions of the specimen. 

A slow continuous change of potential of a few microvolts between the Hall 
contacts was found to occur after the field was reversed. The change was so very 
small and very slow, however, that by making instantaneous readings of the 
potential difference when the field is reversed this source of potential, which is 
due to the Ettinghausen effect, was virtually eliminated. 
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§4. METHOD OF MEASURING THE THERMOELECTRIC POWER 
| The thermoelectric power was measured by inserting the specimen 
‘efi unsymmetrically in a horizontal tubular furnace, so that a temperature 


ili difference of 20-40° c. existed between the ends of the specimen, as recorded by 
‘oj thermocouples held in firm contact with the specimen beneath the platinum 


potential probes, across which the thermoelectric E.M.F. was measured by the 
} potentiometer. In obtaining reproducible results it was found to be extremely 


@ important that all air currents were excluded from the specimen by enclosing 


the furnace, and before any measurements were made, time had to be allowed 
) for complete thermal equilibrium to be attained. 


§5. RESULTS AND COMMENT ON THE METHOD 


The essential differences between the apparatus described above aad that 
described in a previous paper (Wright and Andrews 1949) are due to the very 
great differences in the specimens which were used. The specimens of CdO 
had resistances of about lohm at room temperature, whilst at similar 
temperatures the specimens of NiO had resistances of the order of 10!°chms. 
Thus the electrostatic instruments which had to be used for measuring potentials 
with NiO could be discarded in favour of the more accurate and easily usable 
+ potentiometer for CdO. 

The object of this investigation was to take a specimen of CdO with a constant 
composition, i.e. constant impurity content of interstitial Cd atoms, and to 
examine its electrical properties. It was known, however, that CdO was unstable 
at about 600°c. and above ; the temperatures of the measurements were 
therefore restricted to below about 550°c. Above this temperature it was 
considered that a gradual change in the composition of the CdO would occur, 
a certain amount of the oxygen content being driven out of the lattice, leaving a 
greater excess of Cd atoms, or impurity centres. Such a change appears to be 
reversible, and upon cooling oxygen is replaced in the lattice from the 
atmosphere; this is indicated by the fact that above 600°c. the CdO changes 
colour to black, and upon cooling the original light brown colour returns. Thus 
no permanent change will be produced, and the variation of the constitution 
will be very difficult to allow for in the interpretation of the results. By limiting 
the temperature to 550°c. this difficulty is not allowed to arise. 

Preliminary experiments upon the mode of conduction, by Mr. E. W. Bastin 
(private communication) verified the fact that the conduction in CdO obeys 
Ohm’s law, so long as the specimen is dry; this was taken to indicate electronic 
conduction in the dry state. To ensure that the specimens were completely 
dry when used, before a series of experiments was begun, the specimen was 
heated to a temperature in excess of 100°c. for about 24 hours. The actual 
series of readings was then performed at temperatures from 0-550°c. and the 
readings were repeated as the specimen cooled down to room temperatures. 
No permanent changes were found to have resulted from an experimental run. 

Two specimens were investigated thoroughly, and the results were so similar 
and easily reproducible, that it was considered unnecessary to repeat upon more 
specimens, particularly since initial measurements made upon two other 
specimens, which were later broken before the series could be completed, 
yielded similar results. 
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The direction of the Hall £.m.F. and the thermoelectric E.M.F. were such as 
to indicate conduction by negatively charged carriers, i.e. electrons. The results 
of the measurements of the Hall constant, specific conductivity, and thermo- 
electric power at various temperatures are shown graphically in Figures 1, 2, and 3 
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§6. DEVELOPMENT OF THE THEORY OF ELECTRONIC 
CONDUCTION IN AN IONIC LATTICE 


(1) Concentration of Free Electrons 


Consider the free electrons in the conduction band of a substance with an 
79 ionic lattice. The number of states with energies between EF and E+dE above 
} the bottom of the conduction band (Seitz 1941) is 


4 2m* 3/2 
N(E)dE = ee EVdE. 


: where m* is the effective mass of the electrons. Thus the number of free electrons 
i) with energy between E and E+dE is given by N(E)f(E)dE, where f(£) is the 
#} Fermi—Dirac function representing the probability of an electron occupying a 
4 state of energy E. The total number of free electrons present is therefore 

| 4n(2m*)3?2 i E\2dE 


h® 9 CBRE 4 1? 


* where £* is the normalizing parameter of the Fermi—Dirac function, and is the 
“%) chemical potential of the free electrons. Let the reduced energies be 7 = E/RT, 
m™ and 7»* =E*/kT, then 


a | 3 N(E)f(E)dE = 


- 4n(2m*RT)32 
he DEST Aa a cater ee (1) 
ea) nd 
where | Bagh | ei nos iui eanar (2) 


(ii) Mean Free Path of the Electrons 


If 1/7(£) is the probability that an electron with energy EF above the bottom 
of the band will suffer deflection in the next unit time, then 7(£) is termed the 
relaxation time, and it has been shown (Frohlich 1937, Mott and Frohlich 1939, 
Davydov and Shmushkevitch 1940) that for electrons with energies E greater 
than hv, where v is the frequency of the characteristic vibrations of the lattice, 


CEM 
7(£) = gi/kT 7? jee lare 6.6 (3) 
2 \12 hDva* 1 1 1 
where C= (=) ok and D = De ai De 


where D., and D_ are the masses of the positive and negative ions of the lattice, 
respectively, and a is the lattice constant. 

The mean free path of an electron with energy F is therefore given by 

KE) = (Zim) Eter(E)=lyBy, sos! = anni) abi! (4) 

2RT/hv 
1+2/{(emr? —1}" 
The term B varies very slowly with the temperature, provided that hv/RT is 
less than unity. For values of T much greater than hv/k, B=1 and is constant. 
When T=hy/k then B=1-08. Thus, for temperatures which are greater than 
hv/k, we may consider B as approximately constant and equal to unity. 

If we further restrict ourselves to cases in which the distribution of free 
electrons is such that the vast majority have energies greater than hv, as may be 


where /) is a constant, and B= 
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caused by a high concentration of free electrons, or high temperatures, or a jj | 
combination of both, then we may consider that equation (4) may be used as a 
reasonable approximation for all the electrons. 


(iui) Electrical Conductivity 

From the Sommerfeld—Lorentz theory of electronic conduction in solids } 
(Seitz 1941) the following expression may be obtained for the specific . 
conductivity : 


o=— 


*oe2 70 
l6nm*e f aE) ap. 


SRE dE 
substituting for / from equation (4) and expressing in reduced energies 
% e2 co 
6 tone fae | pED) gio. BTE Ge eee (5) 
3h 0 dn 


(iv) Thermoelectric Power 


The Sommerfeld—Lorentz theory also leads to the following expression for 
the thermoelectric power © using reduced energies: 


" edfidn)dy | 
Cg tee he. | season Bos ia eee (6) 
| Inlaflaayer 


Substituting from equation (4) for the mean free path we may show that 


@= Efe -1}- Bea 


(v) The Hall Constant 


The Sommerfeld—Lorentz theory leads to the following expressions for the 
Hall constant R and the conductivity o(H) in a magnetic field H, 


Obit \ es ae "et fF 2 \N2 at 

pee) Ces eee ae 2S (| eae 

e (7) HL2+L2 Bus Ge (=) Le 
_ 41r(2m*)8 ¢° El(df/dE)dE _ 4n(2m*)8? (° Els( df/dE)dE 
where T=" |, — tal ey eee 
‘and s=—elH/m*v where v is the velocity of the electrons. Therefore 


Ro(H)=—L,/L,H. For weak magnetic fields s? is much less than 1, and o(H) 
approaches o. ‘Thus neglecting terms in s®, using reduced energies, substituting 
from equation (4) for J, from equation (5) for o, and introducing the number 
of free electrons n, from .equation (1) 


— 115 Fin(n*)Fae(n®)_ 4 |. (8) 
nel6  [F(n*)]}? Mee” 


(vi) Conclusions 


Equation (1) gives the number of free electrons in any electronic semi- 
conductor. Equations (5), (7), and (8) for the conductivity, thermoelectric 
power and Hall constant are applicable to conduction by electrons in ionic 
lattices, with the restrictions that must be applied to the use of equation (4) for 
the mean free path, viz. that the concentration of free electrons and the 
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# exceed Ay, and that the temperature is equal to or greater than hv/k, so that the 
‘ factor B is approximately constant. No derivation of the equations has been 
| previously given, but Shifrin (1944) has quoted similar results. 


$7. DISCUSSION OF THE HALL CONSTANT RESUETS 
AND THE THEORY 
Values of Fy/.(n*) and F'3).(n*) have been tabulated by McDougall and Stoner 
(1938) for various values of 7*, and the author has tabulated F(7*) and F,(n*) 
{ over the same range and at the same intervals. Values of the latter integrals for 
| integral values of 1* are shown in Table 1. 


Table 1. Calculated values of F,(y*) and F,(y*) for integral values of 1*, with 
the calculated values of the thermoelectric powers as predicted for ‘ high 


temperatures’ from equation (8), and for ‘low temperatures’ from 
equation (11). 


* * * Ohigh Olow 

u) Fy(n*) F(7*) (uv /°K.) (uv/°K.) 
—4 0:0182 0-:0366 = 558:0 
—3 0:0492 0:0990 515-0 474:2 
—2 0-1310— 02662 431-9 393-2 
—1 0-3387 0:7050 352°8 SI 
0 0-8225 1-8030. 2788 ASS 
1 1-8062 4-3282 222:0 209°9 
A 35135 silts) 17227 178-0 
3 60957 18-9686 142:8 157°8 
4 9-6267 34-529 118-5 144-0 
5 14-1382 58-130 99-77 135-0 
6 19-6424 91-744 86-24 128-5 
7 26°144 137-364 75°39 123" 
8 33-645 196 -987 67°25 120-0 
2 42-145 272-609 60-04 —= 
10 51-645 368 232 54:59 116°5 
11 62:145 479-856 49-72 — 
12 73-645 615-480 46-20 = 
13 86-145 774-1 42-70 111:4 
14 91-645 960:7 39-98 = 
15 114-145 1174-4 S122 => 
16 129-645 1418-0 34°65 107-0 
17 146-145 1693-7 32-84 “= 
18 163-645 20033 31-02 = 
iS) 182-145 2348°8 29:33 a 
20 201-645 ASG 28-13 = 


From these integrals the values of the constant A in the equation (8) may be 
calculated for various values of * (Table 2). 


Table 2 


n*  —00 8 0 5 10 20 +00 
A 11045 1-1021  1-0858  1-0221 1:0074 1:0020 41-0000 


It will be seen that the values of A vary only by 10°% for the whole range of * 
and the variation is extremely gradual. Furthermore, it will be shown that over 
the range of temperatures 0-550° c., the value of y* appears to vary between 5 and 
10 for CdO. Treatment of A as a constant and equal to unity will therefore be 
likely to incur an error of not more than about 2%. 

The values of the Hall constant measured, and the values of the concentration 
of free electrons m,, as calculated from equation (8) assuming A =1, are shown 
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graphically in Figure 1. It will be seen that the concentration of free electrons 
in the samples is constant. ‘This experimental result will be considered below, 
together with its probable explanation. 

It has been suggested previously (Peierls 1937, Wilson 1939, Shifrin 1944) 
that the activation energy in a semiconductor should vary with the impurity 
content. As the concentration of impurity centres increases, they will become 
closer together, and interaction will occur so that the potential lattice no longer 
consists of wells below the conduction band at the centre, but a series of potential 
~ humps are formed between the centres. The region between the tops of the 
humps and the bottom of the conduction band consists of a band of discrete 
energy levels due to the excited states of the impurities. Thus the height of the 
conduction band above the ground state levels is decreased. Eventually, with a 
sufficiently large impurity content, the conduction band will be lowered down to 
the ground state, when the impurity electrons will be continuously free. Such a 
material will have a constant concentration of free electrons, independent of 
temperature, as appears to be the case in cadmium oxide. This simple picture 
has, however, several difficulties which more rigorous treatments have as yet 
failed to overcome. 

Such a variation of activation energy has been observed by Pearson and 
Bardeen (1949) and quantitatively explained, using specimens of silicon with 
varying amounts of impurities. The activation energy decreases with increasing 
impurity content, until eventually it becomes zero, the electrons being 
continuously free in the ‘impurity band’. Thus it seems reasonable to suppose 
that a similar state of affairs exists in the CdO specimens investigated, there being 
such a concentration of excess interstitial cadmium atoms, that the ground state 
‘orbits’ of the centres overlap, and the electrons may pass freely through the 
crystal lattice. 


§8. COMPARISON OF THEORY WITH EXPERIMENTAL RESULTS 


In this comparison the constancy of the electron density will be taken as 
proved. The comparison will not depend upon the correctness or otherwise of 
the explanation of this fact given in the previous section, but merely upon the 
interpretation of the Hall constant results. 

Let the effective mass m* of the electron be equal to Gm, where m is the mass 
of a free electron, and Gis a constant. Then we may rewrite equation (1) 

hé Neg 
Fiy(q*) = 5s ser S. 
TAG ) 4n(2MRT)32 G32’ (9) 
so that from a knowledge of either 7* or n,/G?? (effective electron density), the 
other may be calculated at any temperature. 


(i) Thermoelectric Power Results 


The thermoelectric powers for the various values of 7* have been calculated 
from equation (7), and the results are shown in Table 1, and Figure 6. From a 
selected experimental point giving the thermoelectric power at a certain 
temperature, we may deduce 7* and hence n,/G?._ From this value of n,/G32 
the thermoelectric power may be calculated throughout the temperature range, 
and the results of such calculations are shown as the continuous lines in Figure 3. 
It will be seen that the theory predicts the variation of the thermoelectric power 
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j2xperimental and the theoretical curve. Once again the agreement is good. 

| It is of interest to note that the relationship between thermoelectric power 
jand temperature in the range considered is practically linear, both experimentally 
Sand theoretically. 


(11) Conductivity Results 


|) Assuming a value for n,/G?’, we may calculate 7*, and thus F,(7*) at two 
different temperatures T, and T,. Hence the ratio of the conductivities at the 
aitwo temperatures o7/o7, may be calculated from equation (6): 


op, /or, = Ty Fy(1*)/ TF (79*). 


i'The values of this ratio for a range of values of n,/G?? and for T, =300°xk. and 
¢ 7T,=800°. are given in Table 3. 


Table 3 
# Log (n¢/G?2) 20:5 20-6. 20-7 20:8 20-9 21-0 24-1 
ae iceoo 2-308 2-368 2-442 2-473 2-520 2-552 2-586 


From the experimental results we may calculate a value of the ratio and thus by 
J interpolation in Table 3 we may obtain the appropriate value for m,/G?. From 
if this value we may calculate the variation predicted by the theory at the intermediate 
‘temperatures ; this is shown as the continuous line in Figure 2. There is seen 
‘to be close relationship between the experimental results and the predicted 
» theoretical variation. 

| (111) Correlation of Results 


From the values of 2,/G??, which were obtained from the conductivity 
) results, the thermoelectric powers at the various temperatures may be predicted, 
i and these are shown in relation to the experimental results in Figure 4, graphs II 
i and III. Good agreement is obtained. From a comparison of the values of 
| n,/G?? derived from the conductivity and the values of ,; derived from the Hall 
constant, values of G may be obtained, as shown in Table 4. The values of G 
. obtained are small, and it would appear that the various approximations in the 
theory have been compensated by giving the electron a small effective mass. 


Table 4 
Specimen I Specimen II 
n¢/G?/? from : 
Wondtctivityes melee ser 7-50'x 107° 3°50 1028 
Thermoelectric power .... TOL x 1020 SOO ange 
Flalliconstantir.eciciaiene ei. 1-69 x 10° IOSKey SMO 
Vea value olny sre eietets ers, 51s els 0-080 0-093 


(iv) Approximate Formula 
From the results which have so far been analysed we may conclude that 7* is 
very large for all temperatures used in the experiments. This indicates that the 
number of carriers is sufficiently large to produce these large values of n*. For 
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these and any lower temperatures (since * increases as the temperature 
decreases provided the concentration of free electrons is constant), we may 
therefore take the equations derived for the electrical properties, and use 
approximate expressions for the functions F’,(y*) which are valid only when 7* 
is large. 


R=tinevw rma ie bot (10) 
DIINO AI POP (11) 

16n2mk? /n\23 
ee 


These relationships are only applicable in ionic semiconductors in which nm, is 
very large and the temperatures are not too high. The linear variation predicted 
for the thermoelectric power by equation (11) is clearly followed by the 
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Figure 5. Graph of specific resistivity against absolute temperature for two specimens of CdO. 


experimental results, as may be seen from Figures 3 and 4. It will be seen 
that it is much easier to calculate the values of n,/G* from this simple formula 
than from equation (7). 

Figure 5 shows a graph of resistivity p against absolute temperature. The 
linearity predicted by equation (12) is quite pronounced.t At the higher 
temperatures, however, deviations should occur, due to the fact that 7* may no 
longer be sufficiently large for the approximations made to be valid. At lower 
temperatures the curve should again depart from linearity, and should approach 
a constant value as TJ approaches zero, due to the effect of the variation of the 
factor B at temperatures less than the Debye temperature. From the constant 
value approached as J-+0, a value for the Debye temperature could be 
obtained. 

+ The values of m¢/G*/ obtained from the conductivity curve will be seen to depend on the 


departure of Figure 5 from linearity. The values obtained will be approximate, and the accurate 
agreement shown in ‘Table 4 is probably largely fortuitous. 
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§9- CONCLUSIONS 
| As a further point of interest, the effect upon the predicted properties has 
( been considered when the temperature is assumed to be low. Here it is assumed 
“| that the temperature and concentration of free electrons is such that a vast 
‘| majority have energies less than Ay, unless they absorb a phonon of energy from 
| the lattice. In this case it has been shown by Shifrin (1944) that 


KE) =Ip{1 — f(E)}e*"(E/RT)'? 
which is very similar to that derived by Mott and Frohlich (1939) except 
that it contains a factor 1—/f(£) to allow for the cases in which degeneracy is 
/ appreciable. By substitution of this expression for the mean free path in 
®% equation (6) it follows that the thermoelectric power is given by 
k 5 Fiio(n*) + 3 Fs") 
= Ef ee 13 
eae 2 Fyj(q*) + 3F_10(4*) ay) 
The values of the thermoelectric power as predicted by this equation for various 
values of 7* are shown in Table 1 and Figure 6, beside those predicted by the high 


Thermoelectric Power (microvolts) 


Figure 6. Calculated values of the thermoelectric power for various values of 7*. 
————~— Values from formula (13), ————— Values from formula (7). 


temperature formula (7). It will be seen that the experimental results which 
have been obtained could not possibly have been explained by the low temperature 
formula (13) since as * approaches infinity the theoretical thermoelectric power 
predicted by this equation approaches its minimum value of 86 uv/deg. c., and 
over a large part of the range of temperatures over which measurements were 
made, values less than this were obtained. 

The agreement shown here between the experimental results, both of the 
author and of Andrews, and the theory developed appears to be very satisfactory, 
and may be taken as good evidence that the lines of argument leading to the 
theory are correct, and that the approximations made are applicable to cadmium 
oxide. The oxide is therefore a good example of a very degenerate electronic 
semiconductor possessing very metallic characteristics. 
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ABSTRACT. Films of some dielectrics become conductive under the influence of 
electron bombardment. This conductivity is investigated both for steady and for 
alternating bombardment in its dependence particularly on temperature, current density 
and the potential difference across the dielectric. The gain, i.e. the ratio of current 
through the film to the bombarding current increases strongly with temperature. 
On arsenic sulphide films, the gain reaches 40,000. The effect is discussed in the light of 
theories of conduction in dielectrics and of recent work on allied subjects. 


§1.- INTRODUCTION 

ONDUCTIVITY in a dielectric can be induced by photons only if 
their energy falls into one of the bands in the absorption spectrum of 

the dielectric; high-speed electrons, however, should induce conductivity 
in any dielectric, by the production of a trail of ions and slow electrons. Early 
experiments indicated that dielectrics became only very weakly conducting 
under the influence of electron bombardment, and it is only recently that a 
number of authors have drawn attention to the existence of a strong electron 
bombardment conductivity (or £.B.c.) (McKay 1948, Rittner 1948, Ansbacher 
and Ehrenberg 1948, 1949, Pensak 1949). We shall report here examples of 
strong E.B.c. which we have observed in thin films of dielectrics. 


* Now at H. H. Wills Physical Laboratory, Bristol. 
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| The accepted view of electronic conduction in dielectrics leads one to expect 
1 the conductivity to be greater the more perfect the crystal and the lower the 
+ temperature, because the number of trapping centres will be smaller in more 
perfect crystals and because the mean free path of the electrons decreases 
with increasing temperature. In fact, neither of these two conclusions is 
borne out by our experiments. A particularly strong effect was found in 
amorphous arsenic sulphide, As,S,, and the effect increases strongly with 
temperature. In §§ 2-4 of this paper the experimental evidence will be given; 
in §5 it will be compared with the results of other authors and its theoretical 
implication will be discussed. 


§2- THE E.B:C. OF DIELECTRIC FIEMS 
We have observed, with electrons up to 50kv., a strong E.B.c. effect if the 
electrons were directed against a film of dielectric sandwiched between two 
conducting electrodes (Figure 1; the galvanometer may be replaced by a 
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Figure 1. Figure 2. 


resistance coupled to a cathode follower). ‘The top layer was always thin enough 
not to retard the electrons appreciably. ‘The current across the film is 
composed of the leakage current 7, which is normally negligible, of zo, the fraction 
of the primary current penetrating the film which is observed when no field exists 
between the electrodes, and the £.B.c. current z,. ‘This was measured as a function 
of the following independent variables: the accelerating voltage V of the 
bombarding electrons, the current density J, or the current 7, of the bombarding 
electrons, the ‘bias’ voltage v across the film, the temperature 7 of the film, 
the thickness d of the film, and the variation with time of the bombarding current 
density. 
24-2 
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The bias voltage v is called positive if it accelerates electrons in the direction 
of flow of the incident electrons, i.e. if the bottom electrode is positive, and 4, 
is called positive if it corresponds to electrons moving from top to bottom. 
Whenever v is different from zero the sign of 7, agrees with that of v. The ratio 
i,/i, =g is called the gain, which we term positive or negative gain according to 
the sign of v. In most tests 7, was constant with time (D.C.-E.B.C.); in a few 
experiments 7, followed a square wave on-off pattern. 

The preparation of the specimens, i.e. of the dielectric films, which must be 
uniform and free from pin holes, depends on the individual dielectric. Films of 
aluminium oxide can be grown by anodizing solid aluminium; the specimens 
are completed by thermal evaporation of aluminium on to the oxide film. Films 
of mica are readily obtained by cleavage, films of glass by blowing thin-walled 
spheres. Specimens of arsenic sulphide and antimony sulphide are obtained by 
thermal evaporation of the sulphides on to previously metallized (gold plated) 
glass. Several crucibles at different angles to the base are used successively in 
order to minimize the risk of pin holes. The top electrode is an evaporated gold 
film and is made to overlap the bottom electrode over a small region only (about 
10 mm?) in order to minimize the risk of including faults in the film. 

Evaporated films of As,S, not only show a very high E.B.c. but also, owing to 
their great transparency, their thickness can be determined very readily by 
counting the contour fringes (the refractive index of these films was found by an 
interference method to be 2:38). These films are amorphous and have the same 
‘liquid’ Debye—Scherrer pattern as the B.D.H. ‘ Analar’ powder from which 
they are produced. ‘Their resistivity varies, between liquid air temperatures 
and 180°c., from about 10% to 10''ohmcm.; the variation corresponds to an 
activation energy of about 1-5ev. ‘The films have a breakdown strength of about 
106 volt/cm. (at room temperature). ‘Their dielectric constant is 10. 

The specimens were inserted into a chamber at the base of a demountable 
cathode-ray tube (Figure 2) so that they could be cooled to liquid air temperature 
or heated to whatever temperature they could withstand. The electrons were 
concentrated by a magnetic lens to a circle of about 0-3 cm. diameter on to the 
region of overlap of the electrodes. 


§3. RESULTS WITH STEADY BOMBARDMENT 


The critical voltage. Figures 3, 4, 5 show the variation of the gain g with 
bombarding voltage V for different values of the bias v of either sign at different 
temperatures. The gain is here, within experimental error, independent of the 
polarity of the bias; and g=0 for all values of v, if V is smaller than a critical 
voltage Vj. (The curves for v=0 represent 7, the fraction of the primary 
current, and not the gain; the gain at zero bias is always zero.) 

The existence of a critical bombarding voltage V) had already been observed 
by Becker (1904) and has also been confirmed by Pensak. V,? is proportional 
to d, the thickness of the film (Figure 6). This observation must be interpreted 
as a verification of the Thomson—Whiddington law, provided we assume that 
the E.B.c. begins when some primary electrons penetrate the whole film. Its 
thickness is equal to the range in the dielectric of electrons of energy V>. The 
straight line in Figure 6 passes through the origin. ‘This indicates that the loss 
of energy of the primary electrons in the top electrode is negligible. 
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V, is independent of the temperature of the specimens and of the bias, in 
agreement with the interpretation given above. 


Hence the Thomson- 
Whiddington law can be used to determine the thickness of films which cannot 
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easily be measured directly (for example, antimony sulphide films are opaque, 
aluminium oxide films when thick are milky and do not show interference 


fringes, and many samples are covered by a uniform oxide film and therefore 
do not show contour fringes). 
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The bombardment voltage for which the gain reaches a maximum. As the 
bombarding voltage is increased the gain at a particular bias and temperature 
rises to a maximum and then falls again. ‘This maximum voltage Vmax is 
independent of bias (Figures 3, 4, 5). Below room temperature, Vyyax 18 also 
independent of temperature, but for higher temperatures Vie _increases 
(Figure 7); it also increases with the thickness of the specimens. In Figure 8 an 
attempt is made to correlate Vy and Viy,x; values obtained from aluminium oxide 
films are included. The aluminium oxide films are usually thinner than the 
As,Sz films and consequently the values of Vp and Vix are much smaller, but 
some values fall well within the range of the values for As.53. All group around 
a common curve.. The error in Vix is large because the maximum is broad, 
but the relation V,,,,=2°3 Vo seems to hold. 

The variation of gain with bombarding voltage and bias. In general the (g, V) 
curves obtained from one specimen at different values of v are all of the same 
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Figure 7. Shift of V,,,, with temperature. As,S, film Figure 8. 
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shape and can be made to coincide by multiplying each by a factor f(v) depending 
on the value of v only. As the temperature is changed, the shape of the curves 
usually changes, as well as the values of f(v). Also above 100°c. for low 
bombarding voltages, the positive gain is, in general, slightly less than the 
negative gain. Over a wide range, the variations of gain with bias can be 
represented by 1/f=a(T)v+ b(T)z. 

At T=80°K., b is negligible, at T=26°c. (e.g. for a film 1-2 thick and for a 
bombarding current of 2-7 x 10-7 amp/cm?), b/a~ 0-013 if v is measured in volts ; 
at T=142°c. for a film 1-32 u thick and the same current, b/a~0-08. Thus the 
gain depends linearly on bias at liquid air temperature and becomes quadratic 
in bias at higher temperature. A saturation of gain with bias has not been 
observed. . 

A comparison of specimens of different thicknesses cannot be made at an 
arbitrary bombarding voltage. But as V,,,, is a function of the thickness alone 
(at least at and below room temperature) and the gain is also a function of the 
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‘ : t 
| thickness and of the bias, it is clear that by confining attention to the maximum 
gains alone, samples of different thicknesses can be compared. In Figure 9 the 
_t2tio gain/field at V,,,, is plotted against field v/d for different samples at room 
temperature. The values for samples with highest gain are grouped around a 
_ common curve; under favourable conditions g,,, depends on v/d only. The other 
walucs display no systematic grouping which could be related to the thickness d. 
One may attribute such loss in gain to effects of age of the specimens, previous 
treatment, and unfavourable conditions of evaporation of the As,S, films (for 
@ instance, some specimens may have been hotter during evaporation than others). 
: 0 The dependence of gain on temperature. The effect of temperature is 
_ complicated: a break in the behaviour occurs in the neighbourhood of 100°c.; 
‘below 100° c. the gain for given v,V and current density varies reversibly with 
temperature, that is, if the temperature returns to a certain value, the gain returns 
zw once to the value corresponding to this temperature. Also, the £.B.c. current 
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on reversal of bias is established instantaneously (within the period of the 

vanometer) and is of the same magnitude for positive and negative bias. 
- Above 100°c. the gain becomes dependent on the history of the sample. 
) At room temperature for a field of about 9 x 10° volt/cm. the maximum gain 
obtained on a new sample is of the order 50. After the sample has been heated 
to 100° c. or more, however, and allowed to cool to room temperature, a gain of 
several hundred is measured for the same field. (It is possible that the initial 
gain may be greater, as many adjustments, such as refocusing, are necessary 
before the first observation.) This value is steadily reduced on bombardment 
and finally reaches a value about 1-3 to 1-5 times the original value. Subsequent 
cooling to liquid air temperature with bombardment and warming up to room 
temperature usually re-establishes the original gain. 

For a particular bombarding current density and bias the gain rapidly 
increases with temperature, e.g. at V,,,, for arsenic sulphide films at 
I,=13 x 10-* amp/cm? and fields of the order 5 x 10° v/cm. the gain at liquid air 
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temperature is around 5, at 20° c. around 50, and at 150° c. around 500. Figure 10 
shows typical curves of log (gain) plotted against 1/7. It looks as if two processes 
were involved, one covering the range —200°c. to about room temperature, 
the other applicable above room temperature. The slope of the (log g, 1/T) curve 
is practically independent of the current density and is also independent of 
the thickness of the samples. 

The secondary E.B.C. effect at high temperatures. At elevated temperatures 
a phenomenon is observed which is reminiscent of the secondary photoelectric 
effect, and which will be called secondary £.B.c. effect. 

The E.B.C. current increases at first while the bombardment is going on until 
after several seconds or minutes a value (the ‘high’ gain) is reached which may be 
several times the value at the beginning (the ‘low’ gain). Once the high gain 
is established for a certain bias, the £.B.c. current follows instantaneously (within 
the period of the galvanometer) any variations of the primary beam. If now the 
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polarity of the bias is reversed, but the bombarding current is left on, the gain 
(in the reverse direction) is at first a small fraction of the previous high gain and 
rises slowly to about the same high value it had before the reversal of bias. The 
high gain is reached more quickly if the bombarding current is repeatedly 
interrupted, the bias being left on all the time. This ‘conditioning’ of the 
specimens does not take place, even slowly, without electron bombardment. 

The difference between low gain and high gain increases with increasing 
bias and increasing temperature, and decreases with increasing bombardment 
current density and with increasing thickness of sample; for example, for 
0-27 pa/cm? (the current density at which most measurements were carried out) 
there is no difference between low and high gain up to about 100°c., while for 
a current density of 0-0014 1a/cm? there is already a small difference at 20°C. — 

Plots of log (low gain) and of log (high gain)—dotted curve—against 1/T are 
shown in Figure 11. Over the temperature range the low gain appears to follow 
- a simple law; the high gain deviates from this at 120°c. 
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The low gain values are somewhat uncertain, because they depend on the 
oreceding conditioning of the specimen for high gain before reversal of bias. 
»/This conditioning may take several minutes for large current densities and for 
“temperatures only slightly above 100°c 
‘| Below 100°c. the gain for positive ama negative bias is usually the same (the 
sjother variables being the same) but above 100°c. the negative high gain exceeds 
“athe positive high gain by as much as 100% at low bombarding voltages, and 
‘gapproaches the negative gain for large V. In some specimens the negative gain 
exceeds the positive gain over the whole range of V. 
| The highest gain ever obtained was in the neighbourhood of 40,000 for a 
‘specimen 2:04 microns thick at 163°c. The bombarding voltage was 24kv., 
‘athe bias 100 v. and the current density 6:8 x 10-" amp/cm? (Figure 12 (a)). This 
“current density was already so low that no increase in gain could be expected by 
‘further reduction. On the other hand, at higher temperatures the leakage current 
would have masked the E.B.c. current. It should be emphasized that the 
)secondary E.B.cC. current is the current under electron bombardment after 
“conditioning of the specimen whereas the secondary photoelectric current is a 
‘current which continues after the illumination has ceased. 
The dependence of V,,,, on temperature. Figure 7 shows that V,,,, increases 
i with increasing temperature. - This displacement becomes noticeable at 60°c. 
After a specimen has been heated and then cooled again to room temperature, 
1 Vinax is very often found to be less than before heating. The original values of 
# gain and V,,,, are restored after the specimen has been bombarded at liquid air 
“temperature. The changes in V’,,,,, and the changes in gain referred to previously 

appear to be closely related. 

The dependence of gain on current density. The gain is independent of current 
) density only at very low temperatures. Figures 12 (a), (b) show log (gain) plotted 
against log (current density) at two temperatures for various values of the bias. 
The gain decreases with increasing current density and is therefore sensitive 
) to focusing ; 1t becomes independent of current density only below a critical value 
7 which is smaller at higher temperatures and for thinner specimens. ‘The curves 
j indicate gocI,-*”; «(T) is very small for low T. 
In Figure 12 (4) values both for log (high gain), full curve, and log (low gain), 
| dotted curve, are given. For large current densities the difference between 
the two becomes small. Of all £.B.c. effects only the critical bombarding voltage 
is, within experimental error, independent of temperature. 

Charges in bombarded films. At room temperature, for small bombarding 
voltages (just above V,) the following ‘ discharge’ phenomenon is observed. 
If the bombardment current is switched off before the bias is reduced to zero 
and then switched on again with the bias at zero, a current is observed in the 
direction opposite to that of the previous E.B.C. current, rapidly decreasing to the 
normal i. The time integral of this current, as recorded by a ballistic 
galvanometer was found to be 8x10-® coulomb for a bombarded area of 
0-058 cm?, for a specimen 2 thick. If the bias is reduced to zero before switching 
off the primary current, no ‘discharge’ current is observed on subsequently 
switching on z,. ‘The magnitude of the ‘discharge’ current decreases to zero 
as the bombarding voltage approaches V,,,.. 
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From the ‘discharge’ current at low V the number of traps in the film is 
estimated to be not less than 3 x 10'* cm~? under the assumption that on discharge : 
the carriers traverse on the average half the film. 

Experiments with films other than As,S,. Antimony sulphide films can be | 
prepared in the same way as As,S, films. Thin films are deep red by transmitted 
light and micro-crystalline. Our experiments indicate that Sb,S, shows the same 
effect as As,S; and under corresponding conditions probably a larger gain. 
Al,O, films grown on Al show also a strong £.B.c. depending on the method of © 
anodizing and the type of Al used. Most of our experiments were carried out on 
films produced by the chromic acid process. The top electrode was evaporated 
aluminium. ‘These specimens have the remarkable property of changing their 
capacity reversibly to one-twentieth of their value in air when placed in a vacuum. 
Their thicknesses were obtained from E.B.c. data as mentioned above. Curves 
corresponding to Figure 3 are shown in Figure 13. They differ from the As,.S, - 
curves in the following respects: (a) the gains are smaller; (5) the factors f(v) are 
up to 30% different for positive and negative bias; (c) the negative gain seems 
to saturate with temperature (Figure 14) : the straight portions of a (log (gain), 1/T) 
curve correspond to an activation energy of 0-35 ev. for positive gains and 0-5 ev. 
for negative gains; (d) no secondary E.B.c. effect and no variation of V,,,, with 
temperature is observed; (e) the specimens show for current densities greater 
than 3 x 10-’amp/cm? a fatigue effect at high temperature: the gain decreases 
with time under bombardment. 


§4. EXPERIMENTS ON As.S,; AND Al,O,; SPECIMENS USING 
MODULATED BOMBARDING CURRENTS AND STEADY 
BIAS VOLTAGES (A.C. METHOD) 

Experimental arrangement. ‘The electron beam is modulated by driving the 
grid of the electron gun with a square-wave generator at from 22 c/s. to 66 ke/s. 
During one half of the cycle the bombarding current is cut off, and during the 
other half rises to a constant value within a small fraction of the duration of the 
total current pulse, so that it varies with time as shown in Figure 15 (a). In order 
to measure the E.B.C, current (A.c./E.B.C.) for this alternating bombardment the 
galvanometer of Figure 1 is replaced by a resistance X and the voltage developed 
across X is, after amplification, applied to the screen of a cathode-ray oscillograph. 
The amplifier had a gain of about 20,000 independent of frequency up to 0-5 Mc/s. 
The a.c./£.B.c. deflection is always compared with the deflection obtained by 
passing the bombarding current directly through X (switch in position 1, 
Figure 1). In addition the average primary and secondary currents are measured 
by a galvanometer. For comparison, a measurement of the D.c. gain is also made 
for a D.c. bombarding current equal to the average a.c. bombarding current. 

The scope of the measurements is similar to that described for the D.c. 
experiments; it includes the variation of gain with bombarding voltage, 
bombarding current density, applied bias, temperature of the specimens and 
frequency of the bombarding current pulse. It includes also a comparison 
between the wave-form of the £.B.c. current with that of the bombarding current. 
As these are, in general, not identical the definition of gain becomes somewhat 
arbitrary. In the following, a.c. gain denotes the ratio of the amplitudes measured 
from peak to peak, and average a.c. gain the ratio of average secondary to average 
primary current. The capacity of the specimens is so large that it is necessary 
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‘to allow for distortion of the signals by the measuring circuit. During the 
‘bombardment the film becomes conducting, so that the sample is equivalent to 
ja condenser of constant capacity C shunted by a resistance R(t) which varies 
periodically with time. 
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Figure 15 (a)-(i). A.c./E.B.c. (a) represents the bombarding current varying with time ; 
(b)-(z) E.B.c. signal from a As,S3 film 1:9 w thick, Iy=13-6 10-8 amp/cm?, pulse length 
0:5 msec. The numbers indicate the a.c. gain; (b), (c) specimen at 26°C. small and large 
bias, (d), (e) specimen at 140° c., small and large bias, (f)-(z) specimen at 120° c., change of 
appearance of cathode-ray oscillograph trail on reversal of bias: (f) high a.c./E.B.c., with 
positive bias, (g) bias just reversed, (h) negative bias on for 1 second (low A.c./E.B.C.), 
(i) negative bias on for some time (high A.c./E.B.C.). 


Figure 15 (k)-(p). a.c./E.B.c. signal from aluminium oxide film at 300° c., pulse length 5 msec. 
The numbers represent the a.c./z.B.c. gain in (J) and (n) (V=20 volts), in (k), (m), (0), (p) 
(v=0) they give the portions of the primary current reaching the bottom electrode. 
Ip=0-27 pa/om?, (0) and (p) Ip=2-7 walem? ; (n), (I), (0) V=5 kv. ; (m), (n), (p) V=10 kv. 
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Let 1/R be represented by a Fourier series: 
| NOMS We ars 
a ; reot); yitanaey ey teeeehe 1 
Re + sash a, exp (jrwt), (1) 


where &’ means that r=0 is excluded from the summation. A simple 
calculation leads to an expression for the E.B.c. current 


+ 00 
i,=v/Rot+v X' a, 1+(rwCX)*}- exp j(rwt—p,) see ae (2) 


where ¢, =tan-!(rwCX). 

The first term in (2) is the average current. Comparison of (1) and (2) shows 
that the current 7, will follow the time variation of 1/R faithfully as long as 
rwCX <1 for all harmonics of large amplitude; an estimate of the highest 
frequency faithfully reproduced may be obtained from the first harmonic by 
setting wCX=1. Owing to the comparatively low gain of the amplifier a value ~ 
of X equal to 10,000 ohms was used; the capacity of an average As,S, specimen 
is500cm. Hence, the limiting frequency is 30,000 c/s. For the highest frequency 
used the distortion is appreciable and the amplitude measured is too small by a- 
factor (1+2?)"?. For v=2,000c/s. the distortion can be neglected. 


Results 


The wave form of the A.C./E.B.C. current. Figure 15(a) shows the 
wave-form of the primary current. The a.c./E.B.c. signal for As,S; specimens 
is shown in Figures 15(4)-(e). It appears that part of the E.B.c. current responds 
instantaneously (within the limits of observation) to variations of the primary 
current, but that a considerable fraction of the total a.c./E.B.c. current takes time 
to establish itself. ‘This is also borne out by Figure 17 discussed below. With 
increasing temperature, this time-dependent fraction of the E.B.c. current 
increases. 

An interesting behaviour (Figures 15 (f)-(z)) is observed on sudden reversal 
of bias at a temperature sufficiently high to show the difference between high gain 
and low gain under D.c. conditions. It is seen that before the state of low gain is 
reached, the specimen passes through a stage when the gain is nearly as large 
as before reversal of bias, but the wave-form is almost that of the primary current; 
the time-dependent fraction of the £.B.c. seems to be absent. This state of affairs 
may correspond to the space charge neutralization, discussed by McKay. The 
A.C./E.B.C. Current then decreases to a minimum (Figure 15(h)) with both 
instantaneous and time-dependent portions of the current present, and finally 
rises to a steady high value (Figure 15(z)) when the wave-form is substantially 
that observed before reversal of bias. 

At low bombarding voltages the ‘ discharge’ phenomenon already observed. 
under D.c. conditions is observed directly on the cathode-ray oscillograph screen. 
When the primary current is switched off, then the bias reduced to zero, then the 
(modulated) primary current switched on, a large a.c./E.B.c. pulse is observed. 
in the reverse direction, which rapidly reduces to the normal fraction of the 
primary current reaching the bottom electrode. This shows that a large 
polarization exists in the insulator under bombardment and application of bias, 
which remains after the bombardment has ceased. At higher bombarding 
voltages this reversal is not observed. 


Electron-Bombardment Conductivity of Dielectric Films 373 


The behaviour of the Al,O; specimens resembles closely that of As,S, 
samples, but on account of their large capacities, the distortion of the wave-form 
of the A.c./E.B.c. current is already noticeable at low frequencies. Owing to the 
greater strength of the Al,O, films very large primary currents could be used. 
igure 15 (0) shows clearly the build-up of internal space charge at zero bias and 


‘to the ones shown by McKay (1948). As the bombarding voltage is raised the 
equilibrium is established almost instantaneously (Figure 15(p)). 

The dependence of A.C./E.B.C. on bombarding voltage and current density for 
yng AsyS3. Figures 16(a), (b) show, for two temperatures the variation of a.c. gain 
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value of V, (coinciding with the value of Vy for D.c. experiments), there exists 
a Viyax and the negative and positive gains are equal. ‘The a.c. gain is, however, 
always smaller than the average a.c. gain; this indicates that some E.B.Cc. current 
flows during the off period of the primary current. The a.c. gain does decrease 
with increasing current density, but not as rapidly as the D.c. gain. 

The dependence of A.C. gain on the frequency of the primary current. 
Figures 17 (a), (b) show the variation of a.c. gain at two temperatures with the 
frequency of the primary beam. ‘The gain is initially constant, then decreases 

‘with increasing frequency, approaching a constant value at about 500c/s. Above 
10,000c/s. the gain drops rapidly. ‘The second decrease must be attributed 
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to the circuit. The D.c. gains and the average a.c. gains are also shown on the 
diagrams. The average a.c. gain is, throughout the frequency range, equal to 
the D.c, gain. . 
The decrease of gain with frequency and the shape of the a.c./E.B.c. current 
pulse can be interpreted as follows. For primary current of unit amplitude 
let I, be the time-dependent fraction of the a.c./E.B.c. current (Figure 18 (a)) 
and let J, be the instantaneous fraction (Figure 18 (b)), so that the total a.c./E.B.c. 
has the form shown in Figure 18(c). It will be assumed that the rise and fall of 
the time-dependent current follows an exponential law with a relaxation time 7, 
so that for the rise dJ,/dt=y/r—JI,/7, where y/r represents the effect of 
bombardment, and for the fall d/,/dt=—J,/r. ‘Then for the steady periodic 
state (Figure 18(c)) at a particular bias and bombardment voltage but for a 
frequency v=1/2T the maxima B and the minima A of J, are related by 


Bexp(—T/r)=A and B=y-—(y—A)exp(—T/r). ...... (3) 


The average A.c. gain as measured by the galvanometer is the average height © 
of curve (c) and is given by 


G-[cr+|_ {y—(y—A) exp(—t/n)hdt+ | Bexp(—t/rdt | /7=C-+y, 
be (4) 


independent of frequency. The a.c. gain is the difference of the maxima and 
the minima of curve (c): 


G=C+B-A=C+ytanh(T/27),  —...... (5) 
Hence for very high frequency (<7) G=C, and for very low frequency 


(T>7r) G=C+y=G, so that G drops by y between low and high frequencies. 


Field 


Anode 
Figure 18. Figure 19. 


Also, for T>7, I simply reproduces a sequence of equally long on and off periods 
so that then, neglecting current density effects, G=G=g. This is confirmed 
by Figure 17. 

Half way between the low and the high frequency regions G=C+y/2. The 
frequency for which G=C++y/2 can be read from Figure 17 so that 7 can be 
found from (5): 


C+y/2=C+ytanh(1/4v'r) eee (6) 
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Mt ir 1/7 =2v'ln3.. From Figure 17 for the 75-volt bias curves 
at 26°c., G=10-5, C=6-8, v’=30sec", 1/7=66 sect, 
at 125°c., G=365 C=125, v'=27sec, 1/r=54 sec). 


With these values of G and C and an average 1/7=60 sec"! equation (5) 
\"teproduces in fact the full curves of Figure 17 below 10*c/s. This shows 


These results indicate that the time-dependent part of the E.B.c. depends 
mm an intermediate excitation. One would expect, however, the relaxation 


following brief discussion, this intermediate process will be left out of account. 


§5.. DISCUSSION OF RESULTS 


The main net effect of the passage of a fast electron through an insulator is the 
}2xcitation and ionization of atoms along its path. Some of the secondary electrons 
“land holes will recombine at once, others will wander through the crystal until 
sthey recombine with opposite charges or are trapped or reach the electrodes. 
iLf no charges are lost by recombination, the total transport of charge through the 
icrystal per ion pair produced is always e, and the gain, as defined above, will be 
Jequal to the number of pairs produced per primary. One may expect that an 
energy of about 5 ev. is required for the production of a pair. Then the number 
Jof pairs produced by an electron of 20 kv. will be 4,000, that is, far higher than the 
igain observed at low temperature and far less than the gain observed at high 
temperature. 

_ But trapping and recombination are well illustrated by the existence of the 
‘critical voltage V,: below V, the pairs produced recombine near their birthplace 
‘because the charges carried beyond the region of ionization are trapped and set 
jup a space charge which prevents any further movement of carriers. Also, the 
discharge current observed must be interpreted as the release by the primaries 
‘of trapped charges which move through the film in such a way that it becomes 
| unpolarized. ; 

|. A general equation between the electric field, the density of electrons and 
) holes, their mobilities, the density of traps and trapped charges and the rates of 
{ionization and recombination is readily established but intractable. ‘The far 
| simpler and closely analogous problem, however, of the ionization chamber in 
) which there are no traps and for which the ionization is constant through space 
§ can be solved by successive approximation (Mie 1904, Thomson 1906), and this 
) solution should still indicate at least qualitatively the electric field in the dielectric 


i} 


/ under £.B.c. 

] J. J. Thomson and G. Mie have shown that the field in the ionized gas is 
) nearly constant and small throughout a large part of the gas and rises 
/ exponentially near the electrodes (Figure 19). In the region of nearly constant 
| field the number of recombinations of charges almost equals the number of charges 
| produced by ionization. The narrow regions of strong field near the electrodes 
_ contribute, therefore, nearly all the current. As the applied field increases, the 
region of small constant electric field shrinks, until, in the limit, the field 
_ throughout the sample is nearly constant and equal to the applied field. When 
that happens all the charges produced in the gas are drawn to the electrodes. 
The strong fields near the electrodes arise from the separation by the applied field 
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of charge carriers of opposite sign. When the carriers have equal mobility, the 
field distribution is symmetrical; for unequal mobilities the field is highest at the 
electrode which attracts the slower carriers. Figure 19 should be representative 
for the field in dielectric films under £.B.c. 

Newton (1949) has suggested an expression for the average field Zp required 
to suppress recombination over a given layer of thickness D. By the use of this 
relation one can arrive at values for the field both in the central region and near 
the electrodes. 

Approximate values for the rate of production of pairs and the coefficient of 
recombination are first obtained as follows: For a primary current density of 
13-6 x 10-Samp/cm? the number of electrons impinging on the film equals 
86 x 10!’ electrons/cm?sec. Let us assume that the electrons ionize uniformly 
throughout the film and completely lose their energy within the film thickness 
d=3x10-*cm., also that the ionization energy per pair is 5ev. and the © 
bombarding voltage is 20kv. Then the rate of production of pairs is 


B =(86 x 102 x 20 x 10)/(5 x 3 x 10-4)~12 x 1038 cm-* sect, 


The rate of recombination will be given by «n,n_ where n,,n_ are the densities 
of positive and negative carriers and « is a factor approximately equal to the 
product of the velocity of the carriers and an atomic cross section, 
a~6x 107 x 10-45 =6 x 10®cm3sec!, corresponding to an_ electron of | 
Ein =0-5 ev. One also requires values for the mobilities u,,u_ of the carriers. 
Usually one reckons these to be in the neighbourhood of 100cm?sec™ volt-t. 
This value refers to carriers in the free state. Here, in view of the presence of 
traps, it should be multiplied by the ratio of the time the carriers are free to the 
time they are trapped. We cannot estimate this ratio. The calculation will 
first be carried out for u,<u_, u,~1cm?sec!volt"!. ‘Then the time required 
for either charge to traverse the layer is 7,,_<D/u, Ep. During this time 
BDz =BD?/u,,_ Ep pairs are produced in the layer. The density of free charges 
approximately equals the number of charges produced in the time 7 divided by 
the thickness of the layer, n,, ~D/u, E) and the number of recombinations 
per second in the layer will be smaller than an,?D~a«($D)?/D(u,Ep)?. 
On the other hand, the number of electrons produced in the layer per second 
is BD. The number of recombinations is negligible if it is much smaller than the 
number of electrons produced, that is if Ep) ~($a)!?D/u,. Hence, even in order 
to suppress recombination over the whole thickness d of the specimen 
(3 x10-4cm.), a field of only 200v/cm. would be required. But in all our 
experiments the applied field exceeded 200 v/cm. by a large factor (10 to 10°). 
If, therefore, the applied field existed throughout the film all liberated electrons 
and holes should have been drawn to the electrodes and the largest gain (in the 
example, 4,000) obtained, instead of the small fraction of this usually found. 
Hence one must conclude either that in a wide central region the field is far 
smaller than 200v/cm., and very large fields at the boundary make up the total 
potential drop across the specimen, or that the value taken for u is far too high. 
Now, the thickness A of these boundary layers and the fields F, there can also be 
estimated : for the bombarding voltage and the current density of the example anda 
bias of 100 v. a gain of 20, i.e. an E.B.c. current 7, of 1-7 x 1018 electrons/cm? sec. is 
observed. If this comes entirely from the boundary layers, 7,=2BA or 
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In A~7x10-?cm. As further, approximately 2AF, + (d—2A)E,=100 volts, and the 
‘field E, in the central region must be smaller than 200 volts/cm., 


E,~[100 — (3 x 10-4—1-5 x 10-®)200]/1-5 x 10-6 6 x 10? v/em. 


“This value for E, is rather high although the macroscopic breakdown field 
“strength can be exceeded in very thin layers (see, for example, Plessner 1948). 
“| But this high field would entail field emission from the electrodes which would 

{invalidate the relation A=7,/28. Besides, field emission can barely play an 


i ,u<1 for the condition under examination; in order to arrive at a more reasonable 
"\S value for E,, say 3 x 10° v/cm., u would have to be about 6 x 10-*cm? sec volt 4 
“| and E,~3 x 10° voltcm™ so that most of the potential drop takes place in the 
‘9 central region. 

But the mobilities as determined by the rate of release of charges from traps 
7 must increase exponentially with temperature and it appears that the resultant 
1 increase in gain is due both to the increasing width of the boundary layers and 
) also to the increased boundary fields. Even if the entry of carriers does not take 
“ place at low gain it should do so when a high gain is observed. The main evidence 
0) for this comes from the high gains observed at high temperatures and high bias 
) which cannot be explained by the formation of ion pairs by the incident beam 
l) alone, especially as even at the highest gains observed no trace of saturation was 
/ found. 

. Even higher gains than those reported here have been observed in cadmium 
sulphide crystals irradiated with light and with «- and f-particles. These crystals 
were first prepared by Frerichs (1946) and have been investigated in detail by 
Frerichs and Warminsky (1946), Kallmann and Warminsky (1948), Fassbender 
(1949) and Frerichs (1949). In particular, Kallmann and Warminsky have 
reported secondary currents equivalent to up to 10,000 times the number of 
ion pairs produced, and both these authors and Frerichs have concluded that 
this effect must be due to the entry of carriers. ‘Iwo others processes are known 
which might lead to extra secondary currents, viz. ionization by collisions and 
“an excess temperature of the electron gas over the lattice temperature (Froéhlich 
and O’Dwyer 1950). The first possibility has also been discussed by Kallmann 
and Warminsky and rejected on account of the comparatively low fields used. 
In the experiments reported here one can exclude ionization by cellisions because 
of the too gradual increase of gain with bias. Again, from the figures given by 
Frohlich and O’ Dwyer it appears that an excess electronic temperature would not 
be sufficient to explain the high gain observed or extrapolated. Hence, one should 
conclude that the highest gains observed in E.B.c. are due to the entry of carriers. 
In many respects, however, the effects in CdS and those observed here are very 
different. 

It is not yet possible to follow up in detail the related processes of formation 
of ion pairs, establishment of space charge and the entry of carriers. We know 
very little about the distribution of ionization in the dielectric, or about the 
variation of effective mobility with temperature and we do not know how the 
mobility and current density of the carriers affect the shape of the field. One 
would expect, however, that in the low gain region—where no entry of carriers 
takes place and which corresponds to the regions investigated by McKay and 
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Pensak—the gain roughly follows the effective mobility which should increase 
exponentially with temperature as it depends mainly on the activation of traps. 
Such traps have been studied extensively by Randall and Wilkins (1945) and 
Garlick and Wilkins (1945). In fact, Figure 10 suggests that there are two types 
of trap present: between liquid air and room temperature thermal release 
only takes place from shallow traps less than 0-01 volt deep; the rapid rise of gain 
at high temperatures is due to thermal release from deeper traps with an 
activation energy of 0-5 ev. 

The variation of gain with current density is a complicated effect which 
could be explained in a number of ways. A displacement of V,,, with 
temperature is not incompatible with Thomson—Whiddington’s law. In view of 
the complexity of the process the optimum distribution of ionization throughout 
the layer will be sensitive to changes in mobility. It appears attractive to associate 
the secondary effect with the entry of carriers into the dielectric in the same way 
as Mott and Gurney (1949) have interpreted the secondary photoelectric effect. 
Also, the slow awakening of the secondary E.B.c. is analogous to the slow 
development of the passage of carriers observed in cadmium sulphide crystals. 
But the explanation given by Mott and Gurney does not provide for the 
instantaneous response of the E.B.c., once the dielectric films are conditioned 
by previous bombardment. 


APPENDAEX 


THE PHOTOCONDUCTIVITY, OF As; 5S. 0h TEMS 


The films show photoconductivity. The currents are of the order of 
10-®amp. at room temperature for an illuminated area of 0-1cm? and fields 
of about 0-5 x 10%v/cm. when a 75-watt lamp is placed at 1cm. away from 
specimen. 

The photocurrent, similar to the E.B.c. current, increases rapidly with 
increasing bias and with increasing temperature. It follows the modulations of 
the light source (ON, OFF) without a time lag within the period of the galvanometer 
below 100° c. always and at higher temperatures if the sample has, for some time, 
been illuminated under steady bias. But for a virgin sample, at a given bias, the 
photocurrent reaches a maximum slowly, and on reversal of the bias starts at a 
much lower value and rises slowly to the value it had previously. 

This behaviour differs from the secondary photoelectric current which decays 
or increases slowly when the light is switched on or off. 
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The Specific Heat of Liquid Helium at Temperatures 
between 0°6° and 1°6°kK 


By R. A. HULL,* K. R. WILKINSON anp J. WILKS 
Clarendon Laboratory, Oxford 


Communicated by F. E. Simon; MS. received 2nd October 1950 


ABSTRACT. A detailed account is given of experiments made to measure the specific 
heat of liquid helium between 0-6 and 1:6° kK. Below 1:4° the specific heat may be expressed 
as C=0-024T*? cal/gm/deg. Some implications of this result are briefly discussed. 


§1. INTRODUCTION 

HE specific heat of liquid helium below 1°x. is of some consequence in 
| the interpretation of the properties of liquid helium and was first measured 
by Pickard and Simon (1939) using a sealed capsule containing 
paramagnetic salt and helium. ‘This work was never fully published as the 
authors subsequently had reason to doubt its accuracy. Keesom and Westmijze 
(1941) announced that they had made similar measurements but their results 
have likewise never been fully published. The brief notices that have appeared 
show that while Keesom and Westmijze found that the specific heat varied as T® 
right down to 0-6°, Pickard and Simon found a T® region below 0-8°; there is 
also a discrepancy in the absolute magnitudes at 1:0°. ‘This paper describes a 
new attempt to measure the specific heat and clear up the somewhat confused. 
position. As no full account of measurements with a capsule has yet been 
published it may be useful to describe the experiments in some detail. The 
specific heat with which we are concerned is the specific heat of the liquid in 
equilibrium with the saturated vapour, though for brevity it will be referred to 

as the specific heat of helium. 


* Dr. Hull was killed in an Alpine accident on 23rd August 1949. 
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§2. THE EXPERIMENT 


A calorimeter containing liquid helium and some paramagnetic salt is 
demagnetized adiabatically from a starting temperature of about 1:0°K. As the 
salt cools, the thermal conductivity of the helium ensures that the calorimeter 
and all its contents cool down with it.* The calorimeter is then heated at a 
known rate and its resulting change in temperature is observed. Provided that 
the thermal capacities of the paramagnetic salt and the calorimeter are known, 
the specific heat of the helium may then be obtained. Alternatively we may get 
the same result by performing experiments in the same calorimeter with a constant 
quantity of salt but with different amounts of helium. ‘The magnetic moment 
of the salt ballistically is measured, and this determines the temperature 
provided the dependence of its susceptibility on temperature is known below 1°. 
As both the specific heat and the susceptibility of iron ammonium alum have 
been measured at these temperatures by Kurti, Simon and Squire (Cooke 1949) 
we used it in all the experiments. The calorimeter was mounted inside a double 
walled cryostat of a Simon expansion liquefier; a general description of this, 
and of the technique of adiabatic demagnetization, has already been given by one 
of us (Hull 1947). 


§3. THE CALORIMETER 


The size of the calorimeter is limited by the requirement that the whole 
apparatus must fit between the poles of a magnet. Thus as the specific heats of 
both helium and salt are small in the region below 1°, the thermal capacity of the 
calorimeter and its contents will also be small; with the calorimeter finally used its 
minimum value at 0-6° was only about 4 x 10*erg/deg. Heat wassupplied atthe rate 
of about 1000 erg/minute so as to raise the temperature 0-1° in about five minutes, 
and the thermal capacity obtained by differentiating a plot of temperature against 
time. ‘The heat leak due to imperfect insulation should therefore be no greater 
than about 100 erg/minute. 

The simplest form of a calorimeter would have been a thin walled vessel 
fitting inside a double walled helium cryostat and communicating with the rest 
of the apparatus by a thin walled German silver tube passing through the 
cryostat. However, the existence of the helium II film leads to a heat inflow 
(Rollin and Simon 1939) which even with a fine capillary may amount to several 
hundred erg/minute. 

Our measurements were made using the sealed-off ‘capsule’ technique first 
used by Kurti, Rollin and Simon (1936). Helium gas at about 120 atmospheres 
pressure together with the paramagnetic salt is introduced at room temperature 
into a thick walled calorzmeter or capsule which is then completely sealed off 
with soft solder. ‘The capsule is mounted on nylon threads in a cradle (Hull 1947) 
and placed in the cryostat (Figure 1) so that the only thermal contact with the 
cryostat is via the threads. On cooling the capsule the helium gas condenses 
producing about 1:5cm of liquid. Using this technique we have consistently 
obtained heat inflows of only about 70 erg/minute. 


* Although the heat conductivity of liquid helium falls off below 1-0° (Kurti and Simon 1938 a) 


it should be adequate for our purpose down to at least 0:2°. Certainly there was almost no sign of 


any overheating occurring. 
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St ee CAPS UIE E 


In order that there is as much helium to work with as possible the capsule 
(Figure 2) is made as large as will conveniently fit into the cryostat. It is machined 
from a rod of phosphor bronze (8% tin) which has a high tensile strength, a low 
electronic specific heat, and is free of superconductivity and magnetic effects at 
helium temperatures. 

One end is closed by a screwed plug which may be removed in order to intro- 
duce the salt in the form of a powder C which is pressed down into the bottom 
of the capsule. The plug is then screwed and soldered in while the other end 
of the capsule is kept cool in a mixture of solid carbon dioxide and acetone so that 
the salt is not decomposed by the heat. Helium gas is introduced into the capsule 
through a well-cleaned copper capillary tube A soft soldered into a hole in the 
screw plug. In the bore of this tube is a wire of soft solder B and the gas is sealed 
in the capsule by hammering the capillary flat and then applying heat so that the 
solder runs. Finally, the capillary is cut above the solder seal. 


7 
43c 
lo8cem 
24 
mm 
platinum - 
glass seal for 
heater wires 
h—+| 
(@) (6) 16-5 mm. 
Figure 1. Figure 2. 


After a set of experiments the quantity of helium in a capsule is determined 
by making a fine cut in the copper capillary tube with a razor blade, collecting 
the gas over water and measuring its volume. An independent result is also 
obtained by weighing the capsule before and after releasing the gas. The 
measurements of mass and volume were consistent to 1%. 

An electrical heater of No. 40 s.w.c. silk-covered manganin wire was wound 
on to the capsule, the whole being coated with ‘ Demarda’ bakelite varnish 
(supplied by Messrs. Bakelite Ltd.) and baked so as to obtain adequate thermal 
contact with the capsule. This was important as, if the heater or the leads to it 
were overheated, any gas adsorbed on them might be evolved and ruin the 
vacuum insulation round the capsule. However, no such trouble was encountered. 
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§5. THE MAGNETIC PROPERTIES OF SOME 
TECHNICAL MATERIALS 

All our temperature measurements depend on accurate determinations of the 
susceptibility of iron ammonium alum which, like all paramagnetic susceptibilities, 
is small even at 4°. Therefore it is important that no other part of the apparatus 
should be magnetic and give effects which are superposed on those due to the 
salt. The parts of the apparatus most likely to give trouble in this respect are the 
capsule and the cryostat which both surround the salt, because a very small 
amount of ferromagnetic impurity in this mass (about 200 gm.) of metal would 
give effects comparable with those produced by our 6 gm. or so of salt. Undesirable 
magnetic effects were, in fact, found when using German silver cryostats and 
constantan resistance wire. 

The reasons for this are not clear although the work of Kaufmann and Storr 
(1943) has shown that copper-nickel alloys can have a susceptibility which is 
both temperature and field dependent. Both German silver and constantan 
have a high nickel content so we turned to materials containing little or no nickel. 
For the cryostat and capsule we used phosphor bronze and for the electrical 
heater and leads manganin wire (about 5% nickel). Both these alloys have been 
entirely satisfactory* although the bronze is somewhat difficult to machine. 


§6. THE RATE OF EVAPORATION OF HELIUM 
FROM THE CRYOSTAT 

The cryostat was suspended below the expansion chamber of the helium 
liquefier by a thin walled tube of German silver which served both to introduce the 
liquid into the cryostat and as a pumping line by which the temperature of the 
helium bath could be reduced. Theevaporation to be expected may, from known 
data, be estimated at about 6 cm? of gas (N.T.P.) per minute. This is made up of 
2 cm® evaporated by heat coming down the tubes supporting the cryostat, 3 cm* 
which flows up the tube as a Hel] film below the A-point, the remainder being due 
to any residual gas in the vacuum insulation and in very small measure to radiation. 
These expectations were only borne out when the temperature of the cryostat was 
above the A-point. As the helium in the cryostat was pumped down from the 
-point to 1-5° the evaporation rose to as much as 20 cm3/min. On the occasions 
when a second tube was introduced into the cryostat to measure the vapour 
pressure above the liquid, the evaporation was greater by a factor approximately 
equal to the increase in the total perimeter of the tubes up which film could creep. 
Moreover the dependence of the amount of evaporation on the temperature of the 
cryostat took the same form as would be expected if it were all due to film flow. 
Thus it appeared that film flow of a greater magnitude than that usually reported 
was taking place. Subsequently, Bowers and Mendelssohn (1949) have shown 
that such an increased flow may be produced if the helium gas condensed into the 
cryostat contains impurities: such impurities are deposited as a thin layer on the 
inner walls of the cryostat and the tubes to it, so producing an increased flow of 
film. As confirmation of this we may mention that each time we condensed 
helium gas into the cryostat after it had been pumped down we observed a 
considerable increase in the evaporation. 


_ .* There does not appear to be a simple test which can be made at room temperature to give an 
indication of the low temperature behaviour. For example the constantan wire which showed 
much remanence was in no way attracted by a powerful permanent magnet at room temperature 
Maiti! the manganin wire whose low temperature performance was satisfactory was slightly 
attracted. 
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§7.. THE THERMAL INSULATION OF THE SPECIMEN 

To obtain accurate measurements, the heat leak to the calorimeter must be 
kept to a minimum, i.e. once the calorimeter has been cooled down, it must stay 
cold and not receive any appreciable heat from its surroundings. The problem 
of thermally isolating specimens at very low temperatures has been very fully 
treated by Cooke and Hull (1942), who showed that for good thermal insulation it 
is necessary to support the specimen on thin fibres and to have a minimum of 
residual gas* around it before demagnetization. This latter condition can be 
fulfilled by using a baking out technique (Hull 1947). One further precaution was 
necessary because of the very peculiar effect that any vibration of the specimen was 
accompanied by an evolution of heat, sometimes as great as 300 erg/minute. This 
effect is still being investigated, but in the present experiments the vibration was 
reduced by mounting the apparatus as rigidly as possible. 

Our final practice was to magnetize the capsule at a temperature of 1:-4°. The 
exchange gas was then pumped for an hour and then just prior to demagnetization 
the cryostat was pumped down as far as possible (to about 0-9°). After demagnet- 
ization the heat inflow was about 70 erg/min., most of which was probably 
arriving down the nylon suspension, as this had to be comparatively strong in order 
to hold the heavy capsule steady. ‘The heat flow remained small to temperatures 
of about 1-1° but then gas was evolved and it was necessary to pump the vacuum 
space continually ; even so the heat leak rose to a few thousand erg/min. when the 
capsule was heated up to 1:5°. This, however, was of no great consequence as the 
thermal capacity of the capsule was also rising rapidly. 

Another possible source of heat inflow to the capsule is the eddy currents 
induced in the metal case each time the measuring field (order of 10 gauss) is 
reversed in taking a temperature measurement. ‘This contribution can quite 
easily amount to several hundred ergs at each temperature determination unless 
precautions are taken against it. To make a temperature measurement the 
measuring current is reversed, by disconnecting the measuring coil from the 
battery and then connecting it up the other way round. As the eddy current 
heating for a given change of flux is inversely proportional to the time constant of 
the primary circuit, nearly all of it is occasioned by breaking the coil circuit prior to 
making it in the opposite sense. ‘Therefore we used a reversing switch connected 
up in such a manner that the coil was always shorted before being disconnected 
from the battery. This ensured that the time constant was not unduly small and 
we were able to reduce the eddy current heating to at most 20 ergs per temperature 
measurement. 


§8. THE MEASUREMENTS 


After the calorimeter had been cooled down, the heat inflow to it was measured 
by watching the temperature drift for about ten minutes. Electrical energy was 
then supplied to the capsule at a rate known from measurements of the voltage 
across and the current inthe heater. ‘The temperature changes were observed by 
making ballistic determinations of the magnetic moment of the capsule at half- 
minute intervals for about ten minutes ; the whole cycle was then repeated. A typical 
set of galvanometer readings is given in Figure 3; by differentiating this curve we 

*In order to magnetize the salt isothermally, helium exchange gas is admitted to the space 


surrounding the specimen so that the heat of magnetization may be carried away to the cryostat. 
In order to carry out the demagnetization adiabatically it is of course necessary first to pump away 


all this exchange gas. 
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were able to obtain values of the thermal capacity of the cryostat as a function of 
temperature. Values of this measured capacity (taken from a smoothed curve) are 
given in the Table. 

As has been previously explained (e.g. Kurti and Simon 1938 b) the temper- 
ature scale 7* on which measurements are made below 1° is not the absolute scale 
T but one defined by assuming Curie’s law to hold for the salt even at the lowest 
temperatures +. To convert 7* temperatures to absolute ones it is necessary to 
know the susceptibility of the salt as a function of temperature. We use the 
results of Kurti, Simonand Squire (Cooke 1949) giving T,* and dT,*/dT as functions 
of T to correct temperatures and specific heats to the absolute scale (see Table). 
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Figure 3. ; 
re Thermal Capacities (cals x 10-° per deg.) C108 
a re dljdT, Capsule Case Salt Evap. Helium  (cal/gm/deg.) 

0-4 0-415 0-970 215 12 203 (0) 0 0 
0-5) 0-512 0-979 166 15 142 0 9 5 
0-6 0-610 0-984 138 18 98 0 22 11 
0-7 0-708 0-988 160 21 73 12 54 Li) 
0-8 0-807 0-991 240 De 57 44 118 59 
0-9 0:906 0-994 450 26 46 116 260 130 
1-0 1-006 0-997 810 29 38 250 490 240 
iVei| iol) 1:00 1470 32 31 430, 980 490 
1:2 1-20 1:00 2400 35 26 820 1520 760 
1-3 1:30 1-00 3800 38 22 1310 2400 1190 . 
1:4 1:40 1-00 5500 41 19 1950 3500 1750 
ey aes) 1-00 7800 44 17 2800 4900. 2400 
1-6 1-60 1:00 10300 47 14 3600 6600 3300 


The measured thermal capacity (column 4 of the Table) is made up of four 
components: the capacity of the liquid helium, of the iron alum, of the 
metal capsule and a capacity associated with the latent heat of liquid helium 
evaporating into the free space above the salt and liquid. The specific heat of the 
salt is known, the small capacity of the metal case may be estimated from the known 
electronic specific heats of its components while the capacity associated with helium 
evaporating is calculated in the Appendix. Values of these quantities for a 
capsule filled at room temperatures with 120 atmospheres of helium are given in the 


+ Because of demagnetization effects the measured susceptibility depends on the shape of the 
specimen. Following usual practice our T* temperatures are corrected to T,* temperatures, the 


7T* temperatures that would be obtained if we were using a spherical sample of salt (Kurti and 
Simon 1938 b). 
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Table. The penultimate column in this Table gives the capacity of the helium 
obtained as the difference between the measured capacity of the capsule and the 
total capacity of all the other constituents. 

Clearly below 0-6° the thermal capacity of the helium is so small a part of the 
total that the method becomes insensitive. As a check on the values below 1° 
experiments were done with the same capsule containing the same amount of salt 
but only 15 atmospheres pressure of helium. The difference in the thermal 
capacities of the capsule filled with different quantities of helium gave another 
set of values for the specific heat which agreed with the earlier results to within 
the experimental error. 

The final values for the specific heat C are shown in the last column of the 
Table, in Figures 4(a) and (0) and as a logarithmic plot in Figure 4(c). The 
accuracy is estimated to be within +5°% or +3 x 10-4 cal/gm/deg., whichever is 
the larger. 
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§9.- THE SPECIFIC HEAT OF LIQUID HELIUM 


Between 1-4° and 0-6° we can represent the results by C =0-0247°? cal/gm/deg. 
which agrees closely with the values obtained by Keesom and Westmijze. 

It may be of consequence to note that our measurements have been made 
not on bulk liquid but on liquid in the interstices of a powdered crystalline salt, 
and there might be some dependence of the specific heat on the degree of 
dispersion of the helium.* However, it has recently been shown (Frederikse 1949) 
that while the specific heat of thin films of helium can be much different from 
that of the bulk liquid the difference will be very small if the film is at least a 
hundred molecular layers thick (say 10-°cm.). As the mean diameter of our 
crystals was at least 10-3cm., it is apparent that these size effects would be of no 
importance. 

It is of interest to discuss whether we may extrapolate the T° relation down to 
absolute zero. Both Landau (1941) and Tisza (1947) in their theories of liquid 
helium postulate that at sufficiently low temperatures the only contribution to 
the specific heat will be due to Debye type waves which give rise to the usual 7° 
law. An estimate of the magnitude of such a term has been made (Landau 1941, 

* This was one of the reasons why Pickard and Simon refrained from publishing their results 


but it now seems that the need to determine the absolute value of the energy supplied by their 
gamma-ray source of heat added an additional uncertainty to an already complicated experiment. 
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Dingle 1949) by using Debye’s theory of the specific heat of a continuum but 
omitting transverse waves as HelII has a vanishing viscosity. Using the 
measured value of the velocity of sound they find that the specific heat is given by 
C=4-4x10-4T%cal/gm/deg. which implies 0-0010cal/gm/deg. at 0-6°. Our 
measured value of 0:0011 cal/gm/deg. at 0-6° does not rule out such a term but if 
one is present there must be an abrupt change in the index of the power law at 
about this temperature. 

According to Landau the specific heat follows a T? law from absolute zero 
to about 0-6° when it increases more rapidly because, to use Landau’s terminology, 
the rotons make a contribution which, though only 1°% of the Debye specific heat 
at 0-5°, increases so rapidly that it is 400°% of the Debye term at 1-0°. Values of 
the specific heat calculated from Landau (1947) are given in Figure 5, there is a 
satisfactory agreement between the theoretical and observed values. 
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Clearly it is desirable to obtain data at temperatures below 0-6° to see if the T? 
region appears but this cannot be done using the present apparatus. ‘The 
difficulty is that while the thermal capacity of the salt is rising as T~*, that of the 
helium is falling as 7° and will soon become a very small part of the total capacity 
of the capsule. In fact in order to make measurements down to 0-3° we should 
have to increase the sensitivity by a factor of about 100. We could obtain some 
of this increase by using less salt and demagnetizing from a lower temperature 
but this would increase the difficulty of obtaining good thermal insulation at 
temperatures above 0:5°. Alternatively we might use salts with lower specific 
heats but these have smaller entropies of magnetization and so a greater quantity 
is needed to produce the same cooling. The best method, however, would be to 
use a much greater quantity of helium in the calorimeter. This does not seem 
possible in a capsule of practical size but Hudson, Hunt and Kurti (1949) have 
recently described a calorimeter into which helium may be introduced as a liquid 
and in which the heat inflow due to film has been reduced to a low value. It is 
hoped that by using such a calorimeter measurements will shortly be extended 
to appreciably lower temperatures. 

Although our measurements were made primarily to investigate the region 
below 1:0° it is interesting to compare our values with those of other workers 
above 1-0°. Between 1-0° and 1-4° both our results and those of Keesom and 
Westmijze are in agreement with three points due to Keesom and Keesom (1932). 
Between 1:4° and 2-0° however, the position is not so satisfactory as, apart from 
some early work by Keesom and Clusius (1932) which show much scatter, only 
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three points due to Keesom and Keesom (1932, 1935) are available. Although 
these three points do not seem self-consistent they give values appreciably higher 
than ours at 1:5° and 1-6°; we hope eventually to make measurements in this 
region to clarify the poouee 

As the entropy diagram for helium has already been etiiten (Keesom, 
Bijl and Monte 1941) using Keesom and Westmijze’s value, no modification to 
this is called for. Our results necessitate, however, a small correction to the 
vapour pressure curve of Bleaney and Simon (1939) which was based on the 
results of Pickard and Simon below 1:2°. As this correction will be of the same 
order as the error introduced by the uncertainty in the value of the specific heat 
in the region of 1:5° we do not propose to treat this until more data have been 
obtained in that region. 


IPP END AQ 
The Thermal Capacity Associated with Helium Evaporating in the Capsule 


As the helium in the capsule warms up, some of the liquid evaporates into 
the space above the liquid (about 8cm) and a latent heat must be supplied 
to it. Although the quantity of helium which evaporates between 0-6° and 1-6° 
is a very small fraction of the total, it absorbs an appreciable part of the heat 
input because at these temperatures the latent heat is so much greater than the 
thermal energy. The vapour will also have a thermal capacity which is small 
and may be taken into account by treating it as a perfect gas. We follow Pickard 
and Simon in obtaining an expression for this important contribution to the 
thermal capacity of the capsule. 

If there are n gramme molecules of vapour in the space above the liquid, then 
the rate Y at which heat must be supplied to evaporate the liquid and heat the 
vapour is given by 


V 
Y= Mal 1 >nR erg/deg. aes pl +3e 


dT 2 dT \ RT a ie 
Biel Abe 2 | 73 Ve 
RTA ee VFS ie 


where A= molar heat of evaporation, p= vapour pressure, V = volume of vapour. 
As we may write A=RT?[O(Inp)/0p], we can evaluate AY from the vapour 
pressure curve of liquid helium (Bleaney and Simon 1939). Values of AY as a 
functicn of temperature are given in column 7 of the Table. 
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The Total Cross Section of Beryllium, Aluminium, Sulphur 
and Lead for Neutrons of Energies from 2 Mev. to 6 Mev. 
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ABSTRACT. The total neutron cross sections of beryllium, aluminium, sulphur and 
lead have been measured in the range of energies from 2 Mev. to + Mev. and from 5 Mev. to 
6 Mev. by a transmission experiment. Resonances were observed in all four elements. 
The results obtained for lead and for beryllium are compared with those expected according 
to the one-body theory for light nuclei. 


$1. INTRODUCTION 

HE technique of producing mono-energetic neutron beams has been 

| developed considerably in recent years. Velocity selectors and filtered 

neutron beams are available for the range of energies below 1 Kev.; the 

energy range from 10 key. to 1 Mey. is ahi by a number of (y, n) reactions, 

and in the energy range from 30kev. to 6Mev. there are such reactions as the 

2C (d,n), the -H(d,n), and the *Li(p,n). The parallel development m the 

technique of neutron detection has made accurate cross section measurements 
possible. 

A considerable number of total neutron cross section measurements have been 
reported, particularly in the energy range below 1 kev. (Goldsmith, Ibser and 
Feld 1947). The energy range above about 1 Mev. has not been so well 
investigated. 

Aoki (1939) studied the total neutron cross section of a large number of 
elements in the energy range from 2 Mey. to 3 Mev. by the use of neutrons from 
the ?H(d, n) reaction. Many elements showed an irregular variation in cross 


* Now with the South African Council for Scientific and Industrial Research at the Atomic Energy 
Research Establishment, Harwell. 
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reported by Freee, Fi. Langit and Wiliams (1950) im the range of arson 
es frome OM wer. tw 1-0 wet. and from 25 ev to 3ecr. with carbon, 


The total mentron cross section measurements were made by means of 2 
reamsmasion experiment similar to that used by other experimenters (sce, for 
example, Allen, Barcham and Wilkinson 1947). 

The detector wos 2 pulse iomzzction Gamba alked w 2 pressme of 
a It wes pleced at a distance of 


Fee 1 Eepereecsed] eccomeemen 


i 300cm. from the neutron source and subtended an angle of 
6 minutes at the target. Scattering from the walls, cling and floor was reduced 
by surrounding the detector on all sides with 120cm. of paraffin wax A 
collimating hole 7cm. im diameter passed through the centre. A study of the 
neutrons from 2 monokinetic source had shown that this method of screenmg 
and collimation introduced mo spurious ciiccts. A scale diagram & given m 
Figure 1. 
(a) Neatron Sources 
| Monockinctic neutrons were obtamed from the followmg reacuons- 
G@) 7H(d,2)7He, O value 3-25 +042 wev. (Livesey and Wilkinson 1948); 
(a) **N(d, n) “O, OQ valuc 5-15 + 0-10-mev- (Gibsoa and Livesey 1948). Thick 
targets mounted at the ead of a brass tube 12cm. in length were bombarded 
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with a 930-kev. magnetically resolved deuteron beam. With each reaction the 
neutron energy was varied over a limited range by changing the angle between 
the detector and the direction of the incident deuteron beam. ’ 


(ui) Scatterers 


The lead and the sulphur were cast into blocks, considerable care being 
taken to avoid air gaps inside the casting. The beryllium was in the form of an 
alloy (75:-4% beryllium, 24:6% aluminium). Separate determinations of the 
cross section for aluminium were made at each neutron energy to correct in 
the case of the alloy for its contribution to the measured cross section. The 
aluminium was rolled sheet of 99-5°% purity. 

The scatterers were made several centimetres larger than the diameter of the 
detector and the collimating hole in the wax. They were mounted on a wooden 
stand half-way between the target and the detector. The stand was so placed 
that neutrons scattered off it could not enter the chamber without passing through 
the wax screening. 


§3. PROCEDURE 

In order to check that the attenuation of the neutron beam in passing through 
a thickness x of scatterer could be represented by =, exp (—nox), where I is 
the intensity of the neutron beam of initial intensity J), the number of scattering 
nuclei per cm’, and o the total cross section, measurements were first made of the 
attenuation for several thicknesses of material. Once the above relation had 
been established the cross section was determined by a careful measurement of 
the attenuation for a single thickness of scatterer. 

By the choice of a suitable discriminator bias voltage it was possible to 
eliminate the contribution to the measured cross section from low energy neutrons 
produced by carbon contamination of the targets. At the same time the 
contribution from scattered neutrons not stopped by the wax was also reduced. 

The residual background was determined by placing 100 cm. of wax between 
the target and the detector. With the most unfavourable geometrical arrangement 
used it was found to be less than 0:5°%. 


§4. RESULTS 

For each neutron energy the cross section measurement reported here 
represents the mean of fifteen independent determinations involving a total of 
more than 10° counts. In the range of energies from 2 Mev. to 4 Mev. the overall 
probable error is estimated at +2°; in the range from 5 Mev. to 6Mev. the 
probable error is +3%. 

Owing to the finite size of both the scatterer and the detector some neutrons 
were scattered into the detector; under the geometrical conditions prevailing 
in this experiment this correction was negligible. 

A correction must, however, be applied for the spread in neutron energy 
introduced by the use of thick targets. This correction factor was given by 


Allen et al. (1947) as =i ee ee 
o(E)=o(E)+(E-E)do(E)idE, = sees (1) 


where £ is the maximum neutron energy, E is the mean neutron energy, 


o(E) is the cross section for monokinetic neutrons of energy E, and o(E£) is the 
observed cross section. . 
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oe) 


{ 

This equation assumes that. the neutron spectrum has a sharp cut-off at 
energy E£ and that the absorption of the inhomogeneous neutrons is exponential 
so that a mean cross section can be defined. 

The results are illustrated in Figures 2, 3, 4 and 5 after the above correction 


has been made to the experimentally determined cross section. 


Maxima and 


minima are observed in the results for each of the four elements used. 
The energy resolution varies considerably with the neutron energy. At 
2:2 Mev. the width at half-maximum is 0-15 Mev. and at 3-9 mev. it is 0-5 Mev. 


(Livesey and Wilkinson 1948). 


Between 5Mev. and 6Mev. 


the width at 


half-maximum is 0-4 Mev. (Gibson and Livesey 1948). 
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§5. DISCUSSION 
(i) Lead 


Lead is an interesting element in that the most abundant isotope of mass 
number 208 is made up to 82 protons and 126 neutrons, two of the so-called 
“magic numbers’ which may correspond to the completion of closed shells in 
the nucleus (Mayer 1948, Feenberg and Hammack 1949, Nordheim 1949, and 
others). A considerable amount of evidence has, in fact, accumulated which 
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reveals that some of the lead isotopes behave as light elements with regard to the 
level spacing (see, for example, Kinsey, Bartholomew and Walker 1950, Gittings, 
Barschall and Everhart, 1949, Feld 1949). 

If it is assumed that the observed change in cross section is not due to a 
number of unresolved narrow levels but to a single broad level and that the 
resonance occurs in only one isotope (the most abundant) it becomes possible to 
attempt a comparison of the experimental cross section with that expected 
according to the one-body theory for light nuclei (Breit and Wigner 1936, Bethe 
and Placzek 1937, Macphail 1940, Feshbach, Peaslea and Weisskoff 1947). 
According to this theory the total cross section can be represented by 


ee A a ns A 
aaa ony TES | EER bry, (G2? Pe Ola (2) 


when an incident particle of spin s is scattered by a nucleus of spin #. J is the 
total angular momentum quantum number associated with the compound state. 

Provided the compound nucleus possesses a semi-stable state the width of 
which is narrow in comparison with both the neutron energy and the distance 
between levels the phase shift of the neutron wave, corresponding to given 
values of / the orbital quantum number and of J, with respect to the unperturbed 
wave is given by $.;—S-are wot (R= BRP. le oho (3) 
5 is the contribution of all other levels of the compound nucleus and is 
approximately constant in the neighbourhood of an isolated resonance; 
E,, is the resonance energy; I’, is the width of the resonance level. 

It is also necessary to assume that the inelastic scattering is small in the region 
of the resonance. For lead it is 1-5 barns (Feld 1949) to be compared with a total 
cross section at maximum of 7-7 barns. If these conditions are satisfied then on 
passing through the resonance the variation in \2/z and in 8 is small and sin? by 
will vary between a maximum value of unity and a migimum of zero and the 
cross section will correspondingly change by an amount 

w2 25+1 


max — min = 7 (as+ laity <8 See (4) 


For lead 208 (52-3°, abundance and spin equal to zero) equation (4) leads to the 
computed values of 0-42 barn or 0-84 barn if 7=1 and 0-84 barn or 1-26 barns if 
I=2. The experimental change in cross-section is 1-1 barn. 

At lower energies Barschall et al. (1949, 1950) have reported three resonances 
which are attributed to the isotope of mass number 208. . The resonances 
occurred at neutron energies of 350kev., 525 kev. and 720kev. respectively, 
and /-values of 0, 1 and 1 were assigned to the three states. 

From equations (2) and (3) it follows that the resonance energy E,, and the 
level width I’, are given by 


Ey =1(E, AR) AE Seo «7 ee 6 
l., =(E.—E,) sin 28, sees (6) 


where £, and £, are the neutron energies for which the cross section is a maximum 
anda minimum. The value of 8 is arbitrarily chosen to give the best “shape fit’ 
to the experimental points. 

Withs equal to — 15° it is found that 3-15 < E,, <3-2 Mev., 0-8 <I, <1-15 Mev. 
Most of the uncertainty is in the position of the mmimum. 
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4 If the resonance is cssumed to be due to D-neutrons, then the reduced 
vidth 37 is given im terms of the natural width T (Feshbach, Peaskea and 
Weisskopl 1947) by 3 
sta [SMe te 
; z& (ay ‘ 
Where k=2z/A and a is the nuclear rads (92x10 “em: Sherr 1945). 
‘The reduced width is found to be 13 «10% <7 <19x10*uevcm This 
is im reasonable agreement with the allowed maximum value of 347/2ma 
58 x 10-2 ev. cm. to be expected if the neutron scattering is treated as a single 

icle imtcraction (Wigner 1949). 

The Kfetime of the excited sate given by ~=8T & coloubied © be 
58x10 <> <82«10 sec. The time required for the neutron to cross the 
deus 1s 7-6 x 10 sec. 

- The results presented above on the basis of the smgle particle model must 
xesonance is duc cither to another isotope or to more than one isotope cannot 
be excluded 
| Ie the carer work of Aoki 2 peak m the total cross section at 24 ev. was 
seported. This docs not agree with the resonance reported here (3-2 wev_) and 
the magnitude of the cross section observed by Aoki is somewhat less (see Figure 2). 


(a) Spier 

f Recent mezsurements of the total neutron cross section m sulphur have been 
| made by Adar, Bockeiman and Peterson (1949) m the energy range from | 
W6kev. to 250kev. and by Freier, Fulk, Lampi and Willams (1950) from 
 Qker. to 19 mcr. Both groups used neutrons produced by the *1i(p, 2) 
©) xeaction. A large number of resonance levels were observed, the density of 
which appeared to be so great that the resonances were not completely resolved 
“with 2 Sikev. Li target. 

In the results reported bere im the energy range from 2 uicv. to 6 wev_there 
#s also a rapid variation m cross section with ewergy. For the 285 Mev. 
“resonance”, the observed variation m cross section 1s 0-62 barn and the expected 
wariation for an S-neutron interacting with a target nucleus of spm zero would be 
| 0-97 bam so that here too the resonances are not fully resolved. 


(ii) Alaminiom 

Dhe tena cross section for alumina has been investigated by a number of 
| observers. The energy range up to 1 mev. has been studied im detail and a large 
f umber of resonances have been reported with an average level spacing of 30 kev. 
7 (sce Goldsmath, Ibser and Feld 1947, p. 271, for 2 lst of references). 

Aoks (1939), Zann ef al. (1939), Macphail (1940) and Allen ef al. (1947) have 
made measurements in the range from 2 Mev. to 4 Mev. The published curves of 
; Allen cf al. showed two resonances, one at 2-55 ev. and the second ai 3-) Mev. 
"The measurements reported here do not provide any evidence for the 3-0 Mev. 
i seagnance and while this agrees with the observations of Macphail and of Aoki, 

| the possibilty of 2 resonance 2t 3-)uev. cammot be ruled out because the 
experimental technique is not well suited to the observation of such sharp 
resonances. 
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All results on the 2:55 Mev. resonance agree very well. In the energy range 
from 5 Mev. to 6 Mev. no further resonances were detected. 

As pointed out by Macphail (1940), the ground state of ?”7Al has a spin of 5/2 
and the scattering cross section for a given state is now determined not only by the 
corresponding values of the orbital quantum number but also by the transition 
probabilities between the various components of that state, and, as these are not 
known, an interpretation of the experimental results is not possible. 


(iv) Beryllium 

The results of the present experiment are illustrated in Figure 5 together 
with the relevant data obtained by other observers. The experimental value 
for the change in cross section from maximum to minimum is uncertain but 
probably lies between 1-8 barns and 2-1 barns. A value of 15 is required to 
explain the magnitude of this change in cross section on the basis of the single 
particle model and assuming the spin of the °Be nucleus to be 3/2 units (Kusch, 
Millman and Rabi 1939, Paul 1941). The large value obtained for / may possibly 
be explained (Flowers, private communication) by the fact that the odd neutron 
in the °Be nucleus and the incoming neutron may pair up and appear to add 
angular momenta. A result of this would be that waves of higher angular momenta 
than expected from the incident neutron energy become effective. Moreover, 
as the single particle model requires strong pairing between the incident neutron 
and the bombarded nucleus the results obtained with this model should be 
accepted with reserve in the case of beryllium. 

The following values are obtained for the resonance energy, level width and 
lifetime of the excited state: 2°70<#,,<2:75Mev., 0-5 <I, <0-7 mev., 
13-lx 10-* sec: 7 29-4 x 10-**sec., ‘The time’ required “Wor the: neutron mite 
cross the nucleus is equal to 3-6 x 10~* sec. 
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New Radioactive Isotopes produced by Nuclear 
Photo-Disintegration 


By F. D.S. BUTEMENT 


Atomic Energy Research Establishment, Harwell, Berks. 
MS. received 1st August 1950, and in amended form 22nd November 1950 


ABSTRACT. Photonuclear (y, n) and (y, p) reactions produced by 23 Mev. x-rays from 
a synchrotron have been used to prepare new radioactive isotopes, the method being 
particularly advantageous for neutron-excess isotopes in the region Z=63 to 79. Whenever 
possible the isotopes were chemically identified. The following Table summarizes the 


results. 


Element dlemient 5 , 
: Chemically Half-Life Assignment 
Irradiated ; 
Identified 

Se As 9 minutes 79AS 

Sm Sm 8 minutes 1 SiO beac ot) 

Eu Eu 15 hours Abo 

Gd — 20 minutes 159Ru or Gd™ 

Dy — 14 minutes Ia Movoye 474 Me) 

Er — 96 minutes 168hr, 187Fio or *°°Ho 

Yb — 19 minutes HE ans WH ant ore 9 raat 

Hf Lu 22 minutes NES Eieh (one 202) jul 

W Ta 48 minutes B29 0h 

W Ta 116 hours HOI BE 

Os Re 17 minutes 189Re or 191Re 

Pt Ir 66 minutes TOT Ones lit: 

Hg Au 27 minutes 201Au or 2°8Au 


Definite mass assignments, hitherto uncertain, have been made for 18-hour ©°Gd_ and. 
48-minute 7°°Au. 


S12 INDRODUCTION 


(y, n) and (y, p) reactions. Inspection of an isotope chart shows that, 
especially in the region Z=63 to 79, many new neutron-excess isotopes of 
elements of odd atomic number Z could be made by (y, p) reactions on the 


appropriate Z +1 targets. Some of these isotopes could not be made by any other 
20-2 


‘|e paper describes a search for new radioactive isotopes produced by 
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practicable means, others might be made by (n, p) reactions, but even here the 
photo-disintegration method frequently has the advantage of ease of interpretation 
when a chemical separation is not practicable, since strong activities produced by 
neutron capture are precluded. Certain neutron-deficient isotopes can be most 
conveniently made by (y, n) reactions. Thirteen new isotopes have been found 
and characterized by their half-lives, and where practicable a chemical 
identification. 
§2. EXPERIMENTAL TECHNIQUE 


Samples for irradiation were contained in a cylindrical tube 1 cm. in diameter, 
placed along the axis of the x-ray beam from a synchrotron and 30cm. from the 
target. The x-ray energy spectrum was a continuum up to a limit at 23 Mev. The 
irradiation time was usually one hour, a limit imposed by the heating of the 
synchrotron magnet. The samples received approximately 2000r. per hour. 
Photo-neutrons, both fast and of thermal energy were present, the flux of each 
at the sample being about 1000sec"tcm-?. The samples were therefore screened 
with 1mm. of cadmium foil to prevent activation by slow neutrons, and no 
half-life attributable to a product of neutron capture was ever observed, except 
for 1®Dy. 

For (y, p) products where chemical separation allows the activity to be 
concentrated into a small mass, 1 to 10gm. of target material was irradiated and 
activities of the shorter-lived products of up to about 1,000 counts per minute 
obtained. 

For (y, p) products of certain rare earths, where a rapid chemical method of 
separation is not available, and for (y, n) products, quantities of 0-1 to 0-2 gm. were 
irradiated. Activities of about 100 counts/min. were obtained. With an 
end-window Geiger counter there is little advantage in using larger quantities 
owing to self-absorption of the radiation by the sample. 

The decay of the activity was followed by means of a thin-window 
(1:7mg.cm~*) Geiger counter. ‘The activities found were nearly all of short 
half-life, since it was seldom possible to produce sufficient activity of an isotope 
of half-life greater than one day. 

In order not to miss short activities, chemical separations were made as 
quickly as possible, using methods based on the assumption that only the products 
of (y, n) and (y, p) reactions were present. Activities with half-lives considerably 
less than the time necessary for the separation would not have been observed. 
When oxides were irradiated and counted without any chemical separation, 
half-lives less than 2 minutes would not have been observed owing to the strong 
126-second ¥°O activity. ‘The activities available were usually too weak to permit 
any determination of the energy of the particles emitted. 


§3. REACTION YIELDS 


When an irradiation produces more than one isotope, the ratio of the yields 
provides a useful indication of the type of reaction occurring. This ratio (i.e. the 
ratio of the number of atoms of each radio-isotope produced) can be calculated 
from the half-lives of the isotopes, their activities measured at a known time 
after the end of the irradiation, and the duration of the irradiation. 

The results will be fairly accurate for the thin samples of beta-emitters 
prepared by a chemical separation, since the unknown but small correction 
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factors for back-scattering and window absorption will be about equal for the 
several activities in the sample. For thick samples, i.e. where no chemical 
separation is possible, there will be a larger error owing to the impossibility of 
correcting for the self-absorption for the new isotopes of unmeasured beta-particle 
or positron energy. 

It has been shown (Price and Kerst 1950) that the value o of the absolute 
cross section for a (y, n) reaction is nearly the same for all isotopes of an element, 
and increases steadily with increase of Z. Also Hirzel and Waffler (1947) found 
that the ratio o (y, p)/o (y, n) has values about 1/30 at Z=49-50 and decreases 
as Z increases. ‘The approximate value obtained in this work was 1/200 at Z=78. 

Consequently, when two radio-isotopes of the same element are produced 
simultaneously (by either two (y, n) or two (y, p) processes) and the mass number 
of one of them is already known, then, if the mass number assigned to the other, 
new, radio-isotope is to be plausible, the ratio of the yields should approximately 
equal the ratio of the abundances of those stable isotopes from which the 
radio-isotopes were produced. 

The results of irradiations are described below. Data quoted without 
reference are from the Table of Isotopes by Seaborg and Perlman (1948). 


§4. ARSENIC 
The object was to produce the unknown “As by ®Se (y, p). The stable 
selenium isotopes are as follows: 


Mass number 74 76 Hil 78 80 82 
Abundance (%) 0-87 9-02 7:58 23°52 49-82 9-19 


The radioactive arsenic isotopes which could be produced by (y, p) are : 


Mass number 76 76 77 oS 81 
Half-life ~ 90 days 26:8 hr. 40 hr. ae <10 min. 

The activity due to “As would be negligible. The half-life of *1As is given 
by the Plutonium Project (1946) list of fission nuclei. 

Ten grammes of potassium selenate was irradiated and the radio-arsenic 
chemically separated. The time taken for the isolation was 30 minutes. The 
activity showed a half-life of 9 minutes and a weak residual activity with an 
apparent half-life of about 31 hours. The latter may be attributed to a mixture 
of 26-8-hour “*As and 40-hour 77As. The yields of the 9-minute and ‘31-hour’ 
activities were approximately equal. The 9-minute arsenic is therefore 
probably 7°As, decaying by beta-particle emission into 7Se whose half-life is 
either very short or very long. It cannot be *4As since it does not produce any 
17-minute or 57-minute *!Se daughter activities. 


§5. SAMARIUM 


The natural samarium isotopes are as follows: 


Mass number 144 147 148 149 150 152 154 
Abundance (%) 3-16 15-07 11-27 13-84 7-47 26:63 IM Blass 


The samarium isotopes 143, 146, 151 (20-year) and 153 (46-hour) would be 
produced by (y, n) reactions, and (y, p) reactions would produce much smaller 
yields of seven promethium isotopes. 
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A sample of 200 mg. irradiated samarium oxide of over 99-9% purity showed, 
after the decay of the !°O, two activities with half-lives of 8 minutes (Figure 1) 
and 46 hours. The yield (not corrected for self-absorption in the sample) of 
the 8-minute activity was about one-quarter that of the 46-hour activity. 

These activities were both identified as samarium by irradiating samples of 
samarium dichloride and chemically purifying the samarium. The 46-hour 
activity must be !%Sm. The long-lived %41Sm would not be observed. Apart 


e Total Activity 


s Activity | corrected 
°Sm 


for” 


Activity (counts per minute) 


Half-life 
8 min. 


Minutes 


Figure 1. Decay of radioactive samarium. 


from possible isomerism the 8-minute activity must be 148Sm or ®Sm._ Since 
one natural samarium isotope is an alpha-particle emitter, absorption 
measurements in aluminium were made of the 8-minute activity, which showed 
that the radiation detected was either positrons or electrons. With the counter 
employed, a samarium isotope decaying only by orbital electron capture would 
scarcely be detected. In view of the relative abundances of 44Sm and 1°4Sm, 
the yield of the 8-minute activity relative to Sm is consistent with its assignment 
to a positron-emitting #8Sm, decaying into 1#%Pm, which would not be observed 
since its half-life is 350 days (Hicks 1949), 


§6. EUVROPIUM 

Europium has two stable isotopes of mass numbers 151 (47:8°%) and 
153 (52:2%). A (y,n) reaction would produce Eu (9-2-hour, and a 5-year 
isomer) and the hitherto unidentified ©°Eu. Pool and Quill (1938) produced 
a 27-hour positron activity by bombarding europium with fast neutrons and 
suggested an assignment to °Eu, but a chemical identification was not made. 
A sample of 100 mg. of ‘Specpure’ europium oxide was ire uiaced and the 

europium chemically purified as europous sulphate. 

The activity showed a nearly logarithmic decay with an apparent half-life 
of 12 hours. Such a decay, resulting from two slightly different half-lives, is 
difficult to resolve accurately unless one component is known. By assuming one 
component to be 9-2-hour Eu which must have been formed, the best value 
for the half-life of the other component was 15 hours, and its yield was nearly 
equal to that of the ©?Eu. The 15-hour activity must almost certainly be ®°Eu, 
and since the abundances of 1°1Eu and Eu are nearly equal, the relative yield 
indicates a high counting efficiency, so that decay is probably by positron emission. 
With the counter employed, a Eu decaying only by orbital electron capture 
would scarcely be detected. 
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§7. LUTETIUM 


The object was to prepare 1”8Lu and !”°Lu by (y, p) reactions on hafnium. 
The stable hafnium isotopes are as follows: 

Mass number 174 176 17 178 179 a PLSo 
Abundance (%) 0-18 a6 18-47 Z7ANO 13-84 35-11 

A sample of 0:8 gm. hafnium oxide, free from zirconium, was irradiated and 
the lutetium chemically separated. 

The half-lives found were 22 minutes, 6 days (evidently the known 1”*Lu) 
and an intermediate period whose decay was not quite logarithmic, but appeared 
to be a mixture of periods of 3-5 hours and about 5 hours. It may therefore 
indicate the presence both of 3-5-hour 1”*Lu and one of the new isotopes sought, 
but no certain conclusion can be reached solely on the evidence of the decay. 

The 22-minute lutetium is a new isotope, and consideration of its yield 
(about 70°, of that of the !7*Lu) indicates that its mass number must be 178 or 179. 


§8. TANTALUM 
The object was to prepare 18°T'a and 1®°T'a by (y, p) reactions on tungsten. 
The stable tungsten isotopes are as follows: 
Mass number 180 182 183 184 186 
-Abundance (°%) 0-12 Doty. 14-24 30-68 ZIAT 


In a preliminary experiment tungsten metal was irradiated. The main yield 
of tungsten isotopes, produced by (y,n) reactions, will be the long-lived '8'W 


e Total Activity 
@ Activity corrected for II6-hr.Ta 


Half-life 
48 min. 


Activity (counts per minute) 


I6-hr. Ta 


) 100 200 300 400 500 
Minutes 


Figure 2. Decay of 48-minute tantalum. 


(140 days) and !°W (74 days) which will show only a feeble activity, so that 
short-lived tantalum isotopes will be observable. ‘The shortest half-lives shown 
were 17 minutes (evidently 16-2-minute 1*T'a) and 48 minutes. ‘To isolate the 
tantalum activities 14 gm. of tungstic acid was irradiated and the tantalum 
chemically separated. 

The decay graph (Figure 2), commencing after the disappearance of the 
17-minute activity, shows half-lives of 48 minutes and a long period, originally 
given as 6 days (Butement 1950 a), but subsequently re-determined much more 
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accurately on a larger sample as 116 hours (Figure 3). A sample of this 116-hour 
activity, isolated after decay of the 48-minute activity, showed after 37 days a 
residual long activity of 2-7°% of the initial activity, which is about that to be 
expected for the yield of 117-day 182T'a. The 116-hour activity cannot be *°Ta, 
since there should then have been also an additional residue of 4-9% of the initial 
activity, due to its daughter 74-day 1°W. 


© Total Activity 
= Corrected Activity 


Half-life 
16 hr. 


Activity (counts per minute) 


Long-lived activity 


0 200 400 600 800 
Hours 
Figure 3. Decay of 116-hour tantalum. 


The yields of 48-minute and 116-hour tantalums were in the ratio of 1: 1-2 
respectively, which is in conformity with the nearly equal abundances of 
184W and 86W. Probably the 48-minute activity is Ta and the 116-hour 
activity 18°'Ta. 

§9. RHENIUM 


The object was to prepare 18°Re and }*Re by (y, p) reactions on osmium. 
The stable osmium isotopes are as follows: 


Mass number _—184 186 187 188 189 190 192 
Abundance (%) 0-018 1:59 1:64 13:3 161 264 41-0 


The only appreciable yield of radioactive osmium isotopes produced by (y, n)} 
reactions will be 97-day 1*°Os and 16-day 1!Os, so that short-lived rhenium 
isotopes, produced by (y, p) reactions on osmium, will be observable without 
chemical separation. ‘Specpure’ osmium metal of 99-997, purity was 
irradiated. ‘The two shortest activities observed were 17 minutes and 
approximately 18 hours. There was also a weak longer-lived activity. The 
18-hour activity is presumably Re. Rhenium chemically separated from 
irradiated osmium showed a short period, approximately 15 minutes half-life, 
and a weak longer-lived residual activity, partly 18-hour **Re. Owing to the 
decay of activity during the time taken for chemical separation, the short rhenium 
half-life was more accurately determined as 17 minutes using irradiated osmium, 
without chemical separation. 

The only radioactive rhenium isotopes which will be formed in appreciable 
yield will be 86Re (92 hours), #88Re (18 hours), and the unknown !8°Re and !*!Re. 
The yield of the 17-minute rhenium was slightly greater than that of the }8Re. 
Apart from possible isomerism, therefore, 2s 17-minute rhenium must have a 
mass number of 189 or 191. 
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; §10. IRIDIUM 


The object was to prepare 1°Ir and !9Ir by (y, p) reactions on platinum. 
The stable platinum isotopes are as follows: 


Mass number 192 194 | 195 196 . 198 
Abundance (%) 0-78 625 33°7 25-4 7:23 


Irradiated platinum metal showed activities with half-lives of 80 minutes 
and 18 hours, both attributable to platinum isotopes. It is unlikely that any 
iridium isotope with a half-life of between 1 and 20 minutes is formed, since such 
should have been observable. 

Five grammes of chloroplatinic acid was irradiated, and the iridium isolated 
chemically. The iridium activity decayed with half-lives of 66 minutes and 
19 hours, the latter period being evidently the known Ir. 

The yield of the 66-minute iridium was about one-fifth that of the Ir. 
Apart from possible isomerism, it must be Ir or 1°Ir, and the yield relative to 
Tr suggests 1°7Ir as the more probable assignment, in view of the relative 
abundances of the platinum isotopes 198 and 195. 


§11. GOLD 


The object was to prepare 2 Au and Au by (y, p) reactions on mercury. 
The stable mercury isotopes are as follows: 


Mass number 196 198 199 200 201 202 204 
Abundance (%) 0-15 10-1 17-0 23°3 13-2 29-6 6-7 


The shortest half-life shown by irradiated mercury metal was 40 minutes, 
evidently due to the known isomer (mass unknown) of one of the mercury isotopes. 

It is unlikely that any gold isotope with a half-life between 1 and 5 minutes 
is formed, since such should have been observable. Longer half-lives would 
have been obscured by the mercury activity. 

Twelve grammes of pure mercury was irradiated and the radioactive gold 
chemically separated. The time taken for the isolation was 28 minutes. The 
activity of the gold decayed with half-lives of 27 minutes, 48 minutes, and 
approximately 3 days. 

The 3-day activity is evidently a mixture of the known 2-7-day 1°8Au and 
3-3-day Au. The 48-minute gold has been made by Maurer and Ramm (1942) 
by the reactions Hg (n, p) and TI (n, a), thus proving the mass number to be 
either 200 or 202. Its production also by Hg (y, p) shows that its mass number 
must be 200. Apart from possible isomerism the 27-minute gold must be either 
201 Au or 2%Au. The relative yields of 27-minute : 38-minute : 3-day activities were 
approximately 1:1:1-5 respectively, so that either assignment appears equally 
probable. 


§12. RARE EARTH ELEMENTS 


With most of the rare earth elements there is no practicable method of 
chemically identifying weak short-lived activities, but evidence for new isotopes 
can be obtained by irradiating samples of the pure oxides, and measuring the 
decay of activity after the 126-second O has disappeared. ‘The oxides were 
ignited at 1,000°c. shortly before irradiation to eliminate possible contamination 
by carbonates which could yield 20-minute “C. The results are given below 
under the headings of the element irradiated. 
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§13. GADOLINIUM 


‘Specpure’ gadclinium oxide showed activities with half-lives of 20 minutes 
and 18 hours. 

The latter period is also that of the beta-emitter produced by neutron capture 
in gadolinium, and is therefore °Gd, produced here by the reaction 1°°Gd (y, n). 
This confirms the mass assignment suggested previously (Butement 1949). 

There was no known contamination to account for the 20-minute activity. 
Its yield relative to the ®Gd, not corrected for self-absorption by the sample, 
was 1/80. This is about the yield to be expected for ™*Eu. The activities of all 
the other and longer-lived europium isotopes produced would not be observable. 
Alternatively it may be an isomer of a stable gadolinium isotope. It is not likely 
to be 4Gd produced by !Gd (y, n) since bombardment of europium oxide 
with 10 Mev. protons, which would produce Gd, did not yield any 20-minute 
rare earth activity. 


§14. TERBIUM 


A (y, n) reaction on terbium would produce !8Tb. Irradiated terbium oxide 
showed no significant activity, so that the half-life of Tb is probably long. 


§15. DYSPROSIUM 


A sample of dysprosium oxide containing 0-3°% Ho showed activities with 
half-lives of 22 minutes and 140 minutes. The latter period is evidently 1®Dy 
formed by neutron capture, since the yield was much reduced by increasing the 
thickness of cadmium screening. The short period, given as 22 minutes in earlier 
work (Butement 1950 a) might be inaccurate owing to contamination with '*4Ho, 
the half-life of which has been found to be 34+1 minutes (Hicks 1949). 
Consequently, a sample of very pure dysprositum oxide was prepared by the 
ion-exchange method, and this showed a half-life of 14 minutes. The dysprosium 
isotopes formed by (y, n) reactions will be the long-lived #°*Dy (Butement 1950 b) 
and negligible yields of *Dy and "Dy from the scarce stable isotopes ®&Dy 
and Dy. The 14-minute activity is therefore probably a terbium isotope of 
mass number 162 or 163. 


§16. ERBIUM 


Irradiated ‘Specpure’ erbium oxide showed an activity with a half-life of 
44 minutes (Butement 1950 a). This period is too long to be due entirely to 
34-minute 1®Ho derived from a (‘y, n) reaction on holmium contamination, but 
could be due to a mixture of a 34-minute and a somewhat longer half-life. An 
enalysis of the erbium, made by using the pile-activation method and measuring 
the 27-hour 1°*Ho activity produced, showed a holmium content of about 0:-5%. 

By adding to the erbium a trace of '6Ho and then eliminating it by an 
ion-exchange separation, a sample of erbium oxide containing less than 0-01°% 
of holmium was prepared. After irradiation this showed activities with half-lives 
of 96 minutes, 10 hours and 9 days. 

The 9-day activity is evidently 9:-4-day !°Er produced by a (y, n) reaction. 
A comparable yield of 1®Er will also be produced, with a half-life of 10 hours 
(Butement 1950c), but this shows only a low activity since it decays by orbital 
electron capture, so that the counting efficiency is low. The yield of the 
96-minute activity was about 1/190 that of the Er, after making a correction 
for the self-absorption of the 0-33 Mev. beta-radiation from the 1®°Er in a sample 


New Radioactive Isotopes 403 


of 25mg.cm- thickness. This activity may be 'Er derived from Er whose 
abundance is 1:5%, or may be 16’Ho or 1®°Ho produced by a (y, p) reaction on 
168Er or 1?°Er whose abundances are respectively 26-9°% and 14:2%. 


Soe ele ER BIE 


Irradiated ‘Specpure’ ytterbium oxide showed activities with half-lives of 
19 minutes and about 100 hours, the latter evidently 1=Yb. The only other 
ytterbium isotope which would be formed in appreciable yield is 33-day 1®°Yb, 
which would not be observed. The yield of the 19-minute activity was about 
one-hundredth that of the 1° Yb, after making a correction for the self-absorption 
only of the 0-5 and 0-13 Mev. beta-radiations from the !75Yb in a sample of 
50mg.cm™ thickness. This 19-minute activity is probably a thulium isotope 
of mass number 172, 173, or 175. 


§18. LUTETIUM 


Irradiated lutetium oxide showed only the known 3-5-hour !7*Lu isomer, 
and the natural radioactivity of 1®Lu. ‘The half-life of 1“4Lu which must be 
produced by a(y, n) reaction on 1”°Lu is therefore probably long. 
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On the Proton Component of the Vertical Cosmic-Ray Beam 
at Sea Level } 
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The Physical Laboratories, University of Manchester 
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ABSTRACT. Approximately 1°% of the vertical cosmic-ray intensity at sea level is shown 
to consist of protons. These form a momentum spectrum of index —2-8 in the range 
103 Mev/c< p< 5X10? mev/c. and have an absorption range in lead, 160+40 gm.cm-?. 
Using existing measurements at high altitude, the attenuation in the atmosphere of protons 
in the spectral region investigated is found to take place with a range 140+10 gm.cm~. 

There is no negatively charged component of strong nuclear interaction in the vertical 
beam of intensity as great as 0:1°% of the whole beam. 


§1. INTRODUCTION 

XPERIMENTS have been carried out to examine the absorption properties 

of charged cosmic-ray particles as a function of momentum, particles of 

measured momentum now being available in the output of the magnetic 
spectrograph described in an earlier paper (Hyams, Mylroi, Owen and Wilson 
1950). ‘The work was in the first place expected to give information about the 
range-momentum relation for u-mesons, and to show the presence of any 
measurable fraction of more absorbable particles in the beam, in a manner 
similar to that described by Alikhanian et al. (1948) and by Lheritier, Peyrou and 
Lagarrigue (1947). Exploratory experiments showed at a very early stage that a 
reasonable number of more absorbable particles was in fact present, and the forms 
of apparatus described in the present paper were adopted for the specific purpose 
of studying these particles and identifying them. 


§2. EXPERIMENTAL ARRANGEMENT 

The experiments consist essentially in recording from an upper counter tray 
the entry of particles of measured momentum into an absorbing layer, and also 
from a further tray below the absorber to determine whether or not the particles, 
emerge. Preliminary experiments indicated that the proportion of absorbable 
particles reaches a value which may be as large as 5 °4of the whole flux at momenta 
near 10° mev/c. and falls rapidly at higher momenta, and that the absorbable 
component shows a strong positive excess. The apparatus was therefore 
dimensioned on the assumption that the main part of the absorbable component 
consists of protons. ‘The absorbing thickness was made small enough for the 
critical momentum of penetration for mesons to be significantly less than the 
momentum of the anticipated maximum of a proton spectrum (set by the increase of 
ionization loss from an inverse power spectrum), but nevertheless larger than the 
collision range corresponding to the geometrical cross section of the absorbing 
(lead) nuclei. A 20 cm. absorbing layer of lead was adopted, which u-mesons of 
momentum greater than 0-37 x 10° Mev/c. should be able to penetrate, and which 
is about 1-5 nuclear collision lengths in thickness. The use of this, and even thicker, 
absorbing layers is straightforward on account of the small angular divergence of 
the measured particles from our instrument. Only a relatively small increase of 


- 
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t 
area of the lower counter tray, which detects the emergence of particles from the 
absorber, as compared with the upper tray is necessary to allow both for the 
divergence of the entrant beam and for scattering of particles in the absorber. 

The greater part of the work discussed has been carried out with the two 
arrangements shown in Figure 1 (a) and (db). In both arrangements the entry of 
measured particles to the absorbing layers was detected by tray A and final 
emergence at tray EFGH. ‘The use of intermediate trays to obtain more detailed 
information about the absorption process is discussed below (§ 3.2). 


(a) (b) 


Bottom Counter Row 


(WWYAY of Spectrograph \. (HH) 
ae @@O® 


pe 
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Som, —— 


Pecannol 
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es ee ee 
0 10 20 30 40 cm. 


Figure 1. The absorption apparatus in relation to the bottom 
counter layer of the spectrograph. 


2.1. Efficiency of Detection in the Lower Counter Tray 

Since effects of the order 1°% of the total beam are to be examined, it is necessary 
for tray EFGH to detect particles emerging from the absorber, after traversing it, 
with highefficiency. It accordingly consists of two tightly packed layers of counters 
staggered relative to each other. With this arrangement, about 75% of all 
particles traversing the tray pass through the sensitive region of two counters, 
while the remaining 25°% pass, in one layer or the other, through the dead space 
between counters. (The insensitive space in each layer, ~12%%, arises to some 
extent because the counter tubes (80 cm. long) are not all completely straight.) 
The dead-time of the individual counters, which was measured to be about 
120 psec., leads to an inefficiency per counter of about }%; the total inefficiency 
of the whole tray is thus about $°%. This figure was checked in an experiment 
without absorber. Of 3,700 particles traversing the lowest tray, six failed to 
operate it. 

It is further necessary to ensure that no significant fraction of particles is 
scattered in the absorber so as to miss this tray. ‘To test this possibility, and to 
maintain a continuous check on the working efficiency of the various counters in 
the tray, the sub-groups, E, F, G, H (Figure 1(a)), were recorded separately. 
Thus the most scattered particles are recorded at H alone, and the next most 
scattered at EH. Of 8,000 traversals (p> 6 x 108 ev/c.), six are of the type EH and 
four of type H. We therefore conclude that the cover of the tray EFGH is 
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adequate and that scattered particles which miss this tray can hardly amount to 
more than 0-05 °% of those traversing it. 

The continuous efficient operation of the various counters is checked from the 
actual records by the fraction of all counts which discharge only one of the 
subgroups G, E, F. 

§3. RESULTS 
The results of the main experiments are summarized in Tables 1 and 2. 


Table 1. Absorption of Particles in 20 cm. Lead 


Arrangement as Figure 1 (a). 


Momentum Positive particles Negative particles 
Category Mean (10° Mev/c.) Total Stopped Total Stopped 
0 12 (967) (3) — - 

1 6 1150 8 888 + 
D 3 983 8 799 3 
3 2 739 18 566 2 
4 125 464 5 360 1 
5 iten 363 10 285 2 
6 1-0 236 14 163 2 
d 0-86 180 14 147 0) 
8 0-75 125 7 94 il 
9 0:67 65 7 67 0 
10 0-60 57 5 42 1 
Total 4358 106 3464 16 
+ (967) = (3) 
Table 2. Absorption of Particles in 10, 15 and 20 cm. Lead 
Arrangement as Figure 1 (0). 
Momentum Positive particles Negative particles 
C Mean Total Stopped in (cm. Pb) Total Stopped in (cm. Pb) 
e187)? mey/onie or 1s 0-10) 0=05 9 45200 OF 10-10 ak Oo )5ap1S=20 
0 12 (496) - - (3) — - - ~ 
1 6 576 - 2 1 469 - - 2 
2 3 494 2 il 1 355 - - 3 
3 2 337 4 1 1 307 ~ - — 
4 5 277 8 1 1 194 - - 1 
5 ey 182 3 - 1 118 1 ~ - 
6 1:0 110 3 2 - 77 2 - - 
a 0-86 76 2 2 1 55 = ~ 1 
8 0-75 57 py) - - 41 - - - 
9 0:67 30 3 ~ Di 1 - 1 
10 0:60 30 4 1 - 15 - - - 
Total 2169 31 10 6 1974 4 0 8 
+ (496) +(3) 


3.1. Instrumental Background 


We first observe that the absorption of negative particles is, in both arrange- 
ments, small and independent of momentum: we tentatively regard this apparent 
absorption of negative particles as representing the instrumental background. In 
the first experiment in which the detecting efficiency of the main discriminating 
tray EFGH was very carefully controlled, this background is about 0-45%. It is 


i tocol 
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thus slightly higher than that calculated from the properties of individual counters 
(absorbing system 0-15°%, spectrograph ~0-1°), and therefore indicates slight 
inefficiency of maintenance of the counters in tray EFGH which is not detected in 
the check processes outlined above. In the second experiment the negative 
background is rather higher, but the particular layers in which the particles appear 
to have stopped are very significant. ‘The majority (0-4°%) appear to stop in the 
last 5 cm. of lead: that is to say, the evidence that they have stopped is based only 
on the failure of tray EFGH to discharge. Ifa particle which in fact traverses the 
whole absorber is recorded as stopping in the first (0-10 cm.) or second (10-15 cm.) 
layers of lead, then it must also fail to discharge counter layers UL and BCD in 
addition to layer EFGH. ‘These layers are of relatively low geometrical efficiency, 
about 85°% only, nevertheless a spurious record of stopping in 10-15 cm. lead is to 
be expected only about one-sixth as often as in the lowest layer, while a correspon- 
ding spurious record should occur with negligible frequency for the first 10 cm. lead. 
That 0:4°% of all negative particles appear to stop in the last 5 cm. of the lead 
absorber, while only 0-2°% appear to stop in the first 10 cm., and none at all in the 
intermediate layer, is entirely consistent with the assumption that this 0-4°% 
arises from instrumental inefficiency, but would represent a most remarkable 
property if it were interpreted as real absorption. 

The four negative particles which appear to stop in the first 10 cm., about 
()-2°% of the whole negative beam, present a different problem and will be discussed 
later. 


3.2. Absorption of Positive Particles 


It has been shown that only a very small part of the negative particle beam can 
be considered to be absorbed in 20cm. of lead. In contrast, considerable 
absorption of positive particles appears in both Table 1 and Table 2, and we 
proceed to examine the absorption and the total flux of these particles. To 
determine the total flux of absorbable positive particles, the efficiency of absorption 
in the experimental arrangement must be estimated, and the records of the inter- 
mediate absorbing layers, BCD in Figure 1 (a), UL and BCD in Figure 1 (8), are 
used for this purpose, and at the same time to give the range of the particles in lead. 
The relevant observations are given in Table 3, which refers only to positive 
particles. ‘The records in layer BCD of arrangement (a) have not been given in 
Table 1; the corresponding intermediate values for arrangement (b) appear in 
Table 2 and are transcribed. For the reasons developed in the previous paragraph, 
the most significant comparison in the second experiment is between particles 
stopping in the first 10 cm. and in the next 5 cm. lead. 

In determining absorption from tne data of Table 3 it is necessary to take 
account of the ionization range of the particles (which we assume to be protons). 
The momenta of protons which in absence of nuclear interaction will penetrate 
respectively 5, 10, 15 and 20 cm. lead are 0-66, 0-80, 0-93 and 1-05 x 108 Mev/c. 
(Princeton Range Tables, E. P. Gross 1947). Although each momentum category in 
fact includes particles of an appreciable range of momentum (Hyams et al. 1950), 
it will be sufficiently accurate to divide the measurements with arrangement (a) 
into two groups: categories 1-5, for which protons would not be stopped by 
ionization in the lower absorbing layer, and categories 6-9, for which all protons 
will be stopped in this way if they do not sooner make a nuclear collision. With 
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arrangement (b) only categories 1-6 yield relevant information. The derivation of 
the absorption range is summarized in Table 4; the value deduced is 


(160 + 40) gm.cm™? (14 cm. lead). 


Table 3. Absorption of Positive Particles in Lead 


Momentum Arrangement (a) Arrangement (5) 
Mean Stopped Stopped 
Category (10° Mev/c.) Total O-5cm. 5-20 cm. Total O0-10cm. 10-15 cm. 
0 12 (967) — (3) (496) - = 
1 6 1150 D, 6 576 - 2 
2 3 983 1 7 494 2 1 
3 2 739 6 12 B37 + if 
4 i125) 464 6 9 ZUG; 8 1 
5 ep 363 + 6 182 3 - 
6 1-0 236 3 11 110 3 2 
7 0-86 180 6 8 76 2 2 
8 0:75 125 2 5 57 2 - 
9 0-67 65 4 3 30 3 
10 0-60 Dy, 2 3 30 4 
Total 4358 3i7/ 73 2202 31 10 
Estimated background -+(967) 20+4:°5 —+(496) 


Table 4. Absorption of Protons (of spectral distribution shown in Figure 3) 


in Lead 
(1) (2) (3) (4) 
, 1— exp (—5/A) 
(@) 1-5 exp (—5/A)—exp (— 20/A) 0:54-0:17 20 
1— exp (—5/A) 
(b) 1-6 MecenD sn) 3-041-4 19 


exp (—10/A)—exp (—15/A) 


(1) Arrangement and momentum categories considered. (2) Calculated ratio of particles 
stopped in successive layers as a function of range Acm. (3) Observed value of ratio 
(in lead) frt+e. (4) Deduced value of range (cm. lead). 

a The mean range A and its standard deviation are given, writing €(Aj) = (0A;/Afie(fi), by 

A=LEHA/PA)N/LE1/PA)}], (A) =3/2H1/e?(Aj)}. 


3.3. Proportion of Protons in the Vertical Beam 

The value 160 + 40 gm.cm, referred to in the previous section, although not 
accurate, allows an estimate to be made of the relatively small number of protons 
not stopped in our absorber. We give in Table 5 the fraction of protons present 
in the accepted beam of our apparatus as a function of momentum. Since we are 
concerned here with a momentum range for which incident electrons are rare, and 
since (§ 4.5) we have in this work used a screen above the spectrograph for the 
purpose of degrading the remaining electrons; this fraction may be taken to refer to 
the hard component in directions within 10° of the zenith (Figure 2). 
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Table 5. The Proton Component of the Vertical Beam (within 10° of 
the zenith) as a percentage of the total hard component within the 
‘same angle 


Note: The uncertainties throughout the table are regarded as purely 
statistical, and for convenience are only stated in the final column. 


Mean 
Sere Arrangement (a) Arrangement (5) Total 
(10? Mev/c.) (1) (2) (3) (4) (5) (Din 4y2e) (1) (6) (7) 
42 967 3 45 —15 — 496 -— — 1463 —1:5 —0-1+0-1 
6 2038 8 54 2:6 0:8 1045 2 0-9 3083 6°3 0-2+0:1 
3 1782 8 4:6 3-4 1:0 849 3 1-4 2631 8-8 0-3+0:1 
2 1305.18 3:5 14:5 4:5 637. 55.9 2:3 1942 26:3 1-4+0°3 
1°5 877 15 2-2 12:8 4-0 471 9 4-2 1348 30:0 2°2+0°4 
1-2 645 10 1:9 8-1 2°5 300 3 1:4 945 15:0 2-4-.0-4 
1-0 399 14 1:1 12:9 2-1 AST eSiales 586" 72123 PF 
0:86- 327-14 0-8 13-2 -— (igh Zh cee 458 17-2 3-340°6 
0-75. 219 7.06. 64-— — 98 2- = 317 8-4 
O67 132-7" 073 6-7 — SI) 189 9-7 5-8-4155 
0-60 99 5 0:3 4:7 — 45a 144 9-7 f. 


(1) All particles. (2) Positive particles appearing to stop. (3) Estimated instrumental 
background. (4) Protons stopped. (5) Estimated number of protons not stopped. 
(6) Number of protons. (7) Percentage of protons. 
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Figure 2. Figure 3. The measured spectrum of protons 


in the vertical direction at sea level. 


@ Measurement on slow protons by Rochester 
and Bound (1940). 
© Sea-level measurements (this paper). 
! + Measurements at 30,000 feet after Adams 
et al. (1948) (at arbitrary intensity). 


§4. DISCUSSION OF RESULTS: DEVELOPMENT 


4.1. The Spectrum of Protons at Sea Level 
The spectrum of the protons is derived immediately by application of the results 
of Table 5 to the known meson spectrum; for this purpose the meson spectrum 
given by Rossi (1948) is sufficiently accurate. In Figure 3 we plot on double 
logarithmic scale the sea-level vertical meson spectrum and the proton spectrum 
PROC. PHYS. SOC. LXIV, 4—A 27 
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deduced inthis way. ‘The slope of the proton curve leads toa differential spectrum 
of the form N(p)dp =Kp-*8 dp, K=1-0 x 10-* cm* sec sterad! if p and dp are 
measured in units 10° Mev/c. If we assume this momentum spectrum to extend to 
infinity, the vertical flux of protons with momentum greater than 10% Mev/c. is 
0-55 x 10-4 cm sec“! sterad“1, with about 98°% of the particles falling in the 
measured region 1-10 x 10° Mev/c. The spectrum derived from our measurements 
does not extend to momenta lower than 700 Mev/c., for at lower momenta the 
results are complicated by the instrumental spread of the meson absorption range. 
We have therefore no direct evidence of the reduction of proton intensity which 
must take place when the ionization loss by protons in air becomes large. Rochester 
and Bound (1940) found a total flux of about 0-5 x 10-® cm”? sec“ sterad-t 
of residual range at sea level between 10 and 35 gm.cm™ lead. This 
spread of range corresponds to a band of momenta about 170 Mev/c. broad at a 
mean value 450 Mev/c. Their result is therefore shown in Figure 3 at this mean 
momentum and at proton intensity 0-3 x 10-§ cm? sec! sterad-1 (Mev/c.)-1. 


4.2. Comparison of Total Proton Flux with Data from Photographic Emulsions 
The vertical intensity of protons with momentum greater than 10? Mev/c. may 
be compared with the total flux of star-producing, minimum ionization, charged 
particles derived from star counts in photographic emulsions and the estimated 
value of the collision mean free path of these particles intheemulsion. Page (1950) 
deduces the total flux of star-producing particles at sea level, assuming for the 
collision cross section of the constituents of the emulsion the geometrical cross 
section of the nuclei to be 28 cm-? day-t.. The fraction of these with charged 
primaries (<1-:5 minimum ionization) has not been measured at sea level, but is 
probably not very different from the fraction 0-17 observed at 3,450 m. The flux of 
star-producing charged particles (p>0-98 x 10% ev/c.) is thus 0:55 x 10-4cm™? sect. 
Since this result is based explicitly on the assumption that substantially all 
stars are produced by strongly interacting particles, it is directly comparable 
with our figure. We may now deduce the zenith angle variation of the proton 
flux. Writing J) for the vertical intensity, /f,(@) for the intensity at zenith angle 6, 
and J for the total flux, 
7/2 
si " Tyf(O)2nsin 6 8. 
J0 
If we approximate f(?)=cos” 6, then 
J 
nly 
Thus, inserting experimental values of I) and J, n=(27I)/J)-—1~5. This result 
is, of course, a rough one; using the attenuation mean free path of the protons in 
the atmosphere deduced later (~140 gm.cm~?), the zenith angle effect arising 
purely from the change of total path length with angle is represented by n~7, 
and the true zenith angle effect is clearly to be expected to be somewhat less sharp. 


0 
=— | cos” @d(cos 0) =(n+ 1)-1. 
1 


4.3. Variation of the Vertical Proton Intensity in the Atmosphere 


The proton spectrum at 30,000 feet may be estimated in a form suitable for 
comparison with our data from results given by Adams, Anderson, Lloyd, Rau and 
Saxena (1948). Of 183 counter-controlled cloud-chamber photographs of isolated 
particles, assuming the positive—negative ratio of the relevant mesons to be 1-25, 
about 20 are probably due to protons of momentum greater than 10? mev/c. 
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(30 to protons of lower momentum and about 20 to electrons). The intensity of 
fast protons is thus about 17% of the meson intensity. We take for the vertical 
meson flux at 30,000 feet (310 gm.cm-?) 4:9 x 10-? cm sec! sterad-!, and 


. thus deduce the flux of protons of momentum greater than 10? mev/c. to be 


0:35:10 cm isecxsterad™: 

The proton intensity varies by a factor (0-85 x 10-?)/(0-55 x 10-*), that is 155, in 
traversing 720 gm.cm of air between 30,000 feet and sea level. The attenuation 
range is thus 143 gm.cm-?, with an uncertainty mainly arising from the relatively 
small number of protons observed by Anderson, and their identification by 
difference. We estimate this number to be 20+8: the range is accordingly 


pao Isegmicni*. 


It is further possible to compare the form of the proton spectrum at sea level and 
at 30,000 feet. We deduce from the data of Anderson and his co-workers the 
results of ‘Table 6, and these are plotted at arbitrary intensity in Figure 3. To the 
rather limited accuracy of the data, Anderson’s points conform to the slope of the 
sea-level spectrum, and to this extent the form of the spectrum appears to be 
unchanged between 30,000 feet and sea level. 


Table 6. Proton Spectrum at 30,000 feet derived from data of Adams et al. 
(1948) (arbitrary intensity) 


Curvature range Proton intensity 
Hx 10® gauss cm. 10? Mev/c. (arbitrary scale) 
0-5-1:5 0-3 0+3 
1:5-3-0 0-7 20+6 

3-4 1-0 8:54 
4-8 1:8 2:°5+1:5 
8-12 3-0 1:0+0°8 


An estimate of the intensity of the proton flux at 3,250 m. can be deduced from 
measurements by Alikhanian and his co-workers (1948), but since these measure- 
ments refer to particles stopped in a rather thin absorbing layer they are not 
particularly suitable for making a comparison of the form of the spectra. The 
estimated proton intensity at 3,250 m. is 7-7 x 10-4 cm™? sec"! sterad“1, an increase 
of 14:0 + 2-0 over the sea-level intensity. 


4.4. Production of Penetrating Showers by Protons 

In the foregoing sections we haye considered exclusively particles of strong 
nuclear interaction (protons) which are detected by the absence, after interaction, 
of forward-moving charged particles of range relevant to the method of detection. 
It is possible that this criterion of selection biases the observations in the sense of 
failing to identify protons which give rise, after interaction, to penetrating ionizing 
particles. This problem is discussed again in a later section: here we describe 
features of the measurements with arrangement (b) which might be capable of 
detecting penetrating showers, that is to say, occasions when two or more fast- 
charged particles result from the interaction of a proton. Since the deductions 
made here do not depend critically on the instrumental background inefficiency, 
more extensive data than those used in Table 2 are available. ‘The occurrence of 
discharges of two subgroups of counters, either in one of trays UL, BCD, EFG of 
arrangement (b) or in two trays simultaneously, is summarized in Table 7. No 


examples of double discharges simultaneously in all three trays were encountered. 
27-2 
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Column (4) of the Table shows the occurrence of double discharges in one tray only 

of the arrangement (b). Substantially all of these will be due to knock-on electrons 

produced by mesons, and in fact the percentage of double discharges per tray 

traversed varies in the way to be expected with meson momentum. We calculate 
directly from column (4) the corresponding number of simultaneous events in 

two trays, and this is tabulated in column (6) for comparison with the observed 

number of simultaneous events (column (5)); while the calculated total number of 

double counts is 13-2, 18 were actually observed. The difference is not significant, 

but, in addition, the figure 13-2 is a lower limit, for each momentum category covers 

a range of momenta, and the mesons leading to detected knock-on electrons in one 

tray will be derived preferentially from the higher momentum part of that range. 
and will thus be more likely to produce a second double discharge than is indicated 

in the Table. 


Table 7. Multiple Discharge of Counters in Arrangement (d) 


(1) (2) (3) (4) (5) (6) 
0 12 705 81 1 3-1 
| : 908 70 1 2:3 
at 739 74 6 1:9 
ae 771 66 2 
=) 3 592 40 g at 
3 543 44 1 
<3 464 26 2 ae 
S24 425 21 0 
—4 aos Rees 21 0 CP 
+ (6-10) 0-9 838 24 0 0-2 
+(>10) 0:5 130 0 0 0 
Total double discharges 18-4 13-2 


(1) Momentum category. (2) Mean momentum (10*Mev/c). (3) Total number of 
particles. (4) Double discharge of one tray only. Double discharge of two trays: 
(5) observed, (6) calculated from column (4). 


There is thus no statistical indication, to the accuracy of observation set by the 
knock-on electrons from mesons, of events attributable to showers penetrating the 
whole absorbing system arising from the stopping of protons. The complete 
absence of double discharges simultaneously in three trays is in agreement with 
this conclusion. Some evidence of multiplicity, but comprising particles of 
rather short range, is, however, indicated by records in which multiplicity is 
detected in the last counter tray which is discharged by a particle stopped in the 
absorber. 

The number of events of this type due to knock-on electrons should be one- 
twentieth of the background number of particles appearing to stop (Table 7); the 
actual record shows twice the expected number. .One negative particle from 
twenty appearing to stop gives a double count in the last tray which it triggers, in 
agreement with the knock-on rate, but six positive particles from fifty give a double 
discharge ; in two of these events a double discharge occurred in both trays UL and 


Proton Component of Cosmic Rays Aig 


BCD, and no particle triggered tray EFGH; the probability of such an event 
arising from the formation of knock-on electrons is about 1/400 of all particles 
stopping. ‘Taking into account the low efficiency with which the production of 
secondaries of moderate range will be detected, we estimate that 5-10°% of the 
particles stopped probably give at least two secondary particles of ranges of the 
order 5-10 cm. lead. 


4.5. Negative Particles stopped in Lead 


In Table 2, of about 2,000 negative particles with momentum greater than 
600 Mev/c., four are recorded as stopping in 0-10 cm. lead and none in the following 
5 cm. lead (10-15 cm.). Since we have shown that the probability of a particle 
traversing the counter system failing to discharge counter layers UL and BCD as 
well as EFGH is about 1/30 of that for failing to discharge layer EFGH alone, and 
since the apparent rate of stopping, if due to counter inefficiency in the absorption 
apparatus, should be greater by a fraction of about 5 in the lead layer 10-15 cm. 
than in the lead layer 0-10 cm., these counts cannot be considered to arise from 
inefficiency of the absorption apparatus. 

It has been shown elsewhere (Hyams et a/. 1950) that a spurious record of the 
spectrograph is probably obtained about once for each thousand particles recorded, 
and that much of this inefficiency arises from the action of extensive electron 
showers. Such records would clearly be expected, as a rule, to lead to apparent 
absorption in the first layer of lead, and probably contribute about one-half of the 
observed effect. ‘The remainder may well represent the contribution of isolated 
electrons and z-mesons; this remainder amounts to only about 0-1°% of the 
negative beam. 

According to Williams (1939) about 1°% of cosmic-ray particles of momentum 
greater than 3 x 10° ev/c. at sea level are electrons, but, to minimize the effect of 
electrons in our experiment, the particle beam entering the spectrograph is caused 
to traverse a layer of lead 1-2 cm. thick placed about 30 cm. above the top counter 
tray of the spectrograph. This layer will strongly degrade in energy electrons 
falling on the apparatus, which will also be partially excluded from the accepted 
records by the occurrence of multiple discharges in the top counter tray of the 
spectrograph. Following the cascade fluctuation calculations of Arley (1949), the 
probability of an unaccompanied electron of energy greater than 10? Mev. emerging 
from 1-2 cm. lead is of the order 0-20 ; thus much less than 0:2°% of the spectrograph 
beam is expected to be electronic. 

The probable admixture of 7-mesons in the beam may be deduced, noting that 
for proper lifetime 1-8 x 10-8 sec. the mean range in air at sea level against decay at 
the relevant momenta would be of the order 3 gm.cm™?, in which less than 5°% of 
existing nuclear interacting particles will make collisions. Since less than 5% 
of the nucleons are of sufficient energy to produce 7-mesons of the observed 
momentum, the resultant intensity of these can clearly only be of the order 0-005% 
of the totalbeam. A contribution of about the same magnitude may be expected to 
arise from collision in the 1-2 cm. lead screen above the spectrograph. 

It is thus probable that the 0-1°% fraction of the negative beam which is observed 
to be readily absorbed is mainly electronic. The measured intensity is essentially 
a function of the apparatus and not of the free electron flux at sea level, for the 
apparatus is strongly biased against recording fast electrons. 
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§5. DISCUSSION OF RESULTS: INTERPRETATION 

5.1. The Measured Range of Protons in Lead 700 Mev/c.<p<2,500 mev/c. 

The probability that protons entering the lead will give rise to detectable 
ionizing particles below the absorber is described by the range 160 + 40gm.cm “2. 
If, as is probable (see, for example, Bethe 1949), the collision cross section of 
protens in this range of momentum in lead is the geometrical cross section, 
oo = 0-063 A?? barn, leading to (1/A),,,=156 gm.cm’, it must be concluded that 
very few of these protons give rise to ionizing products of range comparable with 
the primary collision range. 

This result is clearly to be expected at the lower momenta, p~103 Mev/c. 
(Exin~420 Mev.), for the scattering cross section per nucleon extrapolated from 
the Berkeley measurements (0-013 barn per nucleon at 350 Mev.) is still large 
enough to yield at least one collision in a single heavy nucleus, and any considerable 
overall loss of energy by the primary proton will leave it with very short ioniza- 
tion range (of the order of 17 gm.cm~? at 100 Mey.). At higher momenta 
(p~2,500 Mev/c., Ey;,~1,700 Mev.) there is no clear information about either the 
total collision cross section or the cross section for meson production. 
Phenomenologically we would expect our measurement of range to be modified to 
values greater than the collision range (a) by the emission, in the collision, of protons 
of energy greater than about 200 mev., (6) by the emission of mesons of energy 
greater than about 100 Mev. The number observed, however, is quite inadequate 
to determine a separate effective range for these more energetic particles. The 
evidence available would suggest that the increase of effective range for these 
particles is not great. 

In a recent paper (Camerini et al. 1950) protons emitted from stars in a photo- 
graphic emulsion have been studied, and their range before collision is found to be 
about 100 gm.cm™ of emulsion, corresponding to a cross section in the emulsion 
only slightly less than the geometrical cross section. The range of momentum 
covered in these measurements is similar to that covered by ourselves, namely 
0-3 to 1:3 x 103 Mev/c. 

In contrast, the reduction of intensity of star production in lead takes 
place (Bernardini 1948, George and Jason 1949) with a range of greater than 
300 gm.cm™*. The difference clearly depends largely on the part played in star 
production by neutrons of energy comparable with the slowest protons which we 
measure. 


5.2. Form of the Proton Spectrum and the Attenuation in Air 


That practically the whole of the proton flux measured in our experiments must 
be of a secondary nature is established by the following results. 

(i) The work of Allen and Singer (1950) shows that the differential 
primary spectrum (2-6 x 103 Mev.<p<15 x 10° Mev.) can best be fitted by 
dN =0-53(pe/Ze)-*1d(pc/Ze). This spectrum is already much flatter than our 
spectrum, and when allowance is made for the ionization loss of energy in 
traversing the atmosphere the spectrum would become even flatter, corresponding 
to an index of about — 1-5, as compared with our value of —2:8. 

(ii) The similarity between the shape of our spectrum at sea level and that 
determined by Anderson at 30,000 feet again suggests that the particles are of a 
secondary nature, for 
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(iii) the Anderson spectrum must be composed entirely of secondary particles 
since the geomagnetic cut-off momentum at Pasadena is about 5 x 10? Mev/c. 

The attenuation of this low energy proton spectrum can be obtained most 
accurately by comparing the intensities measured by ourselves, Alikhanian, and 
Anderson at three different heights in the atmosphere. These intensities are 


0-55 x 10-* cm sec sterad under 1,030 gm.cm™® air, 
7-7 x10-4cm™ sec! sterad4 under 690 gm.cm™ air, 
85° x10 cm sect sterad under 310 gm.cm™ air. 


Assuming a relationship between the intensity and depth of the form J=J,e-, 
we obtain the attenuation range 1/A, 140+ 10 gm.cm™. 

It is rather difficult to make an adequate comparison of this range with that of 
energetic, penetrating shower-producing particles in the atmosphere or of the star- 
producing component, for both of these ranges are determined from measurements 
for which an application of the Gross transformation is necessary. The range 
which we have derived, on the other hand, to a first approximation is not subject to 
the transformation, since measurements at all levels are made with apparatus with a 
strong bias for the selection of vertical trajectories, this bias being, as is required, 
strongest for the sea-level measurements. A weak transformation should probably 
in fact be applied, since our measurements concern particles of such low energy (in 
terms of the particle rest mass) that diffusion of direction, away from the direction 
of the primary particles from which our protons are derived, is appreciable. 

The exponential range of the particles producing penetrating shower 
phenomena is 118 + 2 gm.cm-? in air (‘Tinlot 1948), and an application of the Gross 
transformation to this value, bearing in mind the limited aperture of the apparatus 
used, probably increases this value to 125 or 130gm.cm-*. The transformed 
range given by Wataghin (1949) is 119+5 gm.cm-?. The transformed range of 
star-producing particles, on the other hand, is given by Camerini and his colleagues 
(1949) as 170gm.cm™ for all stars, ranging down to 105gm.cm™ for stars 
producing six or more fast secondaries. The exponential range equivalent to 
that given by Camerini for all stars, 135 gm.cm™, is in accord with other star 
measurements (Bernardini 1948, George and Jason 1949, Page 1950). 

The attenuation of our protons in the atmosphere thus follows more closely that 
of the energetic primaries of penetrating showers than that of the star-producing 
component, although the measured protons approximate in energy to the latter. It 
is probably correct to consider the penetrating shower primaries to be the direct 
parents of our protons, and it is then necessary to determine whether any large 
difference of attenuation of the primaries and of their product protons in the 
atmosphere is to be expected. 

Assume, as a simplified model, that the primaries of our protons traverse the 
atmosphere without change of spectrum (Heitler and Janossy 1949) with 
an attenuation coefficient 4,. The secondary spectrum arising between depths 
x and x+dx of atmosphere is then I,(p, x) dp =I) exp(—A,x) S(p) dp dx, where 
I, defines the intensity of the primary beam and S(p) the resultant secondary 
spectrum arising from the undefined but constant primary spectrum. The 
observed secondary spectrum at depth X will then be 


rx 
Trotan (DP, X) dp = Lp | ee Pgs iol a(X —x)| dx dp, 
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where A, is the absorption coefficient of the secondary protons in air and « their 
ionization momentum loss in unit thickness. The total number of secondaries 
(p>po) at depth X is accordingly 


co 7X 
N,(PoX)=Iyexp (-AX) [| exp Qy—A)(X—%) S[p + a(X—a)] deedp 
Do~ 0 
and the apparent attenuation of the beam, Aq, is given by 
xy 
| i expSA(X, — x) S[p +0(X, —x)] dx dp 


Ni(Poy Xs) _ 
Ns(Po, Xo) 


> 


exp NC) 
: { | exp dA(X_—x) S[p +a(X,—x)] dx dp 
=exp [—Ay(X,—Xy)]- 
Ignoring the effect of ionization, this clearly becomes 
1 —exp X,6A 
1—exp X,d.° 
Substituting the values at sea level and 30,000 feet for X, and X, respectively 


and taking 1/A, =120gm.cm™, we get the following table of results connecting A, 
and A, : 


exp[ —A,(X, — X)] =exp[—A,(X1 — X2)] 


1/A, (gm. cm-2) 75 100 125 150 
1/A, (gm.cm-2) 125 135 155 175 


The range of the secondary particles required to give 1/A, =140gm.cm™ is 
thus a little greater than 100 gm. cm, a value which is necessarily rough but which 
is quite reasonable. Assuming the above model to give the range for star-producing 
particles also, we find that in this case the secondary particles must have a range of 
about 145 gm.cm-?. 

Since most of the low energy stars are produced by non-ionizing particles, these 
particles can be expected to have a longer range than the 100 gm.cm™ given 
above since, over the momentum range discussed, the ionization energy losses will 
reduce the effective range of charged particles. Also any particles of the third 
generation produced will, if charged, have an extremely short ionization range, 
whereas, if neutral, they will continue until they make a nuclear collision and so. 
help further to increase the effective range of star-producing particles to a value 
greater than that of the charged nucleons which we have observed. 


5.3. The Form of the Proton Spectrum at Formation 


The above analysis also allows us to estimate the formation spectrum of the 
protons for which we have given the sea-level spectrum in the form S(p)~p-?°. 
Assuming for this purpose A, =As, 


aX 
Tsou (PX)~|, [pal X— a)? dump OY (pak), 

Now «X(sea level) =2 x 10% Mev/c.: fitting to the observed spectrum at 

p=10? mev/c. and p =3 x 10° Mev/c., 


Tstotar (1 ? x) ae 


Tia Gek) 17 (experimental) and so y~3°3. 
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Thus Stormation (D) db ~ p** dp. This result can be compared with that measured 
directly by Camerini and his co-workers (1950) in photographic emulsion exposed 
at an altitude of 3,500m. Transforming Camerini’s results from a kinetic energy 
scale to a momentum scale, Spomation (p) dp ~ p21 dp; the agreement is then 
satisfactory. 
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Further Measurements of the Charge Ratio of u-Mesons at 
Sea Level 


By B. G. OWEN anp J. G. WILSON 


Physical Laboratories, University of Manchester 
Communicated by P. M.S. Blackett; MS. received 23rd November 1950 


ABSTRACT. Extended measurements of the charge ratio (positive/negative) of near- 
vertical 4-mesons at sea level show that it increases with meson momentum from about 
1:17 at 10? Mev/c. up to about 1:26 at 4x10 mev/c. At higher momenta the statistical 
accuracy of the measured charge ratio is limited by the instrumental rate of collection, and 
it cannot yet be definitely decided whether the charge ratio decreases with increasing 
momentum, although the evidence now available favours a slow decrease. 

The general implications of the form of variation of charge ratio with momentum are 
indicated. 
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§1. INTRODUCTION 

REVIOUS measurements of the charge ratio of vertical sea-level mesons 
Pos with the Manchester magnetic spectrograph have been reported by 

Owen and Wilson (1949) and Wilson (1949). A further discussion of this 
work is now given because the data available are considerably more extensive, 
and particularly because the parallel investigation of the proton component in 
the vertical beam (Mylroi and Wilson 1951) has now yielded adequate 
information for the reduction of the crude measured charge ratios at the lower 
momenta to values strictly applicable to u-mesons. 

A considerable volume of work by other groups has now been reported, and 
it is also useful to review existing information, to determine whether serious 
discrepancies exist among the various values, and to give as far as possible a 
unified picture of present knowledge. 

No important modifications of the technique described in an earlier paper 
(Owen and Wilson 1949) have been made, and the procedure of measurement, 
field reversals and other measures, will not be discussed again. In particular, 
the angular distribution of particles considered is unchanged. 


§2. EXPERIMENTAL RESULTS 
The crude ratios are summarized in Table 1. These require correction for 
two factors: the presence of the proton component, which is of appreciable 
intensity at the lower momenta, and the contamination of measurements in 
momentum category | with particles of the opposite sign. ‘The detailed results 
of Table 1 are collected into six momentum ranges and corrected accordingly. 


Table 1. The Crude Ratios of Positive to Negative Particles in the Vertical 
Beam at Sea Level 


Deflecting field Momentum Poe ie Observed no. of particles ; 
(gauss cm.) categories pore) +ve —ve are 
2 (10° mev/c.) 
2-9 x 10° 1 0-81 3248 2600 1-249 
2 0-54 1687 1388 1215 
Zor 0-46 1959 1523 1-286 
3+4* 0-30 2068 1639 1-262 
or O-11 1896 1522 1-246 
>4* 0-02 2201 1802 22d 
SHS) IMO? 1 0-85 3382 2815 1-201 
2 0-62 3204 2544 1-259 
3+4 0-35 4287 3348 1-280 
>4 0-07 4553 3639 (£233 
yas Se 1KOY 1 1:06 2624 2161 1-214 
2a 3 0:76 5544 4333 1-279 
44546 0:48 . 5178 4181 1-238 
>6 0-21 4492 3548 1-238 


* Separate sections of the records were summarized in different ways; there is no 
overlap of data. 


The proton component is taken as given in Figure 2 of the paper by Mylroi and 
Wilson (1951). The correction for high momentum contamination is not quite 
so straightforward, for it must be based on the charge ratio of the contaminating 
particles, which are of decidedly higher momentum (mean about 80 x 10? Mev/c.) 
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than the main measured group. ‘The mechanism by which the contaminating 
particles appear in category 1 measurements is explained in Figure 1 (and its 
legend). There seems little indication that the charge ratio is falling off rapidly 


(A) and Sy «l (A) 


48 32 1-6 0 16 


Figure 1. Illustrating the way in which the measured momentum category 1 includes some 
particles (about 3%) of opposite sign to the nominal sign of the group. (Notation as in 
Hyams et al. 1950.) S(A) is the deflection spectrum of mesons falling on the apparatus, 
while P(A— Anz) is the collecting efficiency function of the spectrograph drawn, relative to 
the deflection spectrum, in the position appropriate to momentum category 1. The 
deflection spectrum of particles accepted by the spectrograph, S,(A) is then given by the 
product S(A)P(A—An) and this curve will be seen to extend through A=O!to include 
deflections in the direction opposite to the nominal deflection. 


Table 2. Collected Groups of Charge Ratio Data for the Vertical Beam at 


Sea Level 

Log (mean Observed no. of Crude Corrected 

momentum) particles charge charge ratio 

(10? Mev/c.) +ve —ve ratio (u-mesons) 
1-06 2624 2161 1-214 1-229--0-036 
0-83 6630 5415 1-224 1-240+0-023 
0-69 10435 8265 1-263 1:263--0-019 
0-44 11424 9052 1-262 1-251-+-0-018 
0-28 10847 8535 1-271 1-241-40-018 
0-06 8650 7017 15733 1-172+0-019 


with increasing momentum, and none at all that it is rising. We have therefore 
made the correction on the assumption that the charge ratio does not change 
significantly between 12 10%Mev/c. and 80x 10?Mev/c., and since in fact 
a slow fall of charge ratio is, on balance, seen to be probable, the corrected ratios 
are probably very slightly too high. The collected values in their crude form and 
after correction are given in Table 2. 


§3. COMPARISON WITH OTHER OBSERVATIONS 

In Table 3 we summarize the measurements of other workers who have 
measured the charge ratio near sea level. These measurements are plotted 
together with our present results (column 5 of Table 2) in Figure 2. The early 
cloud chamber observations of Blackett (1937), Jones (1939) and Hughes (1940) 
which were quoted separately in the earlier paper are now combined. 

The majority of the results reported elsewhere are in good agreement with 
our measurements, and taking these into account according to their statistical 
weights does not significantly alter our values. It is particularly satisfactory to 
note the confirmation which we give to Nereson’s low value near 10% Mev/c. 
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This is of comparable weight to our result but is obtained by a completely 
different method. 

Some comment is necessary, however, on the measurement by Brode* and 
on the high momentum group of the values reported by Caro, Parry and 
Rathgeber. 


1640-2 
C.PR 


ae 
S 
B 13p 
So 
8 Com SOP eel Aeaie [epee | mee 
oO 

b2- 

Ba 
I 05 


Pp ( 10° Mev/C) 


Figure 2. Charge ratio of vertical u-mesons at sea level as a function of momentum. @ the results 
reported in this paper, © measurements of other workers. Co, Conversi (1949) Ba, Bassi, 
Clemente], Filosofo and Puppi (1949); N, Nereson (1948); Br, Brode (1949), C.P.R., 
Caro, Parry and Rathgeber (1950). The light horizontal broken line is the average value 
given in cloud chamber measurements. adopted form of the variation of charge ratio. 
---- obtained by including results of Caro et al. at their normal statistical weight 
(momenta >5 x 10° Mev/c.). 


Table 3. Values of Charge Ratio at Sea Level reported by other Workers 


Reference Method ie Charge ratio 
(momentum) 
Blackett 1937 ia 
Jones 1939 Cloud chamber 165133 1-225-+-0-049 
Hughes 1940 ; 
Conversi 1949 Relative decay rates in C and S LES U°24 026 
1:74 1:17 +0:13 
Bassi et al. 1949 Deflection in air 1-93 PASSES sh) 
0:25 1:29 +0-10 
Nereson 1948 Detection of decay electrons 1-95 1-170-L0-020 
(deflection in iron) 
Brode 1949 Deflection in iron 0-23 1-320-+-0-020 
1:48 0:97 +0-23 
1:78 1:24 +0:-20 
1-90 1:07 +0-14 
0 0-95 -+-0-10 
tCaro et al. 1950 Deflection in air ee ' vis an | 
0-40 1:32 0-09 | 
0-54 1:23 +0-10 | 
0:95 1-38 +0-10 | 
(cae 1:60 --0-20 | 


* We were fortunate to have an opportunity of discussing this measurement with Professor 
Brode in September 1950. 

+ The figures quoted are derived from a published graph. In particular the mean momenta 
tabulated have been estimated from the block limits of momentum indicated in the original and from 
the known meson spectrum. 
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The result quoted by Brode is obtained using an iron core deflection system 
in which, unlike the earlier apparatus of Bernardini and his co-workers (1941), 
the selecting counters are out of alignment, the measured effect therefore not 
being diluted by the slightly deflected particles of high momentum. The crude 
measurement may be deduced from data in this paper to be 1-22; this is increased 
to 1-32 by a correction based on the counting rate of the apparatus when the field 
in one half of the iron traversed by the recorded particles is reversed. All such 
counts are attributed to field-insensitive phenomena. Accordingly, this counting 
rate is subtracted from the main measurements (of either sign) and the ratio of 
the remaining counts which are accepted is increased. This correcting count, 
with deflecting fields opposed, is, however, likely to contain a field sensitive 
element arising because the scattering of particles traversing a piece of iron of the 
dimensions used is not negligible compared with the magnetic deflection. That 
part of the count with opposed fields which arises in this way should not be 
applied as a correction to the main measurements, for the probability that 
scattering will reverse the direction of deflection in both halves of the iron system 
simultaneously is negligible. It is therefore likely that Brode’s measurement is 
overcorrected, and that the charge ratio is really smaller than the quoted value. 
This view is supported by the fact that the collected cloud chamber data 
(Blackett, Jones, Hughes) while covering a large range of momenta draws a large 
fraction of particles from the region 10° Mev/c.<p<5 x 10%mev/c. The charge 
ratio so measured (1-225 + 0-049) is in excellent agreement with our measurements 
but does not fit well with any large particle group of high charge ratio. 

We therefore consider that the measurement given by Brode does not 
provide grounds for adopting a value of charge ratio, in the region 
103 mev/c.<p<2 x 10? Mev/c., different from that given in our work. 

The results (Table 3) given by Caro and his co-workers (1950) are made with 
an apparatus similar in principle to our own. They are characterized by generally 
low values of charge ratio up to about 1-5 x 10? Mev/c. and then much higher 
values at greater momenta. The statistical weight of the results is, however, 
generally rather low, and up to 4 x 10° Mev/c. there is no significant discrepancy 
with our measurements. ‘The two values at mean momenta about 9 x 10? Mev/c. 
and 35 x 10*Mmev/c. are of importance, for the second occupies a momentum 
range which we do not reach in our work, while the two together suggest a rather 
rapid increase of charge ratio with momentum continuing at least to 35 x 10? Mev/c. 

However, the detailed account of the experiments upon which these points 
are based is not yet available, and it may be that freedom from systematic errors 
has not been completely achieved in the relatively small number of particles 
recorded; we therefore again regard the curve drawn through our points as the 
best available. Even if we assume the measurements of Caro and his co-workers 
to be free from systematic error and draw the best linear variation of charge 
ratio with logp by the method of least squares, fitting to the low momentum 
part of our curve at 5 x 10° mev/c. (broken curve, Figure 2), the general form 
of the variation of charge ratio with momentum at the higher momenta is not 
greatly altered. 

§4. DISCUSSION 

Our conclusions refer to the full line of Figure 2, which, for the reasons given 
in the previous section, we regard as that most consistent with existing information. 
The increase of charge ratio from about 1:17 at p=10%Mev/c. to about 
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1-26 at p=5 x 10° Mev/c. seems well established, but although the results favour 
rather strongly a decrease of charge ratio at higher momenta, this cannot be 
regarded as decisively established. The decrease, if real, does not appear to be: 
very rapid. 

The data necessary to clarify this important region of the curve (Figure 2) 
cannot be expected quickly. To halve the statistical uncertainty of our point 
at 11 x 108 mMev/c. would require continuous operation (24 hours per day) of the 
Manchester spectrograph for one year, and for a longer period with the Melbourne 
instrument, the operating rate of which is not as great. From neither apparatus 
will results in fact be available so quickly, for other work will require the 
spectrographs to be operated for considerable periods under conditions unsuitable 
for high momentum charge ratio determinations. 

Interpretation. ‘The interpretation of the form of the measured charge ratio 
variation with momentum, in terms of the encounters at which meson formation 
takes place, has been examined by many authors* and no comparison of the 
models set up and the conclusions reached will be attempted here. It will be 
useful, however, to examine the general features which are indicated by the 
curve shown in Figure 2. 

We may assume that all nucleon-nucleon collisions at the relevant energies 
are free and independent whether the bombarding particle is a single nucleon 
or a small nucleus («-particle). The first collision of a cosmic-ray primary reaching 
the earth is then made by protons and neutrons in a proportion determined by 
the relative numbers of protons and «-particles in the primary beam (the 
proportion of heavier nuclei in the primary beam is small, and their composition 
is substantially the same as that of the «-particles). No evidence has yet been given 
of any change in the relative numbers of primary protons and «-particles as a 
function of momentum. 

Since at moderate energies the ratio of protons to neutrons in the primary 
particle beam is roughly 3/1 (80% protons, 20° «-particles), it follows, because 
the meson charge ratio arising is of the order 1-2 to 1-25, that overall, several 
charged mesons are produced for each incident primary cosmic-ray particle. 
In the simplest model we assume that the whole charge excess of the incident 
beam goes to mesons; there is then clearly an overall multiplicity of meson 
production of the order 7. ‘This multiplicity may arise either because about 
seven mesons are formed in the single nucleon-nucleon encounter of a first 
collision, or because a smaller number of mesons, perhaps one only, is formed 
in the first collision but that others are formed by subsequent collisions of one or 
both of the nucleons concerned in the first collision, and of nucleons in their 
turn arising in the later collisions. *We examine the form of the charge ratio 
variation to determine whether it distinguishes in any way between these 
alternatives. 

There is abundant evidence that nucleons formed in first and later encounters 
do take part in further meson forming collisions. We refer to these events as 
secondary collisions and in them the bombarding particles are mainly free nucleons, — 
with roughly equal numbers of protons and neutrons. It is to be observed that 
the first collision is the first nucleon-nucleon collision of a primary nucleon, and 
that it may well be that secondary collisions leading to fast mesons often take 
place within the same air nucleus as the first collision. 


* In particular, Clementel and Puppi (1948), Caldirola (1949), Bassi, Clementel, Filosofo and 
Puppi (1949), Cini and Wataghin (1950), Caldirola and Loinger (1950). 
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We therefore expect first collisions to lead to mesons with a charge ratio 
determined by the relative number of protons and neutrons in the primary 
cosmic-ray beam and by the average number of mesons formed in each first 
collision, while secondary collisions will lead to mesons with a small charge ratio 
because those mesons arise from a component containing similar numbers of 
protons and neutrons. The contribution from secondary collisions thus dilutes 
the charge ratio arising in first collisions. 

There is no reason to expect that the momentum spectrum of primary and 
secondary nucleons will be the same, and there is in fact evidence that the 
secondary nucleon spectrum contains a relatively larger number of low energy 
particles. For example, Mylroi and Wilson (1951) show that the differential 
momentum spectrum of protons between 1 and 5 x 10° mev/c. at sea level is of 
the form p-*® and a similar conclusion, for protons produced in emulsion stars 
at high altitude, may be drawn from the results of Camerini and his co-workers 
(1950). On the other hand the spectrum of primary protons does not, according 
to Winckler, Stix, Dwight and Sabin (1950), approach this index at momenta 
less than about 25 x 10° Mev/c. 

It is accordingly plausible to identify the well-established decrease of charge 
ratio at low momenta with strong dilution of the ratio by mesons formed in 
secondary collisions, and to assume that the products of secondary collisions 
become relatively less important at the high momenta. If this view is correct, 
the reduction of charge ratio at momenta greater than 5 x 10? Mev/c. must depend 
upon variations of charge rtaio in first collisions. While we cannot at present 
exclude the possibility that the reduction arises because the proportion of 
«-particles (and so of neutrons) in the primary cosmic-ray beam increases with 
momentum, it is more probable that it is variation characteristic of the detailed 
production of mesons in nucleon-nucleon collisions, and as such must be taken 
as evidence that the average number of mesons formed in a first collision increases 
with the energy of the mesons produced. 

The limitations of measurement in the spectral region for which the 
high-momentum reduction of charge ratio takes place are sufficiently serious to 
make any more detailed deductions very uncertain. 
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LETTERS FO" LAE cULLOR 
The Electron Voltaic Effect 


We have observed that photovoltaic cells are sensitive to electron bombardment. Ina 
particular series of experiments a Westinghouse selenium cell was bombarded .with a 
stream of electrons accelerated by a potential difference of up to 80 kv., after the varnish 
protecting the cell had been removed. The electrons were slightly concentrated so that 
they impinged on to a circle of 14mm. diameter and both the bombarding current and the 
circulating current were observed. In the Figure the ratio of circulating current to 
bombarding current (gain) is shown for a primary current of 1-5 x 10-§ amp., corresponding 
to about 10-§ amp/cm?. The load of the circulating current was about 700 ohms. It is 
seen that the gain has a maximum at a value of 500 in the neighbourhood of 50 kv. At 
higher current densities the gain is lower. 
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At low voltages the gain rises rapidly with voltage, thus appearing to be proportional 
to the number of ion pairs produced in the selenium. ‘The maximum: may be understood 
on the basis that the sensitive region of the cell is near the surface : with increasing voltage 
the electrons penetrate to greater depths, and the number of ion-pairs formed in the 
sensitive region decreases. 

On copper oxide cells, an effect was observed only after the oxide layer had been reduced 
in thickness by grinding and re-metallized by evaporated gold. Then a small gain was 
first observed at 60 kv.; at 90 kv. this had risen to about 150. The high threshold voltage 
confirms the view that the sensitive layer in copper oxide cells is situated at the oxide— 
copper boundary. 

The experiments are being continued. They appear to have theoretical interest and 
they indicate a new way of measuring electron beams. 


Birkbeck College Research Laboratories, W. EHRENBERG. 
University of London. Cui-Sur Lance. 
16th January 1951. R. WEsT. 
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Calculation of Oscillator Strengths for certain Electronic 
Transitions in B, and Na, 


In this Letter the oscillator strengths for certain molecular electronic transitions in 
B, and Nag are computed using the quantum mechanical one-electron approximation 
(see Bates 1949, Shull 1950, Stephenson 1951). 

The constants of the transitions are given in Table 1 (see Herzberg 1950) where R(a.) 
represents the equilibrium internuclear distance in Angstrém units and R(ay) the same 
distance in terms of the radius dy of the first Bohr orbit. ‘The wavelengths appropriate 
to the electronic transitions are also given in Table 1. 


Table 1 
Transition Xe (A.) R(a.) R (ao) 
Bs 3 2n 3271 1-589 3°00 
Nag a 6812 3-078 5+83 


Expressing these transitions in the form of molecular orbitals (Mulliken 1932, 1939), 
it is found that the boron transition reduces to the form (o42s)—(¢g2p) and the sodium 
transition to (oy3s)—(og3p). ‘The wave functions of these reduced transitions may be 
written as suitable linear combinations of the atomic wave functions in the usual way. 
‘Table 2 gives the values of the appropriate integrals (as given by Stephenson 1951). These 


: Table 2 
nt=3;Z=2-02, 2s 72-42. 
Was 40g Dy B*hg="s, 
P(3, 0,0; 3,1,0) 2-99 P2:0/0592, 115, O)ulnwsllel 
1G;0, 0; 3,1,:0). 2-26 B25, 05:05) 25.450) ogee Oc48 
S(3, 0, 0) 0-46 S(2, 0, 0) 0-52 
S(3, 1, 0) —0-15 S(2, 1, 0) —0-28 


integrals were evaluated using Slater type atomic wave functions with an effective nuclear 
charge Z and principal quantum number n* (see Slater 1930). 
Making use of the expression for the dipole moment factor D 


~ [P(, 0,0; n, 1, 0)—T(m, 0, 0; n, 1, 0)? 
a [1—S(n, 0, 0)][1—S(n, 1, 0)] 


and the relation f=304D/A., where f is the oscillator strength and A, the wavelength in a., 
we finally obtain the following results. 


D 


Table 3 
Transition de (A.) D i Lifetime (sec.) 
+ Naz IDyt+—1yt 6812 0-87 0-039 1-8 <x 10-7 
Bz 9g "Dy 3271 0-83 0:077 ZAK 1058 


+ The value given by Stephenson, Nature, 1950, 166, 191, is in error owing to the use of an 
inaccurate value for the internuclear distance. 


It is to be noted that in these calculations hybridization of the states has been neglected. 


Department of Physics, G. STEPHENSON. 
Imperial College, London S.W.7. 
24th January 1951. , 
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On the Elimination of Divergencies from Quantum 
Electrodynamics 


It has now been well established that all divergencies in any arbitrary process in quantum 
electrodynamics may be reduced to three quantities: the self-energy of the electron, the 
self-charge of the electron, and the self-energy of the photon. According to Schwinger 
(1949) and Tomonaga (Fukuda, Miyamoto, and Tomonaga 1949) the self-energy and the 
self-charge of the electron may be removed by a renormalization of the mass and the charge 
of the electron, while, following Schwinger (1949), the self-energy of the photon may be 
taken as zero on account of the gauge invariance of electrodynamics. Further, it has been 
shown by Dyson (1949) that this procedure can be applied to the collision problems in any 
approximation. ‘Thus, the work of Tomonaga, Schwinger, Dyson and others has produced 
a scheme by which all divergencies in quantum electrodynamics may be removed, and a 
finite result obtained for any physical process. In this way it has also been found possible 
to give a satisfactory explanation for the observed hydrogen-line shift and the anomalous 
magnetic moment of the electron. 

The aim of the present note is to show that the above proprariite may also be carried 
out in a somewhat different way, which involves a more symmetrical treatment of the matter 
and the radiation fields. We choose the Lagrangian density for an electrodynamic system as 


L=Ly4+L1,4+L’+L,4+L,, SIAN. (1) 
where Ly and Lp are the usual Lagrangian densities for the uncoupled electron and 


radiation fields respectively, L’ is the interaction term given by L’=A,j,, and L, and Le 
are Lorentz-invariant counter-terms given by 


=dmirf, Lester 0 ete ee eee (2) 
the constants 6m and df to be determined by an investigation of the self-energy of the 
electron and the polarization of vacuum. After obtaining the Hamiltonian density from (1) 
in the usual way, and passing over to the interaction representation by transforming away 
the parts belonging to the uncoupled electron and radiation fields, we obtain for the 
Hamiltonian density in the ‘natural’ coordinate system : 


a Sf)? 
Agi Smypoys 1 Ft ea Fat 
= — A, j,—Smyob —B8fFy?— HOBf)P?-+ (OPP +... Mat wee (3) 


For practical purposes, it is most convenient to use the Feynman—Dyson form of the 
invariant perturbation theory, and it can be shown by an argument due to Matthews (1949) 
that instead of (3) we can take the Hamiltonian density in the interaction representations as 


H=—Ajyj,—Smbp—-h8fFiy?, ln ee ee (4) 


provided we always take effectively 


/ 0A,(x) 0A,(x’) 
(P( )),- 


OX 0x3" 


. (P(Ay(x), Ag(x’)))p. ws + =e - (5) 


In order to prove that the present approach will remove all divergencies from the 
S-matrix in quantum electrodynamics, we can follow the treatment of Dyson (1949) with 
the only difference that the factor Z, should now be put equal to unity, for the contribution 
of the last term in (4) for a photon of momentum k, is given by 


—6fA,(—k)[Syyk?—Ryuky|Ay(h), 555.04.0 (9) 
which cancels the entire divergency in the ‘photon self-energy part’ after the requirements 
of gauge invariance have been satisfied. Moreover, since the contribution of Z, is cancelled 
by that of Z, (Ward 1950), no charge renormalization is necessary in the present treatment. 

In conclusion, we add a few remarks about the correspondence in the elimination of 
divergencies from classical and quantum electrodynamics. It is well known that in order 


to carry out a mass renormalization, it is necessary to show that the self-energy W of an 
electron with. niomentum p can be expressed in the form 


W=dm{m/(m? +p), se eee) 
where 6m is a Lorentz-invariant quantity. Since it does not seem possible to express the 


classical self-energy of the electron in the form (7), a mass renormalization may not be 
possible in classical electrodynamics. However, this does not matter, for the mass 


OX, Chane 
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t 
renormalization may not be taken too literally, but we may regard this procedure simply as 
a covariant subtraction of the physically meaningless self-energy of an isolated electron. 
It can then be shown that a covariant subtraction of the self-energy of the electron is also 
possible in classical electrodynamics by means of a suitable counter-term (Gupta 1951). 
I would like to express my thanks to Dr. N. Kemmer and Mr. A. Salam for helpful 
discussions. 


Cavendish Laboratory, “Suraj N. Gupta.* 
Cambridge. : 
22nd January 1951. 


* State Scholar of the Government of India. 
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On the Elimination of Divergencies from Classical 
Electrodynamics 


In a recent paper with the above title Gupta (1951) obtains for the force density at the 
world-line of the rth electron in an ati of 7 electrons the expression 


POA toi ic Aw — lees: (1) 


where F,, is the sum of the retarded etd of all the electrons (including the 7th), ae sym 1S 
the mean of the retarded and advanced fields of the rth electron, and j is the current density 
due to the 7th electron, all quantities being evaluated at the world-line of the 7th electron. 
From (1) Gupta derives an equation of motion for the electron which, he observes, is just 
that obtained by Dirac (1938). 

This note is to point out that the equivalence with Dirac’s theory can be seen directly 
from (1), bata wees deriving the eee. of motion. 'To do this we write F,, as 


Ferret 2 fiw, ®t and SO om as Ef 0) ret tf aav. Substitution in (1) then gives 


petty 
eae ar Afi, ret fa, adv) + © Fy, whi [torte meelibhias! on (2). 


The first term in the braces is recognizable as the Dirac self-field of the vth electron, and 
the second is the sum of the retarded fields of all the other electrons at the world-line of the 
rth electron; in other words (2) is exactly the force density at the world-line given by 
Dirac’s theory. This argument is unaffected by the inclusion in F,, of an arbitrary external 
field which Gupta denotes by fy», ¢x and which for simplicity has here been taken as zero. 
The procedure of substracting out the symmetrical field (which diverges quadratically at 
the world-line of the electron) to obtain a finite self-field is not of course limited to Dirac’s 
theory. Eliezer (1947) has reviewed a number of theories which retain this feature, while 
modifying the self-field in the attempt to obtain an equation of motion for the electron which 
does not give the non-physical behaviour of Dirac’s equation. However, Gupta’s choice 
of the form of F,,, ensures that his formulation corresponds unambiguously with Dirac’s 
theory; for in the alternative theories the force exerted by an electron on other charges in its 
field would be determined by a mixture of its retarded and advanced fields. Gupta’s view 
that this force must not involve the advanced field, since the latter cannot be a measurable 
quantity, leads immediately to his adoption of the Lorentz—Dirac field tensor for Fy». 


School of Theoretical Physics, J. W. GarDNER. 
Dublin Institute for Advanced Studies. 
5th February 1951. 
Dirac, P. A. M., 1938, Proc. Roy, Soc. A, 167, 148. 
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Radioactive *’Dy 


In an earlier paper (Butement 1950) the preparation of a long-lived radioactive “*Dy 
was mentioned. This letter describes a detailed investigation of this isotope. 

The stable isotopes of dysprosium are given in the Table (Inghram, Hayden and Hess 
1949): 

Mass number 156 158 160 +=: 161 162 163 164 
Abundance (%) 0:0524 0:0902 2-294 18-88 25:53 24-97 28-18 

The shorter-lived activities produced by neutron capture are 140-minute ™°Dy and, 
by capture of a second neutron by this isotope, some 82-hour 1®*Dy. 

A sample of 10 mg. of dysprosium oxide was irradiated in a neutron flux of 104% cm~*sec™* 
for 20 days and the dysprosium purified by means of an ion-exchange resin column. After 
waiting 60 days for complete decay of the 1**Dy, the decay of the long-lived activity was 
followed for 400 days, and showed a half-life of 132 days (see Figure). 
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Decay of **Dy (a) from pile irradiated dysprosium, and (b) from deuteron irradiated terbium. 


A sample of 14mg. of very pure terbium oxide was given 70 microampere-hours 
irradiation with 9 Mev. deuterons, and the ®*®Dy produced by the reaction ®Tb (d, 2n) 
isolated by the ion-exchange method. Its activity decayed with a half-life of 136 days 
(see Figure). 

Both the above dysprosium samples emitted only electromagnetic radiations with 
energies (determined by absorption measurements in copper and aluminium) of 6-6 and 
44 kev., corresponding to the L and K x-rays of terbium. Positron emission, if any, was 
less than one per thousand disintegrations. It is concluded that the only long-lived product 
of neutron capture by dysprosium is *®Dy. Ketelle (1949) has independently found that 
neutron irradiation of dysprosium gives a 140-day activity, which is evidently this isotope. 

The author thanks Mr. R. L. Evans for assistance with the counting, and Professor 
J. Rotblat for arranging the cyclotron irradiation. 


Atomic Energy Research Establishment, .F. D. S. BUTEMENT. 
Harwell, Didcot, Berks. 
23rd January 1951. 


BUTEMENT, F. D. S., 1950, Proc. Phys. Soc. A, $3, 532. 
INcHRAM, M. G., Haypen, R. J., and Hess, D. C., 1949, Phys. Rev., 75, 693. 
KErTELLE, B. H., 1949, Phys. Rev., 76, 1256. 


Isotope Effect in Hydrocarbon Flame Bands 


Though there have been reasonable grounds to believe that the hydrocarbon flame 
bands have to be attributed to HCO, there has been so far no direct proof. An attempt 
has now been made to determine the emitter by utilizing isotope effect through deuterium 
substitution. "The source used for obtaining the isotope shifts was the glow of deutero- 
acetylene interacting with atomic oxygen at a pressure of 3 mm. Hg. 
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In the original analysis (Vaidya 1939) the hydrocarbon flame bands were divided into 
three groups, A, B, A’, A”. The A group contained most of the bands, while the B bands 
appeared under special circumstances. The v’, v” numbering of the A bands was 
provisional. The stronger bands belonged to v’=0 and v’=1 progressions. 

In Table 1 are shown the isotopic bands and the shifts. From the order of diminution 
of the shifts, the position which would be occupied by the (0,0) band was determined. 

An independent check on the location of the (0,0) band was made by giving longer 
exposures in order to determine the extent of the spectrum in the ultra-violet. This 
resulted in better photographs of the spectrum near A2550 and from a careful examination 
of the spectrum it was concluded that the (0, 0) band was missing. That the system did 
not extend farther was also concluded from the absence of the band which was expected 
to appear at v’=2 and v’=0 on the Condon intensity parabola. 

Thus, by fixing the position of the (0, 0) band, the vibrational constants of the A system 
could be determined. The A bands are represented by the formula 


v= 39849 + [1075(v’ +4) —6(v’ +4)?]— [1883(v" + 3) —9-4(v” +4)? ]. 
The isotopic bands can be represented by 
vi = 39843 +-[1062(v’ + 4)—5-7(v’ + 4)?] —[1855(v” + 4)?—8-7(v" + 8)? ]. 


There were two sets of extra bands which did not belong to the main A system; they 
were denoted B bands and A’, A” bands. By locating shifts corresponding to these it has 
been possible to understand their significance. Unfortunately, owing to overlapping, 
observations were possible on only a fewofthem. The isotopic bands of these two sets are 
shown in Table 2. 


Table 1 
Ordinary Isotopic Shifts Old New 
(cm~) (cm™*) (cm7*) (v’,v") (a, vo”) 
v 39476 0,3 0, 0 missing 
37600 37630 30 0,4 al 
35741 35793 by, 0,5 0, 2 
33909 33978 71 0, 6 0, 3 
* * * 
30302 30423 121 0,8 ONS 
28541 28682 141 0,9 0, 6 
26798 26955 157 0, 10 0,7 
25071 25249 178 sah 0,8 
Table 2 
B bands Associated A bands 
Ordinary Isotopic Shift Normal Isotopic Shift 
(cm7") (cm) (cm~*) (cm™) (cm™") (cm7*) 
vy 29762 29892 130 29500 29734 134 
26290 26447 5/7 26137 26298 161 
A’, A” bands 
Ordinary Isotopic Shift 
(cm~*) (cm™*) (cm7*) 
AGE 30333 36420 87 
A’ 27498 27764 266 


The near equality of the shifts of the B bands and the associated A bands suggests that 
the B bands are the R heads. The shifts of the A’ and A” bands are, however, much larger 
than for the A bands in the same region. . ; 

In extending the spectrum in the ultra-violet a set of new bands appeared which have 
been interpreted as being due to the bending frequency of the emitter. Their wavelength 
and the quantum numbers are given below. 

The small change in the vibrational frequency of the A band system suggests that though 
H is present in the emitter, the vibration which gives rise to the hydrocarbon flame bands 
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is such that H is not directly taking part in it. As regards the likely nature of the other 
vibration concerned in the emission of the band system, the recent experiments of Dyne and 
Style (1947), where they obtained the A bands by fluorescent excitation of formaldehyde, 
suggest that the vibrator must be C=O. Thus, the emitter has to be HC=O. The 
other combination C—OH is very unlikely since the vibrational frequency of single bond 
C—O is very much smaller. 


The discovery of bands corresponding to the bending frequency whose value is of the 


right order, assuming that the emitter is HCO, also supports the assignment to this emitter. 


Table 3 


Quantum numbers 


Ss 
3 


A v U1 Vo V3 Uy Ve 
2420-5(1) 41301 
2436-5(2) 41030 
2447-8(1) 40841 
2463 -8(1) 40575 
2480:0(1) 40310 


is 
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In an ordinary diatomic molecular spectrum, one can calculate the mass factor (p—1) 
from the isotope effect and this enables the reduced mass of the emitter to be determined. 
A similar calculation has been made assuming that HC moves as a unit against O; this 
gives (p—1)=—0-015, using ground frequencies while theoretically it is —0-0198. 

Further details will be published in due course. 


National Physical Laboratory, W. M. Varpya. 
New Delhi, India. 
17th January 1951.- 


Dyng, P. S., and STyLe, D. W. G., 1947, Disc. Faraday Soc., 2, 159. 
Varipya, W. M., 1934, Proc. Roy. Soc. A, 147, 513. 


Emission Band of the Cd, Molecule at 12212 


Certain diffuse bands and continua in the visible and near ultra-violet regions, ascribed 
to the Cd, molecule, have been variously observed in emission, absorption and fluorescence 
by a number of workers. Relevant data have been summarized by Cram (1934), who made 
new observations and interpreted the collated material in terms of the potential energy 
curves to be expected when account is taken of van der Waals interactions. It is worth 
mentioning that Cram acknowledges that he has taken the general form of the potential 
energy curves from theoretical predictions of Van Vleck, to appear.in a later paper of the 
latter; the paper concerned is by King and Van Vleck (1939), and in it Cram’s work is 
cited as lending support to the theory. 

Important features of Cram’s interpretation depend on the behaviour of a diffuse band 
at about A2212. This band had been observed to occur readily in absorption by several 
workers, notably Mohler and Moore (1927) and Winans (1929). An emission band having 
the same location was reported by Robertson (1932) as present when cadmium vapour at 
low pressure was excited in a ring discharge. Both Winans (1929) and Cram (1934) 
reported the absence of the band in emission, using Tesla coil excitation. Cram considers 
Robertson’s band as probably due to an impurity. It is obvious that in interpreting a 
spectrum of this kind importance attaches to the case of a band appearing strongly in 
absorption but absent in emission. 

Recently, plates exposed by members of an undergraduate class at Imperial College were 
being examined. On a picture of the spectrum of the cadmium arc in air, obtained on an 
Ilford Q/1 plate with the Hilger medium quartz spectrograph, there was noticed a pro- 
minent diffuse emission feature in the 2200 A. region. Investigation has shown this band 


| 
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to be centred on 2212 a., and subsequent exposures have proved that it appears regularly. 
On the reproduction of a typical plate, here shown, the A2212 emission band is plainly 
visible at XX. The strongly reversed singlet resonance line at 42288 can also be seen with 
its associated continuous emission band. The A2212 band is unmistakably separated 
from the latter. There seems little doubt that it is identical with Robertson’s band. 
Examination of the plate shows impurity lines to be few and weak, the strongest being the 
Cu doublet AA3274, 3247, and a few feeble Pb lines. 
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Ultra-violet spectrum of cadmium arc. 


That the A2212 band has not been previously reported in these conditions is perhaps in 
part due to the relatively infrequent use of Ilford Q emulsions with medium-quartz spectro- 
graphs. ‘The failure of Cram to observe the band in emission may be partly due to the fact 
that he was working at high vapour pressures, at which it is apparent in Robertson’s 
pictures that the band had disappeared, possibly either through change in excitation conditions 
or re-absorption in the vapour, or both. 

It thus appears that the 42212 emission band is readily excited in emission in cadmium 
vapour, and Cram’s interpretations need reconsideration, perhaps in the light of further 
experiments. y 


Department of Astrophysics, W. R: S. Garton 
Imperial College, London S.W.7. 
9th February 1951. 


Cram, S. W., 1934, Phys. Rev., 46, 205. 

Kine, G. W., and VAN VLEck, J. H., 1939, Phys. Rev., 55, 1165. 
Mouter, F. L., and Moors, H. R., 1927, 7. Opt. Soc. Amer., 15, 74. 
ROBERTSON, J. K., 1932, Phil. Mag., 14, 795. 

Winans, J. G., 1929, Phil. Mag., 7, 555. 


Free Energy of the Double Layers of Two Plates at Large 
Separations 


In a recent paper (Levine and Suddaby 1951) we developed an expansion for the 
interaction of the electric double layers of two parallel (colloidal) plates at large separations, 
when immersed in an electrolyte of the binary symmetrical type. In this note we shall obtain 
the next term in this expansion. Although the derivation of this term is very laborious, 
the final result is not complicated, and is suitable for numerical computations. Let R be 
the separation of the plates, the density of positive (or negative) ions in the interior of the 
electrolyte, v the valency of a positive ion, D the dielectric constant of the electrolyte 
(aqueous) medium, ¢ the electronic charge, k Boltzmann’s constant, T the temperature, 
kK? = 87nve?/DkT the Debye-Huckel parameter, ~m and % the potentials at the median 
plane (between the plates) and at the plate surfaces respectively, y)=eo/RT, Ym=€{m/ RT, 
k=exp(—vnm) and sin d5=k? exp(—}vum). Then the solution of the Poisson—Boltz- 
mann equation yields the relation 

4KR=k!!2{B(k, 7/2)—F(R, Po}, = = §  cecees (1) 
where F’(k, 4p) is the elliptic integral of the first kind. ‘The free energy 2F of the double 
layers is obtained from the force equation 
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assuming that the surface potential jp is independent of R. 
We first expand Fk, $0) in a Taylor series about k'2=0, where k’2=1—k?, namely 


F(k, $9) = —$ In y—}R? In y+ —6- 9 In y)k’ 


= 5 In y)k’6+ ——— (3356+ 6u—1225 Iny)R'8+..., «see (3) 


1 
(128) 
cosh 4ump ee 5 


where y=tanh fvy, d= Further details concerning the 


sinh? Zum sinh? duno 
method of determining this power series in k’ are given in Appendix I of the earlier paper 
(Levine and Suddaby 1951). The corresponding expansion of F(k, 7/2) is given, for 
example, by Jahnke and Emde (Tables of Functions, 3rd edition, p. 73). Making use of 
these series and adopting the method described in the above-mentioned Appendix, we find 
that (1) can be written in the convenient form 
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an analytic function of k’?, and 
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X=Xo exp | & Re (1.0 : kit... : In x | » Xow 1Oy exp.(—4cR)o sie nary (5) 


The relation (4) can be inverted by Lagrange’s theorem, so that k” is obtained as a power 
series in x. Substituting this latter series into (5), we use the method of successive 
approximations to ee X as aserieS iN Yo. In this way we derive that 


R= y9— 4X0? + <= (22 —KR—28) x? 4 an 6 +«KR+28)x0! 


i 


ay ane {359+ 5(kR)?+ (208 —135)«R—2708+ 2082—12p}xyg5 +... ..- eee (6) 


Finally, we expand (2) as a power series in k”, substitute (6) for k” and then integrate (2) 
with respect to R, obtaining for the interaction 


V(R)=2F—2F,,=ViR\1+@it@at.--), sen aee (7) 
where* VR LRET RAMEN Les WO eye eee (8) 
G@.=- a VER = oho) So We ee (9) 

Cs= 5 ae Xo!{3(KR)?+ 3 R(43— 8) +3 158+ 128°— 61} 
iy ye = ORS a (LelzAs0)\dves eden nae A a ee (10) 


and 2F,, is the free energy at infinite separation. We have checked the expression (10) 
for @, by expanding 2F in a double series in powers of m9? and exp (— 4«R) and comparing 
with the same double series which can be derived from the formula for 2F at moderate 
potentials, as given by Levine and Suddaby (1951). A different method of determining (7) 
will be published later. 

In the Table, we compute the values of @, and @,» for various vy) and «R and also 
the values of «R at which @,=0 and @,=0. (The equation @,=0 is quadratic in«R, but 
we need only consider the larger root.) For «R=2 the three terms in (7) yield an accuracy 
of about 1 or 2% for all vj. When vn)<5, say, then we can still use (7) for values of xR 
as low as 1°5, the maximum error at kR=1°5 being about 4%. However, when «R=1-0, 


*We wrote @, as © (vm, «R) in the earlier paper. 
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(7) only gives a qualitative answer, and for still smaller «R it is of little use. Indeed, the 
question of convergence arises at small «R. According to Verwey and Overbeek (1948), 
the region «R >2 is the important one in determining stability properties of colloids, and 


> €. 1 €» —-€1 2 —-€@1 €. €1=0 €2=0 
K 
PB 0:5 0-5 1 1 1:5 1:5 2:0 2:0 KR KR 
UNo 

QO —0-607 0-368 —0-368 Ons 5a 022511 0:0498 0:1353 0-0183 — — 
0-5 —0-588 0-345 —0-:362 07135 0-2231 0:0527 0:1374 0:0206 — — 
1:0 —0-531 0-275 —0-344, 0-131; 0:2223 0:0603 0:1430 0-0268 -— — 
1-5 —0-440 0-158 —0-314 0-117 00-2193 0:0683 0-1504 0-0356 — — 
2:0 —0-:320 0-003 —0-273 0-087 0-2129 0:0721 0-1581 0-0449 — 0-488 
2:5 —0-:176 —0-:171 —0-220 0-036; 0:2019, 0:0689 0:1641 0:0527 0:028 0:844 
3-0 —0-:018 —0-339 —0-160 —0:030 0-1868 0:0562 0:1679 0:0573 0:462 1-102 
335 0-143 —0-477 —0-:095 —0:106 0-1683 0:0355 0:1692 00-0584 0-739 1-293 
4-0 0-301 —0-569 —0-031 —0-183 0-1480, 0:0093 0-1683 0-0563 0-928 1-456; 
5-0 0-580 —0-609 0:087 —0:316 0-1076 —0-:0472 0-1626 0-:0460 1-165 1-671 
6-0 0-794 —0-519 0-180; —0:408; 0:0733, —0:0966 0:1554 0:0334 1:299 1-801 
7-0 0:945 —0-387 0:247 —0-464; 0:0479 —0:1330 0-1489, 0:0225 1-379 1:879 
8-0 1-045 —0-266 0-292 —0:-497 0-0304 —0-:1576 0:1442 0:0146 1-427 1:927 
9-0 1:109 —0-174 0-321 —0:515 0-:0190 —0-1734 0-1410 00-0092 1-456 1-956 
10-0 1-149 —0-111 0-339 —0:526 0-:0117 —0-1833 0-1388 0:0057 1-473 1-973 
ere) HDS 0 0-368 —0-541 0 —0:1991 0-1353 0 1-500 2-000 


therefore series (7) should be quite useful. 


Department of Mathematics, 
The University, Manchester. 


Sir John Cass College, London. 
13th February 1951. 


In a later paper we shall apply it when developing 
further the incorporation of the Stern theory of the double layer into a theory of interaction 
of colloidal plates (see Levine 1951). 


LEVINE, S., 1951, 7. Colloid Science, in the press. 
LEVINE, S., and SuppaBy, A., 1951, Proc. Phys. Soc. A, 64, 1. 


Verwey, E. J. W., and OverBEEK, J. Tu. G., 1948, Theory of the Stability of Lyophobic Colloids 
(Amsterdam and New York: Elsevier), p. 129. 


REVIEWS 


The Meaning of Relativity, by ALBERT EINSTEIN. 


Appendix on 


the 


OF BOOKS 


S. LEVINE. 


A. SUDDABY. 


4th Edition, containing an 
Pp. v +145. 


(London: Methuen and Co., Ltd., 1950.) 7s. 6d. 


This well-known book, originally published in 1922, is based on the text of Einstein’s 
In 1946 the 


Stafford Little Lectures delivered at Princeton University in May 1921. 


ee } 


author added an appendix ‘ On the ‘‘ Cosmologic Problem 
since 1921, and in particular rejected the cosmological constant, which he had himself 


originally introduced in 1917. 


appendix was made by Sonja Bargmann. 


Generalized Theory of Gravitation. 


in which he reviewed advances 


The present edition is noteworthy for a second appendix 
on the author’s new generalization of General Relativity, and it is to this appendix which 
most readers will immediately turn when they obtain this edition. 


The translation of this 


Newtonian physics was founded on the associated concepts of ‘particle’ and ‘action 


at a distance’, concepts which until recently were considered fully adequate for a theory 
of gravitational phenomena, although in order to explain electromagnetic phenomena 
Faraday and Maxwell introduced the ‘ field’ concept. Einstein, who has more than once 
paid tribute to the genius of Maxwell, has striven to extend the field concept so that it 
embraces both gravitation and electromagnetic phenomena in a ‘ unified field theory’. 
A first version was published by him in 1929, and this latest theory may be regarded as the 
fruit of thirty years’ reflection by the most famous mathematical physicist of the century. 
For this reason alone, it must be read with respect even by the most critical. 
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At the outset a fundamental limitation of the new theory will strike most readers, for 
it is clear that in no sense is it a quantum theory. The most that the author sets out to 
achieve is a new unification of macroscopic phenomena. It is too soon to say whether any 
new experimentally verifiable effects will be predicted by the new theory. The General 
Theory of 1916, however, was much less fruitful in this respect than the Special Theory 
of 1905. Of the three famous tests of the General Theory only the first is now regarded as 
fully confirmed, and then strictly speaking only for a single planet. Despite the startling 
agreement with observation in the case of Mercury, the appeal of the General Theory is 
primarily epistemological and aesthetic. The derivation of a more accurate law of 
gravitation than that of Newton from axioms of ‘minimum arbitrariness’ was in itself a 
major intellectual achievement. That the whole of physics can likewise be contained in 
some such compact and ‘ epistemologically primitive’ form is at best still a dream and 
some may think a mirage, but any extension of General Relativity which covers a wider 
class of phenomena without introducing more than the bare minimum of new arbitrariness 
would be full justification for the time which Einstein has devoted to the problem. 

The core of General Relativity is comprised in the standard equation relating the tensor 
&iz, which represents the metric properties of space-time, and the material-energy tensor 
T;. In this theory the tensor gj is symmetrical, and Einstein’s new theory is his answer 
to the problem: to find a natural generalization of this tensor and an associated system of 
field equations which will be the most logically satisfactory generalization of the equations 
of pure gravitation, the only clue from experience being ‘ the vague perception that some- 
thing like Maxwell’s electromagnetic field has to be contained within the total field ’. 

The main steps in Einstein’s answer to this problem may be summarized very briefly 
as follows. First, the symmetrical gj, (ten components, corresponding to the gravitational 
field) is replaced by a general non-symmetrical g;, with six additional components (corres- 
ponding to the electromagnetic field). It is then shown how an absolute differential 
calculus can be constructed for such a non-symmetrical field. The notion of affine parallel 
displacement is introduced into the new theory, and by means of it a general curvature 
tensor is defined and a contracted form Rj, obtained. The analysis gives rise to three 
Hermitian tensors which can be combined to form a new tensor Uj with certain special 
properties. A variational principle d{[g*U;,57}=0 is then shown to imply a system of 
equations (laws) ‘of utmost simplicity’, including Ry,=0, which may be regarded as 
“the natural extension’ of the equation G;,=0 of (pure gravitational) General Relativity. 
In so far as the new laws can be compared to Maxwell’s electrodynamics, they involve a 
special postulate which corresponds to. the assertion that the magnetic current density 
vanishes, but they automatically imply ‘ conservation of electricity’. Although Einstein 
believes that his latest solution of the unified field theory problem is ‘ highly convincing ’, 
he admits that formidable mathematical difficulties will have to be overcome before it is 
possible to consider testing the theory experimentally. G. J. WHITROW. 


Classical Mechanics, by HERBERT GOLDSTEIN. Pp. 399. (Cambridge, Mass. : 
Addison-Wesley Press Inc., 1950.) $6.50. 


There has long been a need for a textbook of classical mechanics written with the needs 
of the physicist in mind. The book under review tulfils this need admirably. 

The field covered is in the main the familiar one of the mechanics of particles and of 
rigid bodies, of tre Lagrange and Hamiltonian equations. Particular attention is paid to 
the latter, three chapters being devoted to Hamiltonian equations, canonical transformations 
and Hamiltonian—Jacobi theory. What is to some extent a novel feature is a chapter on 
relativistic mechanics and on the Hamiltonian and Lagrangian formulations for continuous 
systems and for fields. Oscillations of a system are dealt with in some detail, including the 
transformation to principal axes, normal modes, and the calculation of eigenfrequencies. 

Some acquaintance with vector algebra is assumed, but special sections are devoted to the 
explanation of the properties of tensors, dyadics and matrices, of which full use is made. 
The illustrative examples, both in the text and at the ends of chapters, are in the main chosen 
from physics instead of being, to quote a phrase of the author, “ the usual pedantic museum 
pieces”’. For instance, the normal modes of a triatomic molecule are used to illustrate the 
theory of small oscillations, and the transformation to laboratory coordinates of the scattering 
cross section is discussed in connection with the theory of the collision of particles. An 
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excellent set of references is included at the end of each chapter, with comments on these 
textbooks. This is a novel and rather intriguing feature. The comments are occasionally 
pungent—“ a formidable treatise which often manages somehow to make the elegant 
simplicity of vector and tensor methods appear quite complicated and repellant’”’. As 
a contribution to this process of brightening scientific textbooks the reviewer would like to 
draw attention to a comment of Tawney (Religion and the Rise of Capitalism) on a mediaeval 
historical treatise—‘‘ a comprehensive but completely unreadable account of the subject’’. 

The book is written in a pleasant style and is a pleasure to read. Occasionally the author 
tends to assume a surprising knowledge of theoretical physics on the part of the reader, though 
this usually occur's in ‘ asides’ or footnotes. Very few students will appreciate the point of 
the remark that the cross section for Coulomb scattering is infinite in quantum mechanics, 
nor, by the way, is the author correct in stating that “‘ all fields that fall off faster than a 
Coulomb field produce a finite scattering cross section in quantum mechanics ”’. 

There is hardly any need to add that all physicists are strongly recommended to lose no 
time in acquiring this textbook. M. BLACKMAN. 


Practical Applications of Spectrum Analysis, by H. Dincie. Pp. x+245. 
(London: Chapman and Hall, 1950.) 40s. 


The main theme is qualitative chemical analysis using arc spectra. On the instrumental 
side the general principles of prismatic and grating spectrographs are given, and the various 
factors influencing dispersion, resolution and light-gathering power are discussed. Methods 
of taking comparison spectra and measuring photographic plates for the determination of 
wavelengths are dealt with. “This is all quite elementary, and intended for readers without 
previous spectroscopic training. ‘The principles underlying the excitation of spectra and 
the various factors influencing the strength of spectrum lines are similarly dealt with in an 
elementary way. 

There are two chapters on qualitative analysis, mainly using arc spectra, in which much 
valuable information gained during thirty years of experience in this type of work at Imperial 
College, London, is included. This is the most valuable part of the book. The most 
sensitive lines of the various elements are considered in detail, and especially the influence 
of other possible constituents in masking or interfering with the sensitive lines. It is stressed 
that lines of the same multiplet usually vary in strength in a similar way, and the various 
multiplets among the sensitive lines are grouped so that they may be picked out for study 
together. The need for careful criticism of results before final acceptance is stressed. In 
Professor Dingle’s words, “‘ It is no exaggeration to say that in my experience more than half 
of the spectrographic analyses which have been published during the last thirty years aie 
untrustworthy ”’. Application of the principles and detailed information supplied here should 
help to prevent a continuation ot such a state or affairs. 

One brief chapter is devoted to quantitative work, mostly using the Foster arc. Professor 
Dingle appears to think that this small amount of space given to quantitative work reflects 
the present state of achievement in this field. There will probably be many who will differ 
from him. 

The book contains useful tabular matter on the sensitive lines of the elements, on Burns 
wavelengths of iron arc lines, and wavelengths of lines from a rates ultimes mixture containing 
practically all the spectroscopically sensitive elements except rare earths. ‘There are 
several good plates, including a set of the spectrum of the raies ultimes mixture. ‘The 
book is written in Professor Dingle’s usual clear style, and is excellently produced, especially 
the plates; in some of the tables there is a little difficulty in following along the table from 
the line to its identification. 

It should be clearly understood that this book does not give a comprehensive account of 
the subject. It is essentially a record of personal experience. No mention is made of 
recent developments in quantitative analysis elsewhere in this country, and in Germany and 
America. 'The new photoelectric methods using photon multipliers are not described, and 
subjects such as the quantitative estimation of vitamins by absorption spectra are omitted. 
‘One may even doubt whether the title has been well chosen, for it is not the applications at 
all but a single method, that of the arc, which is the theme. However, within these limits 
the book does give a new approach to the subject, and includes so much valuable experience 
that it may be recommended. A. G. G. 
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An Introduction to Luminescence of Solids, by H. W. LEVERENZ. Pp. xv+569. 
1st Edition. (New York: John Wiley and Sons; London: Chapman and 
Fall, 1950.) 90s: 


This book is published in the ‘ Structure of Matter ’ series, and the author’s intention is 
to provide ‘‘ an orderly introduction to a growing field’ for scientists of other interests. 
He also attempts to systematize the terminology ot the subject and to provide a general 
coordination of luminescence phenomen:, using specific materials as examples. The book 
begins with two introductory chapters, the first of which covers the elementary concepts of 
modern physics but includes rather irrelevant topics, such as Dirac’s theory of positrons and 
pair production. Chapter II is more useful as it provides an introduction to the crystalline 
state with some specific reference to the crystal structure of typical luminescent solids. 
Chapter III is concerned with syntheses of well-known types of phosphors and the system of 
notation for defining the final products. Chapter IV discusses phosphor constitution, 
crystal structure, impurity activators and the electronic energy states of phosphors ; it 
includes a very brief discussion of photoconductivity. 

Chapter V consists of more than two hundred text pages and covers all the physical 
phenomena encountered in the study of luminescence in solids. The space given to each 
topicis rather inadequate, and better treatment of such important aspects of the subject might 
have been possible by affording them the space occupied elsewhere in the text by more 
speculative discussions. It will be difficult for those’‘ not skilled in the art’ to distinguish 
the processes of luminescence treated in this chapter with respect to their relative importance, 
as the author’s style is not very lucid. The next chapter provides a summary of 
phosphor characteristics, or, as the author terms it, ‘a capsulized description of phosphors’, 
while the last chapter is concerned with the various applications of phosphors. Five 
appendices and a glossary of symbols follow. Only the first appendix, dealing with the 
preparation of pure zinc and cadmium sulphides and selenides, and the third, giving useful 
conversion factors for photometric units, appear to be of definite usefulness. Some six 
hundred text references and a supplementary list of references to very recent publications 
are given which, apart from one or two important omissions, provide a valuable guide to the 
literature of luminescence since the titles of individual papers are included. The book 
concludes with adequate formulae, author and subject indexes. 

Taken as a whole, the book contains a large amount of information on phosphors, 
particularly those of commercial interest, but its arrangement does not appear to be the best 
that could have been chosen. The reader may well object to the many cumbersome phrases 
suggestive of the jargon of luminescence laboratories but not helpful to the average scientific 
worker. Such phrases as ‘consciously added activator’, ‘difficultly eradicable complex 
absorptions ’ and ‘ energy transmittal in solids’ may be included with the above. For the 
industrial scientist the book may be of considerable interest as it represents the contribution 
of a recognized worker in applications of luminescence who has set very high standards of 
purity and control in the production, and faithful reproduction, of phosphors of high quality 
to meet the stringent demands of commercial users. G. F, J. GARLICK. 


Photons and Electrons, by K. H. Sprinc. Pp. 108. 1st Edition. (London: 
Methuen and Co., Ltd., 1950.) 7s. 6d. 


‘This monograph gives a simple account of the theoretical results which have been 
obtained for the main high energy processes involving photons and electrons, viz. the 
photoelectric effect, Compton scattering, bremsstrahlung, and pair production. The 
theoretical cross sections are in every case compared with the experimental values, and 
particular emphasis is placed on the angular distributions of the various outgoing particles. 
The last chapter contains an elementary account of the results of the cascade theory of the 
soft component of cosmic radiation, including a discussion on the lateral spread of the 
showers. 

The discussion of the physical ideas is restricted to brief accounts of the radiation 
emitted by an accelerated (classical) electron, Einstein’s explanation of the photoelectric 
effect, Compton’s theory of scattering, and the Dirac ‘ hole’ theory. The mathematical 
development of these ideas is avoided, and there is no attempt to introduce the concepts 
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of a quantized field theory. It is perhaps a pity that the idea of a collision parameter is 
not used in discussing the effect of atomic shielding in the bremsstrahlung process. 

The variation of the cross sections with energy and angle is fully illustrated, and there 
are useful tables of data. The book should be of considerable value in presenting the 
results of quantum-electrodynamics to those who are not specialists in theory, or in cosmic- 
ray work. 

It is regrettable that in a monograph which is intended for a wide public many errors 
in the formulae have escaped the proof reader. In some instances, also, economy of words 
has been carried too far. For example, on page 95 the Rutherford scattering formula is 
given, without any reference to the fact that the scattering potential is of the Coulomb 
type. The use of the word ‘ gaussianly ’ (page 97) could have been avoided. 

J. HAMILTON. 


Faraday’s Discovery of Electromagnetic Induction, by THOMAS MartTIN. 
Pp. 160. (London: Edward Arnold & Co., 1949.) 9s. 


As is well known, Faraday kept a detailed diary covering his years of mature research 
at the Royal Institution, that is, from 1820 to 1862. This was first published at intervals 
between 1932 and 1936, edited by Thomas Martin from the manuscript in the possession 
of the Royal Institution. So voluminous a work—it is in seven volumes—is, however, 
accessible to comparatively few, so that all who are interested in Faraday’s work and in 
the early history of electromagnetism will welcome the appearance of Mr. Martin’s new 
book. 

Faraday’s published papers and the Experimental Researches in Electricity give in a 
masterly manner the results of his researches ; but the steps by which he was led to his 
epoch-making discoveries, the little revealing touches that show how his mind was 
working, these are to be found only in the Diary itself. Mr. Martin, by a judicious choice 
of extracts, has succeeded in conveying much of the charm and intimacy of the Diary. 
His book is not, however, merely an abridged version of the relevant portions of the larger 
work, but gives a connected account of the early history of electromagnetism from the 
announcement in 1820 of Oersted’s discovery to the conclusion of Faraday’s electromagnetic 
researches in the spring of 1832, the period covered by the first volume of the Diary in 
its printed form. 

A brief survey of the main developments in electricity arising out of Galvani’s discovery 
provides the historical setting for Faraday’s work. This is followed by an account of 
Faraday’s own early experiments in electricity—those on electromagnetic rotations, made 
in the latter part of 1821. After the successful conclusion of these experiments ten years 
elapsed before Faraday was again free to give his full attention to electricity, though 
scattered references in the Diary show how his thoughts dwelt on the subject. Then, 
between the end of August 1831 and the beginning of March, 1832, came that amazing 
series of experiments and observations ‘ which may be said to constitute the discovery 
of electromagnetic induction’. Mr. Martin, in describing these researches, makes 
frequent use of Faraday’s own words ; the diagrams. are, with few exceptions, repro- 
ductions of Faraday’s freehand sketches. At the same time much is to be found here that 
is not to be found in the Diary itself. For instance, we read (p. 59) that the art of covering 
copper wire with silk or cotton was apparently a monopoly in the hands of the opticians, 
who charged exorbitant prices ; that experimenters were advised to take their own copper 
wire to a bonnet-wire coverer living in Spitalfields, who worked at much lower rates; 
iron wire could be obtained cotton-covered and was probably bonnet-wire used for 
stiffening ladies’ bonnets. 

We have nothing but praise for this modest little volume, and hope that Mr. Martin will 
deal in a similar manner with other aspects of Faraday’s work. 

N. H. DE V. HEATHCOTE. 
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ABSTRACTS COR SEC BONS 


The Role of Convergence in Stereoscopic Vision, by W. D. Wricut. 


ABSTRACT. An apparatus for testing stereoscopic acuity has been constructed in which 
the distance of a reference aperture P is compared with that of a test aperture Q selected 
from one of a tapering series of aperture plates. The apertures are seen as illuminated discs 
in a black field and all monocular aids to depth perception are eliminated. ‘The observer has 
to state whether Q is in front or behind P, and the acuity is derived by plotting the percentage 
estimates against the distance of Q. Obsenvarane have been made for condition A, in which 
the observer is allowed to look first at P and then at Q, and for condition B, in which the 
observer fixates on P and views Q extra-foveally. From a comparison of the acuity for 
the two conditions for various angular separations 9 between P and Q, deductions can be 
made about the relative contributions to stereoscopic vision of convergence and retinal 
disparity. The results suggest that convergence makes a significant contribution and 
has become the predominant factor at 9=20°. For 0=14°, the acuity for condition A is 
due entirely to convergence, since the blind spot interferes with the retinal disparity. Tests 
with flash exposures suggest that small eye movements do not play any important part in 
depth perception. ‘The presence of depth perception under conditions of gross diplopia was 
confirmed, indicating that fusion is not an essential item for three-dimensional vision. 
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Penetration by High- Velocity (‘ Munroe’) Jets: I, by D.C. Pack and W. M. Evans. 


ABSTRACT. By means of certain simplifying assumptions a formula is developed for 
the penetration into a ductile target by a high-velocity (‘ Munroe ’) jet. The action of the 
jet is divided into two stages, each making its contribution to the total penetration. In the 
first stage a hole is formed by the lateral compression of the target as the jet penetrates it; 
the second stage begins when the last particle of the jet has ceased to act, the hole continuing 
to deepen until the residual energy in the target has been spent. At the high pressures 
set up by a Munroe jet the strength of the target plays only a subsidiary part in the 
phenomenon. 


Penetration by High-Velocity (‘Munroe’) Jets: II, by W. M. Evans and D. C. Pack. 


ABSTRACT. Work is described in which the two stages of penetration by a high-velocity 
(‘Munroe’) jet were separated experimentally. The large penetrations measured in lead 
targets are shown to result from the flow which takes place in the metal after the jet itself 
has been consumed. For a given jet at a given stand-off it is possible to predict the 
penetration into a combination of targets from the results of a very small number of standard 
experiments. ‘The method depends upon the calculation of a certain quantity which is 
constant for a given jet at a given stand-off, and examples are given of the determination of 
this quantity from experimental data. 


An Electromechanical Source of Elastic Waves in the Ground, by F. F. Evison. 


ABSTRACT. In order to extend the study and use of the elastic properties of geological 
bodies to include effects involving the amplitude and frequency of elastic radiation, a new 
source of energy is needed. 

A large moving-coil transducer is described which, when attached rigidly to the rock, 
radiates sinusoidal vibrations of controllable amplitude, frequency and duration. The 
operation of the vibrator is analysed in detail by the method of electromechanical analogy. 

For experiments in a chalk-mine, energy at each of the frequencies 300, 600 and 1,000 c/s. 
was radiated in the form of a square-enveloped pulse, of 20 msec. duration. Vibrations were 
detected in the rock by means of a piezoelectric microphone, amplifiers and an oscillograph, 
with an overall magnification of about 109. Oscillograms show the complex disturbances 
produced by the arrival of direct and reflected energy, at distances out to the limit of the mine 
at 450 ft. Measurements of amplitude at 130 ft. are used to examine the mechanism of the 
transfer of energy from the vibrator to the rock. 

The power radiated was about 0:05 watt at 600 c/s., of which about one-sixth was in the 
compressional mode and five-sixths in the shear mode; but this power may readily be 
increased. 


Contour Modes of Square Plates excited Piezoelectrically and Determination of 
Elastic and Piezoelectric Coefficients, by R. BECHMANN. 


ABSTRACT. The suitability of square plates, vibrating in contour modes, for the 
determination of piezoelectric coefficients is considered. Piezoelectric measurements on 
such plates provide checks for the theoretical solution of the various modes and in particular 
for the distribution of displacement. The recently published solutions for square plates 
form a basis for this work and the reliability of these solutions is discussed. "The materials 
used were sodium chlorate, sodium bromate, quartz, ammonium dihydrogen phosphate and 
ethylene diamine tartrate (EDT). New determinations of the elastic and piezoelectric 
coefficients of these materials have been made, with the exception of EDT, the data for 
which were published in a previous paper. 


An Absolute Measurement of the Susceptibility of Tantalum and other Metals, by 
F. E. Hoare and J. C. WALLING. 
ABSTRACT. ‘The magnetic susceptibility of tantalum has been determined absolutely, 


the value obtained for the mass susceptibility being yr, x 10°=0-8490+0-0006 at 20° c. 
Small specimens of tantalum have been used to determine, by a comparative method, the 
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susceptibilities of specimens of platinum, palladium and rhodium at 20°c. ‘The values obtained 
are: yptX 10°=0-9712+0-0007, xpqg x 10°=5-231+0-004, yRn x 10°=0-9903+0-0008. All 
specimens were prepared from spectrographically standardized materials. 


Thermal Instability of Contact Rectifiers: The Effect of the Constituent Materials 
on the Efficiency of a Rectifying Junction, by E. BILuic. 


ABSTRACT. If due allowance is made for the lowering by the image force of the height of 
the potential barrier in contact rectifiers, current—voltage characteristics can be calculated 
which are in good agreement with observed characteristics. On application of a voltage in 
the blocking direction, electrons (or holes) pass over the barrier owing to their thermal 
energy, the number of successful passes, i.e. the leakage current, being exponentially 
dependent on temperature. As the applied voltage is increased and the barrier lowered, the 
power loss corresponding to ‘the increased leakage current will eventually raise the temper- 
ature in the thin barrier layer sufficiently to increase the current in turn until the cooling 
arrangement provided cannot cope with the increased power loss and the rectifier becomes 
thermally unstable. The conditions for this instability to arise are investigated quantita- 
tively. Amongst other parameters such as mode of operation, cooling, etc., the conductivity 
of the semiconductor, and in particular the height of the potential barrier, have a paramount 
effect on the maximum voltage that can be withstood. The simultaneous requirements for 
large resistivity and barrier height severely restrict the choice of materials which can be 
used for efficient rectifiers. ’ 


The Role of the Self Magnetic Field in High Current Gas Discharges, by P. C. 
THONEMANN and W. T. Cowuic. 


ABSTRACT. It is shown that the self magnetic field of a high current discharge plays 
an important part in determining the radial electron density distribution. The effect 
of this magnetic field can be treated by adding a ‘magnetic potential’ term, Vm=f(wH dr 
to the electrostatic potential V; in the Boltzmann equation for the radial electron density 
distribution, w being the longitudinal electron drift velocity. The drift of electrons and 
ions to the tube walls is equalized at low are currents by a radial electric field. Since the 
magnetic field reduces the drift of electrons to the wall whilst having little influence upon 
the ions, increasing the arc current leads to a progressive reduction of the radial electric 
field. Probe measurements in an arc at low pressures show that at currents over 
50 amperes, the radial electric field is zero except for a narrow sheath at the wall. The 
axial density is then mainly determined by the self magnetic field. A theoretical 
expression is derived for the equilibrium density distribution assuming a balance between 
the diffusion rates due to the concentration gradient and the magnetic forces. Experiments 
in a mercury vapour arc at high currents support this theory. 


The Relationship between the Radial Concentration Variations of excited Atoms 
and Electrons in a Discharge in Thermal Equilibrium, by H. EpDEts. 


ABSTRACT. The excitation and Saha equations applicable to gaseous electrical 
discharges in thermal equilibrium are used to obtain convenient expressions for calculating 
the radial concentration variations of electrons and excited atoms when the radial intensity 


variation of one spectral line is known, and when the temperature on the axis of the discharge 
can be estimated. 
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Figure 5. Cu foil, nearly normal transmission. Figureé Annealed rolled Cu; reflection 
after etching. 


Figures 50, 51,52. 5.6 Cyclo-penteno-1.2 Benzanthracene. 


Figures 53, 54, 55. 1.2 Benzpyrene. 
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Figure 15. 
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Figure 16, CdI, quintuplet, Figure 17. Cdl, tailed arcs. 
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Figure 20. Pd leaf (arced doublet). 


Figure 23. Al leaf, recrystallized. 


Figure 25. NaNOs, ‘ tailed” arcs, 
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Figure 26 a. Cs3.H¢, doublet. Figure 27. Paraffin wax. Figure 28. C,,H4. (+ Au rings). 


Figures 29, 30, 31. 


Figures 32, 33. n-C,,.H.,0H. Figure 34. n-C,;HsgsI. 


Figures 35, 36, 37. n-C,;H3;,OH. 
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Figure 41. Elaidic acid. Figure 42. Fet** stearate. Figure 4 


Figure 44. Anthracene. Figure 45a. 1.2 
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Figure 47.. 5Me.1.2 B.A. 
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@)miG@nitials 25.) (a) “(Plan x3) 


(b) (Initial 2-5.) (b’) (Plan <3.) (d’) (Side view x3.) 


Figure 56. Rotational slip in potassium ferrocyanide trihydrate; four crystals, a—d. 
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Figure 57. Rotational slip in gypsum (plan views of three crystals). 
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